matlaB作业
matlab作业
上机实践练习五实验五：MATLAB 数据可视化实验目的：掌握MATLAB 二维、三维图形绘制，掌握图形属性的设置和图形修饰；掌握图像文件的读取和显示。
学时：2学时实验内容：(1)二维图形绘制。
(2)隐函数绘图(2) 三维曲线和三维曲面绘制。
(3) 图像文件的读取和显示。
1.写出实现下列图形的M文件。
t=[0:0.1:pi];y1=sin(2*t);y2=cos(3*t);plot(t,y1,'b-o');hold onplot(t,y2,'k:o');xlabel('ʱ¼ä')ylabel('·ùÖµ')text(1,sin(2*1),'\fontsize{16}\leftarrowsin(2t) ');text(2.4,cos(3*(2.4)),'\fontsize{16}cos(3t)\rightarrow ',... 'HorizontalAlignment','right')2. 将窗口分割成4个区域，并且在[]0,2π区间上绘制1sin 2,6y x π⎛⎫=+⎪⎝⎭()224cos 3,y x =()3cos y x =，()35sin 1y x =+四条曲线，并且给每一个图形添加标题和标注。
clearx=0:pi/10:2*pi; y1=sin(2.*x+pi/6); y2=4.*cos(3.*x.^2); y3=cos(x); y4=5.*sin(x+1); subplot(2,2,1); plot(x,y1); title('y1');text(1,sin(2*1+pi/6),'sin(2*x+pi/6)') subplot(2,2,2); plot(x,y2)title('y2');text(2,4.*cos(3*2^2),'4*cos(3*x^2)') subplot(2,2,3); plot(x,y3)title('y3');text(2,cos(2),'cos(x)')subplot(2,2,4);plot(x,y4)title('y4');text(2,5.*sin(3),'5*sin(x+1)')3.已知y1=x2，y2=cos(2x)，y3＝y1*y2，其中x为取值-2π~2π的等差数列(每次增加0.02π)，完成下列操作(1)在同一坐标系下用不同的颜色和线型绘制三条曲线，给三条曲线添加图例。
MATLAB课程作业
注意：11月15日是交作业的最后期限。
（交word和MATLAB程序清单，word以“班级-学号-姓名”命名，每个MATLAB程序以”zuoye_题号”命名后，放在一个文件夹下（文件夹命名为“作业程序”），如zuoye_1_1.m），将word和“作业程序”文件夹放在一个文件夹（文件夹以“班级-学号-姓名”）交到班长处，班长统一交到教学办公室。
第一部分程序设计1.1、用MA TLAB可以识别的格式输入下面两个矩阵：(1) 矩阵A的维数；(2) 矩阵A中的元素a41的值；(3) 修改矩阵A的元素，使a41 =3.0；(4) 矩阵A中最后2行和最后3列交汇形成的子矩阵的值。
(5)求出A和B的乘积矩阵C，并将C矩阵的右下角2X 3子矩阵赋给D矩阵。
1.2、已知111121111,131,111214A B⎡⎤⎡⎤⎢⎥⎢⎥=-=-⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求(1)AB-2A，(2)A*B，(3)A﹒*B，(4) AB-BA 1.3、解线性方程：1.4、解方程组：1.5、用MA TLAB语言实现下面的分段函数：1.6、已知x=[1 2 3 ],y=[4 5 6 ]，试计算z=x.*y 、x.\y 和x./y.1.7、分别用for 和while 循环语句编写程序，求出1.8、 已知在平面坐标中两点(x 1, y 1)和(x 2, y 2)之间的距离计算公式为 ()()222121y y x x L -+-=(1) 利用命令文件的形式，编写求解该距离的M 文件dis1.m ；(2) 利用函数文件的形式，编写求解该距离的M 文件dis2.m ；(3) 给定两点坐标的值(2,3)和(8, －5)，试分别调用命令文件dis1.m 和函数文件dis2.m 求解该两点间距离的值。
1.9、求解方程x 5+6x 3一3x 2＝10的5个根，并将其位置用五角星符号标记在复平面上，要求横纵坐标袖的刻度等长，注明虚轴和实轴，在title 位置上写㈩方程。
matlab作业
实验三某厂利用a、b、c三种原料生产A、B、C三种产品，已知生产每种产品在消耗原料方面的各项技术条件和单位产品的利润，以及可利用的各种原料的量，是制定适当的生产规划使得该厂的总利润最大。
对于这个问题，可先将该问题转化为数学模型。
Z=2x1+4x2+3x3X1+3x2+2x3<=802x1+x2+2x3<=403x1+4x2+2x3<=60X1,x2,x3>=0据此可以画出图形，可依靠图形来求出Z的最大值。
c=[-2 -4 -3]A=[1 3 2;2 1 2;3 4 2]b=[60 40 80]lb=[0 0 0][x,f,exitflag,output]=linprog(c,A,b,[],[],lb,[])将上述程序输入matlab。
就可得出结果。
结果如下：x =4.000012.000010.0000f =-86.0000exitflag =1output =iterations: 7cgiterations: 0algorithm: 'lipsol'由上述结果可以知道当A、B、C分别生产4,12,10 公斤时，总利润最大，可达到86。
某工厂生产甲、乙、丙三种产品，单位产品所需工时分别为2、3、1个；单位产品所需原料分别为3、1、5公斤；单位产品利润分别为2、3、5元。
工厂每天可利用的工时为12个，可供应的原料为15公斤。
为使总利润为最大，试确定日生产计划和最大利润。
对于这个问题，可先将该问题转化为数学模型。
Z=2x+3y+5z2x+3y+z<=123x+y+5z<=15据此可以画出图形，可依靠图形来求出Z的最大值。
c=[-2 -3 -5]A=[2 3 1;3 1 5]b=[12;15]lb=[0 0 0][x,Z,exitflag,output]=linprog(c,A,b,[],[],lb,[])将上述程序输入matlab。
就可得出结果。
Matlab大作业
Matlab的基本操作一、使用函数实现对下列矩阵的左旋和右旋以及反转已知答案：如图：矩阵如图矩阵的左旋如图一矩阵的右旋如图二矩阵的左右反转如图三矩阵的上下反转如图四二、已知A=[8 9 5 ] B=[-1 3 -2][36 -7 11] [2 0 3][21 -8 5] [-3 1 9 ]计算：1)A+5*B;2)A*B和A.*B3)A^3和A.^34)A/B和B\A5)[A,B]答案：A=[8 9 5；36 -7 11；21 -8 5];B=[-1 3 -2；2 0 3；-3 1 9];A+5*BA*BA.*BA^3A.^3A/BA\B[A,B]程序运行结果如下：第二章、Matlab程序设计一、已知S=1+2+22+23+…+263求S的值答案：代码如下：ClearClcs=0,j=2for i=1:63s=s+j^iends程序运行结果如下：二、第三章、Matlab绘图一、在同一坐标系中画出下列函数的图像：x2，-x2，xsin（x）在[0,2π]上的函数图像程序代码如下：clearclcx=0:pi/100:2*pi;y1=x.^2;y2=-x.^2;y3=x.*sin(x);title('同一坐标下的函数图像')plot(x,y1,':',x,y2,'h',x,y3,'--')程序运行结果如下：二、绘制极坐标图像：程序代码如下：clearclcx=0:0.01:2*pi;y=2*cos(2*(x-pi/8));title('极坐标图像')xlabel('x')ylabel('y')polar(x,y)程序运行结果如下：第四章、Matlab符号运算一、求函数y=e-x2 的傅立叶变换及其逆变换程序代码如下：clearclcsyms x t;y=exp(-x^2);Ft=fourier(y,x,t)fx=ifourier(Ft,t,x)程序运行结果如下：二、求下列极限值：答案：1、程序代码：clearclcsyms x;s1=sin(2*x)/sin(5*x);limit(s1,x,0)s2=(1+1/x)^(2*x);limit(s2,x,inf)程序运行结果如下：第五章、Matlab数值运算一、建立一个3*4阶随机矩阵，求矩阵的最大值、最小值、方差和标准差、极差、协方差，和自相关阵程序代码如下：clearclcA=rand(3,4)B=var(A)C=std(A)D=range(A)E=cov(A)F=corrcoef(A)程序运行结果如下：二、求函数f（x）=x3-2x+1在x=[-1 1]之间的极小值和x=-1附近的零点程序代码如下：clearclc[x,y]=fminbnd('x.^3-2.*x+1',-1,1)[x,y]=fzero('x.^3-2.*x+1',-1)程序运行结果如下：第六章、Matlab图形用户界面一、建立一个具有三个输入框的窗口对话框程序代码如下：clearclcprompt={'姓名','年龄','班级'}；DTitle='注册学生信息';line[1;1;1];def{'乔阳'，'18','100412105'};info=inputdlg(prompt,DTitle,line,def,'ON')程序运行结果如下：二、设计一个表现下载进度的进度条程序代码如下：clearclch = waitbar(0,'正在下载，请稍等……');for i=1:1000waitbar(i/10000,h)endclose(h)程序运行结果如下：第七章、Matlab Simulink 仿真一、仿真信号x（t）=sin(t)sin(10t)的波形仿真步骤：1.建立模型窗口：生成无标题（intitled）的模型窗口；2.添加信号源模块（sine）、输出模块（scope）、数学模块（Dot product）3.设置模块参数：如下图：4.编辑模块即将各个模块连接起来。
matlab1-8章课后作业
MATLAB基础教程1~8章作业Matlab第一章1.阐述Matlab的功能Matlab作为一种高级计算软件，是进行算法开发、数据可视化、数据分析以及数值计算的交互式应用开发环境，已被广泛应用于不同领域。
Matlab的基本功能包括：数学计算功能、图形化显示功能、M语言编程功能、编译功能、图形用户界面开发功能、Simulink建模仿真功能、自动代码生成功能。
Matlab第二章1．创建double的变量，并进行计算。
(1)a=87，b=190，计算a+b、a-b、a*b。
(2)创建uint8 类型的变量，数值与(1)中相同，进行相同的计算。
>> a=87,b=190a =87b =190>> a+bans =277>> a-bans =-103>> a*bans =16530>> c=uint8(87), d=uint8(190)c =87d =190>> c+dans =255>> c-dans =>> c*dans =2552.计算(1)sin(60)(2)e^3(3)cos(3π/4)>> sind(60)ans =0.8660>> exp(3)ans =20.0855>> cos(3*pi/4)ans =-0.70713.设u=2，v=3，计算：(1)(2)(3)>> u=2;>> v=3;>> 4*u*v/log(v)ans =21.8457>> (exp(u)+v)^2/(v^2-u)ans =15.4189>> sqrt(u-3*v)/(u*v)ans =0 + 0.4410i4.计算如下表达式：(1)(2)>> (3-5*i)*(4+2*i)ans =22.0000 -14.0000i>> sin(2-8*i)ans =1.3553e+003 +6.2026e+002i5.判断下面语句的运算结果。
MATLAB作业
MATLAB作业⼀、必答题：1. MATLAB系统由那些部分组成？答：MATLAB系统主要由开发环境、MATLAB语⾔、MATLAB数学函数库、图形功能和应⽤程序接⼝五个部分组成。
2. 如何启动M⽂件编辑/调试器？答：在操作界⾯上选择“建⽴新⽂件”或“打开⽂件”操作时，M⽂件编辑/调试器将被启动。
在命令窗⼝中键⼊“edit”命令也可以启动M⽂件编辑/调试器。
3. 存储在⼯作空间中的数组能编辑吗？如何操作？答：存储在⼯作空间的数组可以通过数组编辑器进⾏编辑：在⼯作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输⼊修改内容即可。
4. 在MATLAB中有⼏种获得帮助的途径？答：在MATLAB中有多种获得帮助的途径：（1）帮助浏览器：选择view菜单中的Help菜单项或选择Help菜单中的MATLAB Help菜单项可以打开帮助浏览器；（2）help命令：在命令窗⼝键⼊“help” 命令可以列出帮助主题，键⼊“help 函数名”可以得到指定函数的在线帮助信息；（3）lookfor命令：在命令窗⼝键⼊“lookfor 关键词”可以搜索出⼀系列与给定关键词相关的命令和函数（4）模糊查询：输⼊命令的前⼏个字母，然后按Tab键，就可以列出所有以这⼏个字母开始的命令和函数。
5. 有⼏种建⽴矩阵的⽅法？各有什么优点？答：（1）以直接列出元素的形式输⼊；（2）通过语句和函数产⽣；（3）.在m⽂件中创建矩阵；（4）从外部的数据⽂件中装⼊。
6. 命令⽂件与函数⽂件的主要区别是什么？答：命令⽂件: M⽂件中最简单的⼀种，不需输出输⼊参数，⽤M ⽂件可以控制⼯作空间的所有数据。
运⾏过程中产⽣的变量都是全局变量。
运⾏⼀个命令⽂件等价于从命令窗⼝中顺序运⾏⽂件⾥的命令，程序不需要预先定义，只要依次将命令编辑在命令⽂件中即可。
函数⽂件：如果M⽂件的第⼀个可执⾏⾏以function开始，便是函数⽂件，每⼀个函数⽂件定义⼀个函数。
matlab大作业参考课题
matlab大作业参考课题
当选择一个 MATLAB 大作业的参考课题时，你可以考虑以下几
个方面来确定一个合适的课题：
1. 兴趣和专业方向，首先，你可以考虑自己的兴趣和专业方向。
如果你对信号处理感兴趣，可以选择音频处理或图像处理的课题；
如果你对控制系统感兴趣，可以选择系统建模和控制设计的课题。
2. 实际应用，考虑选择与实际应用相关的课题，这样可以增加
课题的实用性和吸引力。
例如，可以选择医学图像处理、金融数据
分析或者工程控制系统设计等方面的课题。
3. 数据和资源，确保你能够获取到相关的数据和资源来支持你
的课题研究。
这些资源可以是公开数据集、实验数据或者模拟数据。
4. 研究难度，考虑课题的研究难度，确保它既不会过于简单以
至于缺乏挑战性，也不会过于复杂以至于无法完成。
举例来说，如果你对音频处理感兴趣，可以选择一个课题是基
于 MATLAB 的音频信号处理和音乐合成；如果你对图像处理感兴趣，
可以选择一个课题是基于 MATLAB 的图像处理和分析；如果你对控制系统感兴趣，可以选择一个课题是基于 MATLAB 的系统建模和控制设计。
总之，选择一个适合的 MATLAB 大作业参考课题需要考虑到自己的兴趣、专业方向、实际应用、资源获取和研究难度等因素，希望这些建议对你有所帮助。
MATLB实验作业
实验一MATLAB运算基础1、用逻辑表达式求下列分段函数的值。
t=0:0.5:2.5y=t.^2.*((t>=0)&(t<1))+(t.^2-1).*(( t>=1)&(t<2))+(t.^2-2*t+1).*((t>=2 )&(t<3))2、求[100，999]之间能被21整除的数的个数。
p=rem([100:999],21)==0;sum(p)3、建立一个字符串向量，删除其中的大写字母。
ch='KdDdfKaWdsfCI',k=find(ch>=' A'&ch<='Z'),ch(k)=[]4、输入矩阵，并找出A中大于或等于5的元素。
A=[1 2 3;4 5 6;7 8 9],[m,n]=find(A>=5),for j=1:length(m)x(j)=A(m(j),n(j))xend5、求矩阵的行列式值、逆和特征根。
a11=input('a11='),a12=input('a12 ='),a21=input('a21='),a22=input('a22 ='),A=[a11,a12;a21,a22],DA=det(A),IA=inv(A),EA=eig(A) 6、不采用循环的形式求出和式的数值解。
sum(2.^[0:63])实验二1、1行100列的Fibonacc数组a，a(1)=a(2)=1，a(i)=a(i-1)+a(i-2)，用for循环指令来寻求该数组中第一个大于10000的元素，并指出其位置i。
n=100;a=ones(1,n);for i=3:na(i)=a(i-1)+a(i-2);if a(i)>10000a(i),break;end;end,i2、编写M脚本文件，定义下列分段函数，并分别求出当（）、（）和（）时的函数值。
Matlab基础及应用第一章上机作业
1、新建一个文件夹（自己的名字命名）
答：
2、启动MATLAB，将该文件夹添加到MATLAB路径管理器中。
保存，关闭对话框。
使用path命令查看MATLAB搜索路径。
3、学习使用help命令，例如在命令窗口输入help eye，然后根据帮助说明，学习使用指令eye。
4、使用clc、clear，观察command window、command history和workspace等窗口的变化结果。
5、编写一段简短的脚本文件，保存并显示运行结果。
6、创建两个double型变量a=32、b=5，并计算a+b、a-b、a*b、a/b、a\b；
创建两个int8型变量a、b，取同样数值，并计算a+b、a-b、a*b、a/b、a\b，对于计算结果与前次计算结果不同的情形请给出解释。
7、查看int16数据类型的取值范围(intmin,intmax)；查看单精度数据类型的取值范围和精度(realmin,realmax,eps)。
8、求下列表达式的值，然后显示MATLAB工作空间的使用情况并保存全部变量：
(1)sin(60o) (2) e3 (3) cos(3/4π)
9、用两种定义复数的方法计算(直接定义，complex(a,b))
(1)(3-5i)(4+2i) (2) sin(2-8i)。
MATLAB基础训练作业
实验一 MATLAB 工作环境熟悉及简单命令的执行一、实验目的：熟悉MATLAB 的工作环境，学会使用MATLAB 进行一些简单的运算。
二、实验内容：MATLAB 的启动和退出，熟悉MATLAB 的桌面（Desktop ），包括菜单（Menu ）、工具条 （Toolbar ）、命令窗口(Command Window)、历史命令窗口、工作空间(Workspace)等；完成一些基本的矩阵操作；学习使用在线帮助系统。
三、实验步骤：1、启动MATLAB ，熟悉MATLAB 的桌面。
2、在命令窗口执行命令完成以下运算，观察workspace 的变化，记录运算结果。
（1）（365-52?2-70）?3（2）>>area=pi*^2（3）已知x=3，y=4，在MATLAB 中求z ：()232y x y x z -= （4）将下面的矩阵赋值给变量m1,在workspace 中察看m1在内存中占用的字节数。
m1=⎥⎥⎥⎥⎦⎤⎢⎢⎢⎢⎣⎡11514412679810115133216 执行以下命令>>m1( 2 , 3 )>>m1( 11 )>>m1( : , 3 )>>m1( 2 : 3 , 1 : 3 )>>m1( 1 ,4 ) + m1( 2 ,3 ) + m1( 3 ,2 ) + m1( 4 ,1)（5）执行命令>>help abs查看函数abs 的用法及用途，计算abs( 3 + 4i )（6）执行命令>>x=0::6*pi;>>y=5*sin(x);>>plot(x,y)（6）运行MATLAB 的演示程序，>>demo ，以便对MATLAB 有一个总体了解。
四、思考题1、以下变量名是否合法为什么（1）x2（2）3col（3）_row（4）for2、求以下变量的值，并在MATLAB 中验证。
matlab课程设计大作业
matlab课程设计大作业一、教学目标本课程的教学目标是使学生掌握MATLAB基本语法、编程技巧以及MATLAB 在工程计算和数据分析中的应用。
通过本课程的学习，学生将能够熟练使用MATLAB进行简单数学计算、线性方程组求解、函数图像绘制等。
1.掌握MATLAB基本语法和编程结构。
2.了解MATLAB在工程计算和数据分析中的应用。
3.熟悉MATLAB的函数库和工具箱。
4.能够使用MATLAB进行简单数学计算。
5.能够使用MATLAB求解线性方程组。
6.能够使用MATLAB绘制函数图像。
7.能够利用MATLAB进行数据分析和处理。
情感态度价值观目标：1.培养学生对计算机辅助设计的兴趣和认识。
2.培养学生团队合作和自主学习的能力。
二、教学内容本课程的教学内容主要包括MATLAB基本语法、编程技巧以及MATLAB在工程计算和数据分析中的应用。
1.MATLAB基本语法：介绍MATLAB的工作环境、基本数据类型、运算符、编程结构等。
2.MATLAB编程技巧：讲解MATLAB的函数调用、脚本编写、函数文件编写等编程技巧。
3.MATLAB在工程计算中的应用：介绍MATLAB在数值计算、线性方程组求解、图像处理等方面的应用。
4.MATLAB在数据分析中的应用：讲解MATLAB在数据采集、数据分析、数据可视化等方面的应用。
三、教学方法本课程采用讲授法、案例分析法、实验法等多种教学方法相结合的方式进行教学。
1.讲授法：通过讲解MATLAB的基本语法、编程技巧以及应用案例，使学生掌握MATLAB的基本知识和技能。
2.案例分析法：通过分析实际工程案例，使学生了解MATLAB在工程计算和数据分析中的应用。
3.实验法：安排上机实验，使学生在实际操作中巩固所学知识，提高实际编程能力。
四、教学资源本课程的教学资源包括教材、实验设备、多媒体资料等。
1.教材：选用《MATLAB教程》作为主要教材，辅助以相关参考书籍。
2.实验设备：为学生提供计算机实验室，配备有MATLAB软件的计算机。
Matlab习题及答案
现代计算方法Matlab 作业答案1.绘出函数f(x)=sin x x ，在[0,4]上的图形解：在M 文件输入：x=0:pi/100:4;y=x.*sin(x);plot(y)运行2. 求3x +2x +5 = 0的根解：在命令窗口输入：>> solve('x^3+2*x+5=0')ans =((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3) - 2/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3))1/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)) - ((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)/2 -(3^(1/2)*i*(2/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)) + ((108^(1/2)*707^(1/2))/108 -5/2)^(1/3)))/21/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)) - ((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)/2 +(3^(1/2)*i*(2/(3*((108^(1/2)*707^(1/2))/108 - 5/2)^(1/3)) + ((108^(1/2)*707^(1/2))/108 -5/2)^(1/3)))/23.321436min x x x z ++=120..321=++x x x t s301≥x5002≤≤x203≥x解：运用单纯形法计算此题，首先把约束条件化成标准形式：,,,,,205030120654321635241321≥=-=+=-=++x x x x x x x x x x x x x x x（1）在M 文件输入SimpleMthd 函数：function [x,minf] = SimpleMthd(A,c,b,baseVector)sz = size(A);nVia = sz(2);n = sz(1);xx = 1:nVia;nobase = zeros(1,1);m = 1;for i=1:nViaif (isempty(find(baseVector == xx(i),1)))nobase(m) = i;m = m + 1;else;endendbCon = 1;M = 0;while bConnB = A(:,nobase);ncb = c(nobase);B = A(:,baseVector);cb = c(baseVector);xb = inv(B)*b;f = cb*xb;w = cb*inv(B);for i=1:length(nobase)sigma(i) = w*nB(:,i)-ncb(i);end[maxs,ind] = max(sigma);if maxs <= 0minf = cb*xb;vr = find(c~=0 ,1,'last');for l=1:vrele = find(baseVector == l,1);if (isempty(ele))x(l) = 0;elsex(l)=xb(ele);endendbCon = 0;elsey = inv(B)*A(:,nobase(ind));if y <= 0disp('不存在最优解!');x = NaN;minf = NaN;return;elseminb = inf;chagB = 0;for j=1:length(y)if y(j)>0bz = xb(j)/y(j);if bz<minbminb = bz;chagB = j;endendendtmp = baseVector(chagB);baseVector(chagB) = nobase(ind);nobase(ind) = tmp;endendM = M + 1;if (M == 1000000)disp('找不到最优解!');x = NaN;minf = NaN;return;endend（2）在命令窗口输入：clear allA=[1 1 1 0 0 0;1 0 0 -1 0 0;0 1 0 0 1 0;0 0 1 0 0 -1];c=[6 3 4 0 0 0];b=[120;30;50;20];[xm,mf]=SimpleMthd(A,c,b,[3 4 5 6])xm =0 50 70mf =4304.计算下面函数在区间(0，1)内的最小值。
matlab1-8章课后作业
MATLAB基础教程1~8章作业Matlab第一章1.阐述Matlab的功能Matlab作为一种高级计算软件，是进行算法开发、数据可视化、数据分析以及数值计算的交互式应用开发环境，已被广泛应用于不同领域。
Matlab的基本功能包括：数学计算功能、图形化显示功能、M语言编程功能、编译功能、图形用户界面开发功能、Simulink建模仿真功能、自动代码生成功能。
Matlab第二章1．创建double的变量，并进行计算。
(1)a=87，b=190，计算a+b、a-b、a*b。
(2)创建uint8 类型的变量，数值与(1)中相同，进行相同的计算。
>> a=87,b=190a =87b =190>> a+bans =277>> a-bans =-103>> a*bans =16530>> c=uint8(87), d=uint8(190)c =87d =190>> c+dans =255>> c-dans =ans =2552.计算(1)sin(60)(2)e^3(3)cos(3π/4)>> sind(60)ans =0.8660>> exp(3)ans =20.0855>> cos(3*pi/4)ans =-0.70713.设u=2，v=3，计算：(1)(2)(3)>> u=2;>> v=3;>> 4*u*v/log(v)ans =21.8457>> (exp(u)+v)^2/(v^2-u) ans =15.4189>> sqrt(u-3*v)/(u*v) ans =0 + 0.4410i 4.计算如下表达式：(1)(2)>> (3-5*i)*(4+2*i)22.0000 -14.0000i>> sin(2-8*i)ans =1.3553e+003 +6.2026e+002i5.判断下面语句的运算结果。
matlab综合大作业(附详细答案)
m a t l a b综合大作业(附详细答案)-标准化文件发布号：（9456-EUATWK-MWUB-WUNN-INNUL-DDQTY-KII《MATLAB语言及应用》期末大作业报告1．数组的创建和访问（20分，每小题2分）：1)利用randn函数生成均值为1，方差为4的5*5矩阵A；实验程序：A=1+sqrt(4)*randn(5)实验结果：A =0.1349 3.3818 0.6266 1.2279 1.5888-2.3312 3.3783 2.4516 3.1335 -1.67241.2507 0.9247 -0.1766 1.11862.42861.5754 1.6546 5.3664 0.8087 4.2471-1.2929 1.3493 0.7272 -0.6647 -0.38362)将矩阵A按列拉长得到矩阵B；实验程序：B=A(:)实验结果：B =0.1349-2.33121.25071.5754-1.29293.38183.37830.92471.65461.34930.62662.4516-0.17665.36640.72721.22793.13351.11860.8087-0.66471.5888-1.67242.42864.2471-0.38363)提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的矩阵C；实验程序：C=[A(2,2),A(2,4);A(3,2),A(3,4)]实验结果：C =3.3783 3.13350.9247 1.11864)寻找矩阵A中大于0的元素；]实验程序：G=A(find(A>0))实验结果：G =0.13491.25071.57543.38183.37830.92471.65461.34930.62662.45165.36640.72721.22793.13351.11860.80871.58882.42864.24715)求矩阵A的转置矩阵D；实验程序：D=A'实验结果：D =0.1349 -2.3312 1.2507 1.5754 -1.29293.3818 3.3783 0.9247 1.6546 1.34930.6266 2.4516 -0.1766 5.3664 0.72721.2279 3.1335 1.1186 0.8087 -0.66471.5888 -1.67242.4286 4.2471 -0.38366)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E；实验程序：E=flipud(fliplr(A))实验结果：E =-0.3836 -0.6647 0.7272 1.3493 -1.29294.2471 0.80875.3664 1.6546 1.57542.4286 1.1186 -0.1766 0.9247 1.2507-1.6724 3.1335 2.4516 3.3783 -2.33121.5888 1.2279 0.6266 3.3818 0.13497)删除矩阵A的第2列和第4列得到矩阵F；实验程序：F=A;F(:,[2,4])=[]实验结果：F =0.1349 0.6266 1.5888-2.3312 2.4516 -1.67241.2507 -0.17662.42861.5754 5.3664 4.2471-1.2929 0.7272 -0.38368)求矩阵A的特征值和特征向量；实验程序：[Av,Ad]=eig(A)实验结果：特征向量Av =-0.4777 0.1090 + 0.3829i 0.1090 - 0.3829i -0.7900 -0.2579 -0.5651 -0.5944 -0.5944 -0.3439 -0.1272-0.2862 0.2779 + 0.0196i 0.2779 - 0.0196i -0.0612 -0.5682 -0.6087 0.5042 - 0.2283i 0.5042 + 0.2283i 0.0343 0.6786 0.0080 -0.1028 + 0.3059i -0.1028 - 0.3059i 0.5026 0.3660 特征值Ad =6.0481 0 0 0 00 -0.2877 + 3.4850i 0 0 00 0 -0.2877 - 3.4850i 0 00 0 0 0.5915 00 0 0 0 -2.30249)求矩阵A的每一列的和值；实验程序：lieSUM=sum(A)实验结果：lieSUM =-0.6632 10.6888 8.9951 5.6240 6.208710)求矩阵A的每一列的平均值；实验程序：average=mean(A)实验结果：average =-0.1326 2.1378 1.7990 1.1248 1.24172．符号计算（10分，每小题5分）：1)求方程组20,0++=++=关于,y z的解；uy vz w y z w实验程序：S = solve('u*y^2 + v*z+w=0', 'y+z+w=0','y,z');y= S. y, z=S. z实验结果：y =[ -1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))-w] [ -1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))-w] z =[ 1/2/u*(-2*u*w-v+(4*u*w*v+v^2-4*u*w)^(1/2))] [ 1/2/u*(-2*u*w-v-(4*u*w*v+v^2-4*u*w)^(1/2))]2)利用dsolve 求解偏微分方程,dx dyy x dt dt==-的解； 实验程序：[x,y]=dsolve('Dx=y','Dy=-x')实验结果：x =-C1*cos(t)+C2*sin(t)y = C1*sin(t)+C2*cos(t)3．数据和函数的可视化（20分，每小题5分）：1)二维图形绘制：绘制方程2222125x y a a +=-表示的一组椭圆，其中0.5:0.5:4.5a =；实验程序：t=0:0.01*pi:2*pi; for a=0.5:0.5:4.5; x=a*cos(t); y=sqrt(25-a^2)*sin(t); plot(x,y) hold on end实验结果：2) 利用plotyy 指令在同一张图上绘制sin y x =和10x y =在[0,4]x ∈上的曲线；实验程序：x=0:0.1:4; y1=sin(x); y2=10.^x;[ax,h1,h2]=plotyy(x,y1,x,y2); set(h1,'LineStyle','.','color','r'); set(h2,'LineStyle','-','color','g'); legend([h1,h2],{'y=sinx';'y=10^x'});实验结果：3)用曲面图表示函数22z x y =+；实验程序：x=-3:0.1:3; y=-3:0.1:3; [X,Y]=meshgrid(x,y); Z=X.^2+Y.^2; surf(X,Y,Z)实验结果：4)用stem 函数绘制对函数cos 4y t π=的采样序列；实验程序：t=-8:0.1:8;y=cos(pi.*t/4); stem(y)实验结果：4. 设采样频率为Fs = 1000 Hz ，已知原始信号为)150π2sin(2)80π2sin(t t x ⨯+⨯=，由于某一原因，原始信号被白噪声污染，实际获得的信号为))((ˆt size randn x x+=，要求设计出一个FIR 滤波器恢复出原始信号。
Matlab 作业1
喻晓磊 200731510103 电子科学与技术习题一：编写M文件，计算以下问题：口口口口×口＝口口口口以上9个口代表1～9这9个数字，不得遗漏或重复要求：给出解算思路和M文件代码注意算法的效率解题理念：1，设以上方框中的数字分别是1~9之间的数字a、b、c、d、e、f、g、h、i。
2，我们可以推出，个位数e肯定不是1或者5 ，否则两个四位数的个位必然相等，即d=i；再者，其不可能等于9，由不重复的1~9组成的四位数，最大为9876，最小为1234，故e的最大可能取值为9876/1234=8.003……，即最大取8 。
3，我们看到，e的最小取值为2，那么a肯定不能大于或等于5 ，否则乘法运算后的结果为5位数，所以a 只能取1~4。
4，在此基础上，我们采取“穷举法”，逐个试验，从 a 开始，每一位与前面出现位的数字不相等，然后检测其是否满足算式。
程序代码如下：global aglobal bglobal cglobal dglobal eglobal fglobal gglobal hglobal ifor a=1:4for b=1:9if b~=a;for c=1:9if (c~=b)&&(c~=a);for d=1:9if (d~=b)&&(d~=a)&&(d~=c);for e=2:8 e~=5;if (e~=a)&&(e~=b)&&(e~=c)&&(e~=d);for f=1:9if (f~=a)&&(f~=b)&&(f~=c)&&(f~=d)&&(f~=e);for g=1:9if(g~=a)&&(g~=b)&&(g~=c)&&(g~=d)&&(g~=e)&&(g~=f);for h=1:9if(h~=a)&&(h~=b)&&(h~=c)&&(h~=d)&&(h~=e)&&(h~=f)&&(h~=g);for i=1:9if(i~=a)&&(i~=b)&&(i~=c)&&(i~=d)&&(i~=e)&&(i~=f)&&(i~=g)&&(i~=h);x=1000*a+100*b+10*c+d;y=e*x;z=1000*f+100*g+10*h+i;if y==z;r=xs=et=yend;end;******end;end;得到的结果如下：我们看到，结果有两组：1738 * 4= 69521963 * 4= 7852这两组结果都显示到了MATLAB软件的命令窗（COMMAND WINDOW）中，我们是使用了一个小技巧：将x、e、y的值赋给r、s、t时并没有在句末加分号，所以计算的过程也显示出来了，但是看工作区（work space）中，r、s、t 的值却只对应后面一组，这是因为，当循环进行到算出第二组结果时，便替换掉了第一组的值。
MATLAB第一章作业答案
第一章M A T L A B概况与基本操作1.选择题：(1)最初的MATLAB核心程序是采用A语言编写的。
(2)即将于2011年9月发布的MATLAB新版本的编号为D。
2011Ra 2011Rb R2011a R2011b(3)在默认设置中，MATLAB中的注释语句显示的颜色是D。
A.黑色B.蓝色C.红色D.绿色(4)如果要以科学计数法显示15位有效数字，使用的命令是B。
long long e long g long d(5)在命令窗口新建变量a、b，如果只查看变量a的详细信息，使用的命令为B。
a a(6)如果要清除工作空间的所有变量，使用的命令为 C 。
all C.两者都可 D.两者都不可(7)在创建变量时，如果不想立即在命令窗口中输出结果，可以在命令后加上D。
A.冒号B.逗号C.空格D.分号(8)如果要重新执行以前输入的命令，可以使用B键。
A.下箭头↓B.上箭头↑C.左箭头←D.右箭头→(9)如果要查询函数inv的相关信息，并显示在命令窗口，应使用命令A。
D.三者均可(10)如果要启动Notebook文档，下列D操作是可行的。
A.在命令窗口输入notebook命令B.在命令窗口输入notebook filename命令C.在Word中启动M-book文档D.三者均可2.填空题：(1)MATLAB是MATrix和LABoratory两个单词前三个字母的组合，意为“矩阵实验室”，它的创始人是Cleve Moler和Jack Little。
(2)在MATLAB的默认设置中，关键字显示的字体为蓝色，命令、表达式、计算结果显示的字体为黑色，字符串显示的字体为紫色，注释显示的字体为绿色，错误信息显示的字体为红色。
(3)在命令窗口中，输出结果显示为各行之间添加空行的命令为format compact，各行之间不添加空行的命令为format compact。
备注：本题布置给大家时有一点小错误，现在予以更正。
山东大学《matlab程序设计》期末作业
《Matlab 程序设计》 期末作业学院： 专业： 班级：学号： 姓名： 成绩：一 、 简答题（每题10分，共50分）1、假设有两个矩阵A 和B ，请用这两个矩阵说明Matlab 数值计算中的“左除”与“右除”有什么区别？2、脚本M 文件与函数M 文件的主要区别是什么？3、设矩阵A=[1 4 3 9；5 30 12 42；2 18 32 15]，试编程找出A 中所有大于或等于5的元素，并把这些元素组成一个列向量，并计算出这些元素的总和（注意附上标注说明）4、利用Matlab 编写程序求取335lim 21x x x →∞++和120()t t dt -⎰，并对求取结果进行分子和分母的提取（注意附上标注说明）。
5、利用Matlab 编写程序求微分方程x e y y y x 2cos 52=+'-''的通解。
二、 综合题（每题20分，共40分）1、在同一张图上分别画出两组坐标)10,,2,1(),4,(),,(3222 =+i i i i i i 的散点图和折线图。
要求只能利用plot 命令进行绘制，X 坐标范围[0,105] ，Y 坐标范围[0,1450]。
2、已知学生的名字和百分制分数。
编写相应的程序要求输入所有学生的百分制分数后，一次性判断每个学生的成绩分别是“满分”、“优秀”、“良好”、“及格”还是“不及格”，并要求最后用单元数组按行分别输出每个学生的名字、分数和成绩判断结果。
学生姓名 Jack Marry Peter Rose Tom成 绩 72 83 56 94 100判断标准：满分(100)，优秀(90-99)，良好(80-89)，及格(60-79)，不及格(<60)三、 论述题（共10分）浅谈一下Fortran 、matlab 、maple 和mathematic 这四种计算语言的特点和区别。
MATLAB第三章和第四章作业
习题三1.选择题(1)已知a=4，b='4'，下面说法错误的为C 。
A. 变量a比b占用的存储空间大B. 变量a和b可以进行加、减、乘、除运算C. 变量a和b的数据类型相同D. 变量b可以用eval命令执行(2)已知s='显示"hello"'，则s的元素个数是B 。
A. 12B. 9C. 7D. 18(3)运行命令>>a=sym('pi','d')则对于变量a的描述D是正确的。
A. a是符号变量B. a显示为10位的数值C. a显示为32位的数值D. a不存在(4)对于符号表达式g=sym('sin(a*z)+cos(w*v)')，如果运行命令diff(g)，则求导自变量是D。
A. aB. zC. wD. v(5)运行命令>>a=double(sym('sin(pi/2)'))，则变量a是 C。
A. 符号变量B. 字符串'1'C. double型的数值1D. 出错(6)运行命令>>y=dsolve('x*D2y-3*Dy=x^2','t')求解微分方程，则D。
A. Dy是指dy/dxB. 得出的y是通解有一个常数符C1C. D2y是指d2y/dxD. 得出的y是通解有两个常数符C1和C2(7)有一个2行2列的元胞数组c，则c(2)是指C 。
A. 第1行第2列的元素内容B. 第2行第1列的元素内容C. 第1行第2列的元素D. 第2行第1列的元素(8)对于一个2行2列的元胞数组c，如果要删除第二列的元素，应该使用的命令是B 。
A.c{:,2}=[]B.c(:,2)=[]C.A、B两种方式均可D.A、B两种方式均不可(9)有一个3行10列的构架数组student，每个构架有name（姓名）、scores两个字段，其中scores是以1×5数组表示的五门课的成绩，那么要查看第4个学生的第2门课成绩，采用的命令中达不到要求的是。
matlab作业
1、猜数游戏。
首先由计算机产生[1,100]之间的随机整数，然后由用户猜测所产生的随机数。
根据用户猜测的情况给出不同提示，如猜测的数大于产生的数，则显示“High”，小于则显示“Low”，等于则显示“You won”，同时退出游戏。
用户最多可以猜7次。
程序文件：nu m=fix(rand(1)*100);n=1;while n<=7a=input('请输入一个1到100的随机整数');n=n+1;if a>nu mdisp('High');elseif a<nu mdisp('low');elsedisp('won');break;endend2、从键盘输入5个数，求其中最大数和最小数。
并在屏幕上提示输入“请输入您要显示数字的序号”，输入后，回显“您要显示数字为……（此处为实际需要显示序号对应的数字）”。
程序文件：disp('请输入5个数字')for n=1:5;x=input('');A(1,n)=x;endAB=sort(A);Bdisp('最大值是');B(1,5)disp('最小值是');B(1,1)disp('请输入您要显示数字的序号');nu m=input('');disp('您要显示数字为');A(1,nu m)3、求Fibonacci数列(1)大于4000的最小项。
(2)5000之内的项数。
程序文件：（1）大于4000的最小项f=[1,1];j=1;while f(1,j+1)<4000f(1,j+2)=f(1,j)+f(1,j+1);j=j+1;endfdisp('大于4000的最小项是');f(1,j+1)运行结果f =Colu m ns 1 through 101 123 5 8 13 21 34 55Colu m ns 11 through 1989 144 233 377 610 987 1597 2584 4181大于4000的最小项是ans =4181（2）5000之内的项数f=[1,1];j=1;while f(1,j+1)<5000f(1,j+2)=f(1,j)+f(1,j+1);j=j+1;endfdisp('5000之内的项数');j运行结果f =Colu m ns 1 through 101 123 5 8 13 21 34 55Colu m ns 11 through 2089 144 233 377 610 987 1597 2584 4181 67655000之内的项数j =19。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第一章m=3880;%总质量f=0.013;%滚动半径r=0.367;%车轮半径CdA=2.77;%空气阻力系数*迎风面积io=5.83;%主减速器传动比lf=0.218;%飞轮转动惯性lw1=1.798;%二前轮转动惯量lw2=3.598;%四后轮转动惯性i1=5.56;i2=2.769;i3=1.644;i4=1;i5=0.793;L=3.2;%轴距a=1.947;%质心至前轴距离hg=0.9;%质心高n=linspace(400,4200);ua1=0.377*r*n/(io*i1);ua2=0.377*r*n/(io*i2);ua3=0.377*r*n/(io*i3);ua4=0.377*r*n/(io*i4);ua5=0.377*r*n/(io*i5);Tq=(-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/1 000).^4);Ft1=Tq*io*i1*0.85/r;Ft2=Tq*io*i2*0.85/r;Ft3=Tq*io*i3*0.85/r;Ft4=Tq*io*i4*0.85/r;Ft5=Tq*io*i5*0.85/r;F1=m*9.8*f+CdA/21.15*ua1.^2;F2=m*9.8*f+CdA/21.15*ua2.^2;F3=m*9.8*f+CdA/21.15*ua3.^2;F4=m*9.8*f+CdA/21.15*ua4.^2;F5=m*9.8*f+CdA/21.15*ua5.^2;plot(ua1,Ft1,'k',ua2,Ft2,'k',ua3,Ft3,'k',ua4,Ft4,'k',ua5,Ft5,'k',ua1,F1,'k',ua2,F2,'k',ua3,F 3,'k',ua4,F4,'k',ua5,F5,'k')title('汽车驱动力-行驶阻力平衡图');xlabel('ua(km/h)');ylabel('Ft(N)');gtext('Ft1');gtext('Ft2');gtext('Ft3');gtext('Ft4');gtext('Ft5');gtext('Ff+Fw');for b=1:3400;n(b)=600+b;Tq(b)=(-19.313+295.27*(n(b)/1000)-165.44*(n(b)/1000).^2+40.874*(n(b)/1000).^3-3 .8445*(n(b)/1000).^4);Ft5(b)=Tq(b)*i5*io*0.85/r;ua5(b)=0.377*r*n(b)/(io*i5);F5(b)=m*9.8*f+CdA/21.15*ua5(b).^2;E(b)=abs(Ft5(b)-F5(b));endfor b=1:3400;if(E(b)==min(E))disp('最高车速')disp(ua5(b))endendfor p=1:3400n(p)=600+p;Tq(p)=(-19.313+295.27*(n(p)/1000)-165.44*(n(p)/1000).^2+40.874*(n(p)/1000).^3-3 .8445*(n(p)/1000).^4);Ft1(p)=Tq(p)*i1*io*0.85/r;ua1(p)=0.377*r*n(p)/(io*i1);F1(p)=m*9.8*f+CdA/21.15*ua5(p).^2;af(p)=asin((Ft1(p)-F1(p))/(m*9.8));if(af(p)==max(af))imax=tan(af(p));C=imax*L/(a+hg*imax);endenddisp('最大爬坡度')disp(imax)disp('最大附着率')disp(C)最高车速99.088864831425923最大爬坡度0.338877068884551最大附着率0.481532745549641ig=[5.56 2.769 1.644 1.00 0.793];lw=lw1*2+lw2*4;for i=1:86n(i)=(14+i)*40;Tq(i)=(-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000).^2+40.874*(n(i)/1000).^3-3.8 445*(n(i)/1000).^4);endfor j=1:5for i=1:86Ft(i,j)=Tq(i)*ig(j)*io*0.85/r;ua(i,j)=0.377*r*n(i)/(ig(j)*io);Fc(i,j)=m*9.8*f+CdA/21.15*ua(i,j)^2;q(j)=1+lw/(m*r^2)+lf*ig(j)^2*io^2*0.85/(m*r^2);a(i,j)=(Ft(i,j)-Fc(i,j))/(m*q(j));if(a(i,j)<0.001)a(i,j)=a(i-1,j);endu1(i,j)=ua(i,j);b1(i,j)=1./a(i,j);endendplot(u1,1./a)axis([0 98 0 20])title('汽车加速度倒数曲线');xlabel('ua(km/h)')ylabel('1/a');for i =1:86A=ua(i,3)-ua(86,2);if(A>0&A<1)j=i;endB=ua(i,4)-ua(86,3);if(B>0&B<2);k=i;endif(ua(i,4)<=70)m=i;endendb=1./a;t=ua(1,2)*b(1,2);for p1=2:86t1(p1)=(ua(p1,2)-ua(p1-1,2))*(b(p1,2)+b(p1-1,2))/2;t=t+t1(p1);endfor p2=j:86t2(p2)=(ua(p2,3)-ua(p2-1,3))*(b(p2,3)+b(p2-1,3))/2;t=t+t2(p2);endfor p3=k:mt3(p3)=(ua(p3,4)-ua(p3-1,4))*(b(p3,4)+b(p3-1,4))/2;endt=t+(ua(j,3)-ua(86,2))*b(86,2)+(ua(k,4)-ua(86,3))*b(86,3)+(70-ua(m,4))*b(m,4); tz=t/3.6;disp('加速时间');disp(tz);加速时间18.982570197107204第二章m=3880;%总质量f=0.013;%滚动半径r=0.367;%车轮半径CdA=2.77;%空气阻力系数*迎风面积io=5.83;%主减速器传动比lf=0.218;%飞轮转动惯性lw1=1.798;%二前轮转动惯量lw2=3.598;%四后轮转动惯性ig=[5.56 2.769 1.644 1.00 0.793];L=3.2;%轴距a=1.947;%质心至前轴距离hg=0.9;%质心高nd=400;%怠速转速Qid=0.299;%怠速油耗g=9.8;%重力加速度m1=2000;%装载质量m2=1800;%整车整备质量n1=[815 1207 1614 2012 2603 3006 3403 3804];pg=7.05;for j=1:5for i=1:69n(i)=(i+11)*50;Tq(i)=-19.313+295.27*(n(i)/1000)-165.44*(n(i)/1000).^2+40.874*(n(i)/1000).^3-3.84 45*(n(i)/1000).^4;Pe(i)=Tq(i)*n(i)/9550;ua(i,j)=0.377*r*n(i)/(ig(j)*io);Pc(i,j)=(m*g*f*ua(i,j)/3600+CdA*ua(i,j)^3/76140)/0.85;endendplot(ua,Pe,ua,Pc);axis([0 120 0 100])title('功率平衡图');xlabel('ua(km/h)');ylabel('p(kw)');gtext('Pe1');gtext('Pe2');gtext('Pe3');gtext('Pe4');gtext('Pe5');gtext('Pc');for i=1:8Tq(i)=-19.313+295.27*(n1(i)/1000)-165.44*(n1(i)/1000).^2+40.874*(n1(i)/1000).^3-3.8445*(n1(i)/1000).^4;Pe(i)=n1(i)*Tq(i)/9550;ua4(i)=0.377*r*n1(i)/(ig(4)*io);ua5(i)=0.377*r*n1(i)/(ig(5)*io);endb(1)=1326.8-416.46*Pe(1)+72.379*Pe(1)^2-5.8629*Pe(1)^3+0.17768*Pe(1)^4;b(2)=1354.7-303.98*Pe(2)+36.657*Pe(2)^2-2.0553*Pe(2)^3+0.043072*Pe(2)^4;b(3)=1284.4-189.75*Pe(3)+14.524*Pe(3)^2-0.51184*Pe(3)^3+0.0068164*Pe(3)^4; b(4)=1122.9-121.59*Pe(4)+7.0053*Pe(4)^2-0.18517*Pe(4)^3+0.0018555*Pe(4)^4; b(5)=1141.0-98.893*Pe(5)+4.4763*Pe(5)^2-0.091077*Pe(5)^3+0.00068906*Pe(5)^4; b(6)=1051.2-73.714*Pe(6)+2.8593*Pe(6)^2-0.05138*Pe(6)^3+0.00035032*Pe(6)^4; b(7)=1233.9-84.478*Pe(7)+2.9788*Pe(7)^2-0.047449*Pe(7)^3+0.00028230*Pe(7)^4 ;b(8)=1129.7-45.291*Pe(8)+0.71113*Pe(8)^2-0.00075215*Pe(8)^3-0.000038568*Pe( 8)^4;b=[b(1) b(2) b(3) b(4) b(5) b(6) b(7) b(8)];B4=polyfit(ua4,b,3);B5=polyfit(ua5,b,3);for j=1:69bh4(j)=polyval(B4,ua(j,4));bh5(j)=polyval(B5,ua(j,5));Q4(j)=Pc(j,4)*bh4(j)/(1.02*ua(j,4)*pg);Q5(j)=Pc(j,5)*bh5(j)/(1.02*ua(j,5)*pg); endplot(ua(:,4),Q4,ua(:,5),Q5);title('4档5档等速百公里油耗图');xlabel('ua(km/h)');ylabel('Q(L/100km)');gtext('4档');gtext('5档');第三章m1=2000;m2=1800;mz=3880;g=9.81;r=0.367;CdA=2.77;f=0.013;nT=0.85;If=0.218;Iw1=1.798;Iw2=3.598;L=3.2;a=1.947;hg=0.9;Qid=0.299;i0=[5.17 5.43 5.83 6.17 6.33];ig=[5.56 2.769 1.644 1.00 0.793];Iw=2*Iw1+4*Iw2;for k=1:5n=[815 1207 1614 2012 2603 3006 3403 3804];ua=0.377*n*r/(ig(5)*i0(k));Pe=(mz*g*f*ua/3600.+CdA*ua.^3/76140.)/nT;B0=[1326.8,1354.7,1284.4,1122.9,1141.0,1051.2,1233.9,1129.7];B1=[-416.46,-303.98,-189.75,-121.59,-98.893,-73.714,-84.478,-45.291];B2=[72.379,36.657,14.524,7.0035,4.4763,2.8593,2.9788,0.71113];B3=[-5.8629,-2.0553,-0.51184,-0.18517,-0.091077,-0.05138,-0.047449,-0.00075215]; B4=[0.17768,0.043072,0.0068164,0.0018555,0.00068906,0.00035032,0.00028230,-0 .000038568];b=B0+B1.*Pe+B2.*Pe.^2+B3.*Pe.^3+B4.*Pe.^4;ua1=25;Pe1=(mz*g*f*ua1/3600.+CdA*ua1^3/76140.)/nT;b1=interp1(ua,b,ua1,'spline');Qs1(k)=Pe1*b1*50/(102*ua1*7.05);ua2=(25+15/16.7/2):15/16.7:(40-15/16.7/2);b2=interp1(ua,b,ua2,'spline');DIt=1+Iw/(mz*r^2)+If*ig(5).^2*i0(k)^2*nT/(mz*r^2);Pz2=(mz*g*f*ua2/3600.+CdA*ua2.^3/76140+DIt*mz*ua2*0.25/3600)/nT;Qs2(k)=sum(Pz2.*b2/367.1/7.05);ua3=40;b3=interp1(ua,b,ua3,'spline');Pz3=(mz*g*f*40/3600+CdA*40^3/76140)/nT;Qs3(k)=Pz3*b3*250/(102*ua3*7.05);ua4=(40+10/14):10/14:(50-10/14);b4=interp1(ua,b,ua4,'spline');DIt=1+Iw/(mz*r^2)+If*ig(5).^2*i0(k)^2*nT/(mz*r^2);Pz4=(mz*g*f*ua4/3600.+CdA*ua4.^3/76140+DIt*mz*ua4*0.25/3600)/nT;Qs4(k)=sum(Pz4.*b4/367.1/7.05);ua5=50;b5=interp1(ua,b,ua5,'spline');Pz5=(mz*g*f*50/3600.+CdA*50^3/76140.)/nT;Qs5(k)=Pz5*b5*250/(102*ua5*7.05);Qs6(k)=(50-25)*Qid/3.6/0.36;Qs(k)=(Qs1(k)+Qs2(k)+Qs3(k)+Qs4(k)+Qs5(k)+Qs6(k))/1000/1.075*100; endfor k=1:5ig0=[5.56;2.769;1.644;1.00;0.793];ua0=0:0.01:120;[ua,ig]=meshgrid(ua0,ig0);n=(ig*i0(k)).*ua/0.377/r;n(find(n<600|n>4000))=NaN;Tq=-19.313+295.27*(n/1000)-165.44*(n/1000).^2+40.874*(n/1000).^3-3.8445*(n/10 00).^4;Ft=ig.*Tq*i0(k)*nT/0.367;Fz=CdA*ua.^2/21.15+mz*g*f;DIt=1+Iw/(mz*r^2)+If*ig.^2*i0(k)^2*nT/(mz*r^2);A=(Ft-Fz)./(DIt*mz);A_=min(min(1./A(2,:),1./A(3,:)),min(1./A(4,:),1./A(5,:)));t(k)=sum(A_(find(ua0<=70&A_>0)))*0.01/3.6;endQs_=14.8:0.01:15.9;t_=polyval(polyfit(Qs,t,4),Qs_);plot(Qs,t,'*',Qs_,t_);title('燃油经济性加速时间曲线')xlabel('百公里油耗')ylabel('加速时间')第四章m1=4080;hg1=0.845;a1=2.100;%空载m2=9290;hg2=1.170;a2=2.950;%满载L=3.950;%轴距b=0.38;%制动力分配系数z=linspace(0,1,21);vf1=b.*z*L./(L-a1+z*hg1);vr1=(1-b).*z*L./(a1-z*hg1);vf2=b.*z*L./(L-a2+z*hg2);vr2=(1-b).*z*L./(a2-z*hg2);v=z;for i=1:21%ECE制动法规if(z(i)<0.3&&z(i)>0.15)v3(i)=z(i)+0.08;endif(z(i)>=0.3)v3(i)=0.38+(z(i)-0.3)/0.74;endendz1=0.15:0.01:0.3;v4=z1-0.08;plot(z,vf1,z,vr1,z,vf2,':',z,vr2,':',z,v,z,v3,'o',z1,v4,'+') axis([0,1,0,1.2])title('利用附着系数与制动强度的关系曲线')xlabel('制动强度z/g')ylabel('利用附着系数v')gtext('空车前轴')gtext('空车后轴')gtext('满载前轴')gtext('满载后轴')gtext('ECE法规')Er1=(a1/L)./((1-b)+vr1*hg1/L)*100; Ef=(L-a2)/L./(b-vf2*hg2/L)*100;Er=(a2/L)./((1-b)+vr2*hg2/L)*100; plot(vr1,Er1,vr2,Er,vf2,Ef)axis([0,1,0,100])title('前后附着效率曲线')xlabel('附着系数v')ylabel('制动效率（%）')gtext('满载')gtext('Ef')gtext('Er')gtext('Er空载')C=0.8;%前轴先抱死ua=30;%行驶车速t1=0.02;%制动系反应时间t2=0.02;%制动减速度上升时间E1=(a1/L)./((1-b)+C*hg1/L);E2=(a2/L)./((1-b)+C*hg2/L);g1=E1*C*9.8;g2=E2*C*9.8;s1=(t1+t2/2)*ua/3.6+ua^2/(25.92*g1); s2=(t1+t2/2)*ua/3.6+ua^2/(25.92*g2); disp('空载时不抱死的制动距离')disp(s1)disp('满载时不抱死的制动距离')disp(s2)运行结果：空载时不抱死的制动距离 6.840557850124176满载时不抱死的制动距离5.331923306045581b1=0;%前部管路损坏b2=1;%后部管路损坏Ekf=(L-a1)/L./(b2-C*hg1/L);akf=Ekf*C*9.8;skf=(t1+t2/2)*ua/3.6+ua^2/(25.92*akf);disp('空载时后部管路损坏时制动距离')disp(skf)Ekr=(a1/L)./((1-b1)+C*hg1/L);akr=Ekr*C*9.8;skr=(t1+t2/2)*ua/3.6+ua^2/(25.92*akr);disp('空载时前部管路损坏时制动距离')disp(skr)Emf=(L-a2)/L./(b2-C*hg2/L);amf=Emf*C*9.8;smf=(t1+t2/2)*ua/3.6+ua^2/(25.92*amf);disp('满载时后部管路损坏时制动距离')disp(smf)Emr=(a2/L)./((1-b1)+C*hg2/L);amr=Emr*C*9.8;smr=(t1+t2/2)*ua/3.6+ua^2/(25.92*amr);disp('满载时前部管路损坏时制动距离')disp(smr)运行结果：空载时后部管路损坏时制动距离： 8.087876141447570空载时前部管路损坏时制动距离：10.006134596695821满载时后部管路损坏时制动距离：13.598568594104307满载时前部管路损坏时制动距离：7.585384718859294第五章m=1818.2;lz=3885;L=3.048;a=1.463;b=1.585;k1=-62618;k2=-110185;i=20;disp('(1)稳定性因数K、特征车速uch为：')K=m/L^2*(a/k2-b/k1)uch=sqrt(1/K)u=0:0.1:180;Wr_delt=u/L./(1+K*u.^2);Wr_delt0=u/L;plot(u,Wr_delt,'r',u,Wr_delt0,'b')axis([0 180 0 25])xlabel('车速u（km/h）');ylabel('稳态横摆增益Wr/delt')title('稳态横摆角速度增益曲线')text(65,20,'中性转向K=0');text(100,3,'不足转向K=0.0024'); disp('(2)车速为22.35m/s时的转向灵敏度为：')22.35*3.6/L/(1+K*(22.35*3.6)^2);disp('(3)静态储备系数S.M')SM=k2/(k1+k2)-a/Ldisp('前、后轮侧偏角绝对值只差a1-a2=')ay=0.4*9.8;R0=15;K*ay*Ldisp('转弯半径的比值R/R0=')delt=L/R0;L/(delt-K*ay*L)/R0disp('(4).固有圆频率')u1=30.56*3.6;W0=L/u1*sqrt(k1*k2*(1+K*u1^2)/(m*lz))disp('阻尼比为')f=(-m*(a^2*k1+b^2*k2)-lz*(k1+k2))/(2*L*sqrt(m*lz*k1*k2*(1+K *u1^2)))disp('反应时间')t=atan(sqrt(1-f^2)/(-m*u1*a*W0/(L*k2)-f))/(W0*sqrt(1-f^2)) disp('峰值反应时间')e=atan(sqrt(1-f^2)/f)/(W0*(1-f^2))+t(1)稳定性因数K、特征车速uch为：K= 0.002355273063120Uch=20.605320186515470(2)车速为22.35m/s时的转向灵敏度为：Ans=1.6247(3)静态储备系数S.MSM =0.157646729257626前、后轮侧偏角绝对值只差a1-a2=ans=0.028141179401843转弯半径的比值R/R0=ans =1.160752707608149(4).固有圆频率W0 =4.703453307131383 阻尼比为f =0.194016603943443 反应时间t =0.053348915482336 峰值反应时间e =0.357241690032395第六章f0=1.5;g=0.25;r=9;u=10;fs=3;gs=0.25;;n0=0.1;df=0.2;N=180;f=0:0.2:36;ff0=f/f0;d=[(1-ff0.^2).*(1+r-1/u*ff0.^2-1)].^2+4*g.^2*ff0.^2.*[r-(1/u+1)*ff0.^ 2].^2;z1q=r*sqrt(((1-ff0.^2).^2+4*g^2*ff0.^2)./d);z2z1=sqrt((1+(2*g*ff0).^2)./((1-ff0.^2).^2+(2*g*ff0).^2));ffs=f/fs;pz2=sqrt((1+(2*gs*ffs).^2)./((1-ffs.^2).^2+(2*gs*ffs).^2));loglog(f,z1q,'*',f,z2z1,'.',f,pz2);axis([0 100 0 10])set(gca,'xtick',[0 0.1 1 10 100])set(gca,'ytick',[0 0.1 1 10])xlabel('激振频率f/Hz')gtext({'人体-车身部分的幅频特性';'车轮部分的幅频特性';'车身部分-车轮部分的幅频特性'});Gqn0=2.56e-4w=2*pi*f;Gz1=w.^2.*z1q.*(2*3.14./w).*sqrt(Gqn0^2*u);loglog(f,Gz1)hold onGz2=w.^2.*z1q.*z2z1.*(2*3.14./w).*sqrt(Gqn0^2*u);loglog(f,Gz2)Gza=w.^2.*z1q.*z2z1.*pz2.*(2*3.14./w).*sqrt(Gqn0^2*u);loglog(f,Gza)title('加速度均方根值')xlabel('激振频率f/Hz')ylabel('|Gz|')gtext({'车轮部分';'车身部分';'人体'})format longGq=2.56e-008;n0=0.1;f0=1.5;%固有频率g=0.25;%阻尼比r=9;%刚度比u=10;%质量比fs=3;%人体座椅系统的固有频率gs=0.25;%人体座椅系统的阻尼比ua=20;%车速Df=0.2;%频率步长N=180;%频率点数f=0.000001:0.2:36;ff0=f/f0;w=2*pi*f;d=((1-ff0.^2).*(1+r-ff0.^2/u)-1).^2+4.*g^2.*ff0.^2.*(r-(1/u+1).*ff0.^2).^2;z1q=r*sqrt(((1-ff0.^2).^2+4.*g^2.*ff0.^2)./d);z2z1=sqrt((1+(2*g*ff0).^2)./((1-ff0.^2).^2+(2*g*ff0).^2));ffs=f/fs;pz2=sqrt((1+(2*gs*ffs).^2)./((1-ffs.^2).^2+(2*gs*ffs).^2));Ez=sqrt(sum(16.*pi.^4.*f.^2.*Gq.*n0.^2.*ua*Df));disp('路面不平度加速度均方根值=')disp(Ez)Ez1=sqrt(sum(z1q.^2.*16.*pi^4.*Gq.*ua.*f.^2*Df*n0^2));disp('车轮加速度均方根值=')disp(Ez1)Ez2=sqrt(sum(z2z1.^2.*z1q.^2.*16.*pi^4.*Gq.*ua.*f.^2*Df*n0^2));disp('车身加速度均方根值=')disp(Ez2)Gza=w.^2.*z1q.*z2z1.*pz2.*(2*pi./sqrt(Gq*n0^2*ua));Eza=sqrt(sum(Gza.^2*Df))disp('人体加速度均方根值=')disp(Eza)Gza=w.^2.*z1q.*z2z1.*pz2.*(2*pi./w).*sqrt(Gq*n0^2*ua);for k=1:180if(f(k)<2&&f(k)>0.5)f1(k)=f(k);Gza(k)=w(k).^2*z1q(k).*z2z1(k).*pz2(k).*(2*3.14./w(k)).*sqrt(Gq*n0^2* ua);a1(k)=0.5^2*Gza(k).^2*0.2;endif(f(k)>2&&f(k)<4)a2(k)=(f(k)/4).^2*Gza(k).^2;f2(k)=f(k);endif(f(k)>4&&f(k)<12.5)a3(k)=Gza(k).^2;f3(k)=f(k);endif(f(k)>12.5&&f(k)<80)a4(k)=12.5^2/f(k)*Gza(k).^2;f4(k)=f(k);endendaw1=trapz(f1,a1);aw2=trapz(f2,a2);aw3=trapz(f3,a3);aw4=trapz(f4,a4);aw=aw1+aw2+aw3+aw4disp('人体加权加速度均方根值=') disp(aw)Law=20*log10(aw*10^6)disp('加权振级=')disp(Law)路面不平度加速度均方根值= 0.350811721645482车轮加速度均方根值=0.239075671885801车身加速度均方根值=0.016788100118782Eza =5.342310914121454e+007 人体加速度均方根值=5.342310914121454e+007 aw =1.527317*********e-004 人体加权加速度均方根值=1.527317*********e-004 Law =43.678586277719972加权振级=43.678586277719972format longGq=2.56e-008;f0=1.5;%固有频率g=0.25;%阻尼比r=9;%刚度比u=10;%质量比fs=3;%人体座椅系统的固有频率gs=0.25;%人体座椅系统的阻尼比ua=20;%车速Df=0.2;%频率步长N=180;%频率点数f=0.000001:0.2:36;ff0=f/f0;w=2*pi*f;fs1=linspace(1.5,6,100);gs2=linspace(0.125,0.5,10);for i=1:100d=((1-ff0.^2).*(1+r-ff0.^2/u)-1).^2+4.*g^2.*ff0.^2.*(r-(1/u+1).*ff0.^ 2).^2;z1q=r*sqrt(((1-ff0.^2).^2+4.*g^2.*ff0.^2)./d);z2z1=sqrt((1+(2*g*ff0).^2)./((1-ff0.^2).^2+(2*g*ff0).^2));ffs1=f/fs1(i);pz2=sqrt((1+(2*gs*ffs1).^2)./((1-ffs1.^2).^2+(2*gs*ffs1).^2)); Gza=w.^2.*z1q.*z2z1.*pz2.*(2*pi./w).*sqrt(Gq*n0^2*ua);for k=1:180if(f(k)<2&&f(k)>0.5)f1(k)=f(k);Gza(k)=w(k).^2*z1q(k).*z2z1(k).*pz2(k).*(2*3.14./w(k)).*sqrt(Gq*n0^2* ua);a1(k)=0.5*Gza(k).^2*0.2;endif(f(k)>2&&f(k)<4)a2(k)=f(k)/4*Gza(k).^2;f2(k)=f(k);endif(f(k)>4&&f(k)<12.5)a3(k)=Gza(k).^2;f3(k)=f(k);endif(f(k)>12.5&&f(k)<80)a4(k)=12.5/f(k)*Gza(k).^2;f4(k)=f(k);endaw1=trapz(f1,a1);aw2=trapz(f2,a2);aw3=trapz(f3,a3);aw4=trapz(f4,a4);awa(i)=sqrt(aw1+aw2+aw3+aw4);Law1(i)=20*log10(awa(i)*10^6);endfor m=1:10ffs=f/fs;d=((1-ff0.^2).*(1+r-ff0.^2/u)-1).^2+4.*g^2.*ff0.^2.*(r-(1/u+1).*ff0.^ 2).^2;z1q=r*sqrt(((1-ff0.^2).^2+4.*g^2.*ff0.^2)./d);z2z1=sqrt((1+(2*g*ff0).^2)./((1-ff0.^2).^2+(2*g*ff0).^2));pz2=sqrt((1+(2*gs2(m)*ffs).^2)./((1-ffs.^2).^2+(2*gs2(m)*ffs).^2)); Gza=w.^2.*z1q.*z2z1.*pz2.*(2*pi./w).*sqrt(Gq*n0^2*ua);for h=1:180if(f(h)<2&&f(h)>0.5)f1(h)=f(h);Gza(h)=w(h).^2*z1q(h).*z2z1(h).*pz2(h).*(2*3.14./w(h)).*sqrt(Gq*n0^2* ua);a1(h)=0.5*Gza(h).^2*0.2;endif(f(h)>2&&f(h)<4)a2(h)=f(h)/4*Gza(h).^2;f2(h)=f(h);endif(f(h)>4&&f(h)<12.5)a3(h)=Gza(h).^2;f3(h)=f(h);endif(f(h)>12.5&&f(h)<80)a4(h)=12.5/f(h)*Gza(h).^2;f4(h)=f(h);endendaw1=trapz(f1,a1);aw2=trapz(f2,a2);aw3=trapz(f3,a3);aw4=trapz(f4,a4);awb(m)=sqrt(aw1+aw2+aw3+aw4);Law2(m)=20*log10(awb(m)*10^6);endsubplot(2,2,1)plot(fs1,awa)xlabel('座椅fs1/Hz')ylabel('人体加权加速度均方根值')gtext('aw1')subplot(2,2,2)plot(fs1,Law1)xlabel('座椅fs1/Hz')ylabel('加权振级')gtext('Law1')subplot(2,2,3)plot(gs2,awb)xlabel('座椅gs2')ylabel('人体加权加速度均方根值')gtext('aw2')subplot(2,2,4)plot(gs2,Law2)xlabel('座椅gs2')ylabel('加权振级')gtext('Law2')人体加权加速度均方根值和加权振级随fs1的增大而增大。