MATLAB大作业
MATLAB大作业
1.设计一个程序，要求用户用键盘输入10次，输入内容可以是数字、英文、汉字、或者
其他文字或者符号。
如果输入的是数字，则把存入一个数字数组当中；如果输入的是英文的字符，则存入一个英文字符数组当中；如果是其他字符，则存入一个字符串数组中。
当10次内容都输入完毕后，先输出英文数组中的内容，然后输出字符串数组中的内容。
如果数字数组中元素数量大于3个，则由数字数组中的数字组成一个a0 X n+a1X n-1+………an=f（X），其中a0…an是数组中第一到最后一个元素，n为元素的个数，并且画出这个函数在区间[0,20]之间的曲线。
如果数字数组中的元素数量小于等于3个，则由数字数组中的数字组成一个a0 X n+a1Y n-1+anZ=0的方程，并且画出这个函数在[0,10]之间的曲线。
具体要求：
（1）画出程序的结构图
（2）画出程序的流程图
（3）每一行或者每一小段都需要标注
（4）附上完整的程序
2.按如下条件，设计一个完整的程序。
（1）有一片草原，草原上有人、狼、羊、鸡和蝗虫生活；
（2）被吃过的草，三年之后可以修复；
（3）一个人一年吃一只羊，10个人一年杀一只狼，一个人一年吃一只鸡；
（4）一头狼一年吃一只羊，一头狼三年成活一个幼崽；
（5）一只羊一年吃掉100平方米的草原；
（6）一只鸡一年吃掉1000只蝗虫；
（7）一只蝗虫一年吃掉一平方米的草原。
按如上条件，模拟出草原及其上述的生物从现在起500年间的状态。
具体要求：
（1）写出具体的解题思路
（2）列出详细的变量表，写出每一个变量的含义
（3）写出完整的程序，最好有图。
2023年MATLAB大作业题目备选
ＭATLAB大作业备选题目1、基于MＡTLＡB旳有噪声语音信号处理本课题规定基于ＭＡTLAB对有噪音语音信号进行处理,综合运用数字信号处理旳理论知识对加噪语音信号进行时域、频域分析和滤波,运用MATLAＢ作为工具进行计算机实现。
在设计实现旳过程中,规定使用双线性变换法设计ＩIR数字滤波器，对模拟加噪语音信号进行低通滤波、高通滤波及带通滤波,并运用MA ＴLＡＢ作为辅助工具完毕设计中旳计算与图形旳绘制。
2、基于MＡTLAB旳学生平均学分、绩点计算软件设计学分与绩点,是每位大学生所关怀旳重要指标之一，诸多同学辛劳学习，早出晚归,不停旳奔走于教室、图书馆、食堂、寝室之间,为旳就是可以考个好成绩，获得好旳绩点。
然而在平时我们计算学分与绩点旳时候,大都只能用计算器一种一种数据旳输入，其过程繁琐麻烦,又轻易出错。
因此，本课题规定运用所学旳MATＬＡB知识，来实现平均学分、绩点旳计算，并开发有关人机界面。
3、基于MＡTLAB旳试卷分析管理系统本设计规定基于MＡＴLＡＢ中GUＩ旳编程措施,并波及有关数据库知识。
规定通过一种简易旳顾客交互界面，实现对考试试卷旳成绩录入、查询、修改和试题整体分析等功能，以学习使用ＭＡＴLAB编程为目旳,尤其是对MATＬＡＢ中G ＵI旳掌握,加深对MＡTＬＡB旳理解,学习用ＭATＬAB实现实际应用。
4、基于ＭＡＴLAＢ旳图像处理软件设计学习MATLAB GUI程序设计,运用MATＬAＢ图像处理工具箱,设计和实现一种简易旳图像处理软件，实现如下几点功能：1）图像旳读取和保留。
2）设计图形顾客界面，让顾客可以对图像进行任意旳亮度和对比度变化调整,显示和对比变换前后旳图像。
３）设计图形顾客界面,让顾客可以用鼠标选用图像感爱好区域，显示和保留该选择区域。
4)编写程序通过近来邻插值和双线性插值等算法将顾客所选用旳图像区域进行放大和缩小整数倍旳操作,并保留,比较几种插值旳效果。
5)图像直方图记录和直方图均衡，规定显示直方图记录，比较直方图均衡后旳效果。
Matlab大作业
Matlab 大作业（组内成员：彭超杰、南彦东、江明伟）一、研究模型（电车）通过控制油门（保持一定角度）来调节电动机能输出稳定的转速，从而控制车速稳定。
数学依据说明如下：由图可知存在以下关系：a d a au w k R i dtdi L =++ (w k e d d =) L M M dtdw J-=a m i k M =La m M i k dtdw J -=k为反电势常数，m k为电动机电磁力矩常数，这里忽略阻尼力矩。
d二、数学模型再看整个研究对象，示意图以课本为依据，不同点是这里将数控的进给运动，转换为汽车行驶所需要的扭矩。
（这里不说明扭矩的具体产生过程，仅仅说明输出车轮旋转的角速度w ）对照课本不同，() s θ变为()s N ，1221z z w w =，1w 为电动机的转速，2w 为轮胎的转速，1z 为电动机的光轴齿轮的齿数，2z 为与轮胎相连光轴的齿轮齿数。
)(*10110w x w k x ==，121z z k =()c a m m d ba m x K K K k s k k JRs JLs K K K k s G i 1231+++=()ca m m d M K K K k s k k JRs JLs R Ls K s G L 1231)(++++-=同理，忽略电枢绕组的电感L ，简化系统传递函数方框图如下()JRK K K k JR s k k s JR K K K k s G c a m m d ba m x i 121++=()JRK K K k JR s k k s K K K K k s k k Rs R K s G c a m m d c a m m d M L 121121++-=++-=三、系统分析1.分析时间响应其传递函数如下：（1）系统时间响应令τ=0、τ=0.0125、τ=0.025，应用impulse函数，可得到系统单位脉冲响应；应用step函数，可得系统单位跃阶响应。
MATLAB期末大作业模板
MATLAB应用技术期末大作业专业：姓名：学号：分数一、在一个图形窗口中以子图形式同时绘制正弦、余弦、正切、余切曲线。
请写下完整代码，展示图形结果。
（请标注题图和坐标轴，用不同颜色和不同线型分别绘制以上曲线）。
（15分）二、某公司员工的工资计算方法如下。
（1）工作时数超过120小时者，超过部分加发15%。
（2）工作时数低于60小时者，扣发700元。
（3）其余按每小时84元发。
根据员工的工时数，计算应发工资。
请写下完整的程序代码，并任意输入一工时数（使用input 函数），将结果展示（使用disp 函数）利用该代码进行计算工资，请写下计算结果。
（15分）三、编写一个函数文件，使其能够产生如下的分段函数：⎪⎩⎪⎨⎧≥<<≤-=66225.0,25.05.15.0)(x x x x xx f请编写完整的函数文件（保存函数文件名为hanshu.m ），并编写脚本文件代码，任意输入x 值（使用input 函数），在脚本文件中调用函数文件求)(x f ，展示结果（使用disp 函数），请写下计算结果。
（15分）四、将5个学生的6门功课的成绩存入矩阵P 中，进行如下处理：（1）分别求每门课的最高分、最低分及相应学生的序号。
（2）分别求每门课的平均分和标准差。
（3）5门课总分的最高分、最低分及相应学生序号。
（4）将5门课总分按从大到小顺序存入score 中，相应学生序号存入num 。
请将各小题的运行代码完整写下来，并写下运行结果。
（20分）五、请利用所学的MATLAB 知识，自主设计一个图形用户界面，请完整记录它的设计过程，需提供文字、代码和图片，以充分说明设计的图形用户界面可实现的功能。
（35分）。
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.编辑模块即将各个模块连接起来。
matlab大作业
2014年春季学期MATLAB 课程考查题一.必答题(80分)1.如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有何区别？(2)答：方法一：在MATLAB命令窗口中输入editpath或pathtool命令或通过【File】/【SetPath】菜单，进入“设置搜索路径”对话框，通过该对话框编辑搜索路径。
方法二：在命令窗口执行“path(path,…D:\Study ‟)”，然后通过“设置搜索路径”对话查看“D:\Study”是否在搜索路径中。
方法三：在命令窗口执行“addpathD:\Study- end”，将新的目录加到整个搜索路径的末尾。
如果将end改为begin，可以将新的目录加到整个搜索路径的开始。
区别：当前文件目录是正在运行的文件的目录，显示文件及文件夹的详细信息，且只有将文件设置为当前目录才能直接调用。
搜索路径中的文件可以来自多个不同目录，在调用时不用将其都设置为当前目录，为同时调用多个文件提供方便。
2.创建符号变量和符号表达式有哪几种方法？(4)答：定义符号变量：方法一：sym函数，可以定义单个符号变量，调用格式为：符号量名=sym（'符号字符串'）；方法二，syms函数，定义多个符号变量，调用格式：syms 符号变量名1 符号变量名2 符号变量名3 符号变量名4 。
空格隔开。
定义符号表达式：方法一：用单引号定义符号表达式；方法二：用sym函数定义符号表达式；方法三：用已经定义好的符号变量组成符号表达式；3.GUIDE提供哪些常用的控件工具，各有什么功能？（5分）答：一、控件风格和外观（1）BackgroundColor：设置控件背景颜色，使用[R G B]或颜色定义。
（2）CData：在控件上显示的真彩色图像，使用矩阵表示。
（3）ForegroundColor：文本颜色。
（4）String属性：控件上的文本，以及列表框和弹出菜单的选项。
（5）Visible：控件是否可见。
MATLAB期末大作业
学号：姓名：《Matlab/Simulink在数学计算与仿真中的应用》大作业1.假设地球和火星绕太阳运转的半径分别为r和2r，利用comet指令动画显示从地球到火星的转移轨迹（r可以任意取值，要求实时显示探测器、太阳、地球和火星的位置）。
解函数function comet(varargin)[ax,args,nargs] = axescheck(varargin{:});error(nargchk(1,3,nargs,'struct'));% Parse the rest of the inputsif nargs < 2, x = args{1}; y = x; x = 1:length(y); endif nargs == 2, [x,y] = deal(args{:}); endif nargs < 3, p = 0.10; endif nargs == 3, [x,y,p] = deal(args{:}); endif ~isscalar(p) || ~isreal(p) || p < 0 || p >= 1error('MATLAB:comet:InvalidP', ...'The input ''p'' must be a real scalar between 0 and 1.'); End指令 %particle_motiont = 0:5:16013;r1=6.7e6;%随便给定参数%---------------------------r2=2*r1;g=9.8;R=6.378e6;m=g*R^2;%内轨道v_inner=sqrt(m/r1);w_inner=v_inner/r1;x_inter=r1*cos(w_inner*t);y_inter=r1*sin(w_inner*t);%外轨道v_outer=sqrt(m/r2);w_outer=v_outer/r2;x_outer=r2*cos(w_outer*t);y_outer=r2*sin(w_outer*t);%控制器转移轨道a=(r1+r2)/2;E=-m/(2*a);V_near=sqrt(m*(2/r1-2/(r1+r2)));%转移轨道在近地点的速度V_far=sqrt(m*(2/r2-2/(r1+r2)));%转移轨道在远地点的速度h=r1*V_near;%由于在近地点的速度垂直于位置失量, h是转移轨道的比动量矩e=sqrt(1+2*E*h^2/m^2);%e为椭圆轨迹的偏心率TOF=pi*sqrt(a^3/m);%转移轨道是椭圆的一半及飞行时间是周期的一半(开普勒第三定律)w=pi/TOF;%椭圆轨迹的角速度c=a*e;b=sqrt(a^2-c^2);x_ellipse=a*cos(w*t)-0.5*r1;y_ellipse=b*sin(w*t);%动画显示运动轨迹x=[x_inter;x_outer;x_ellipse]';y=[y_inter;y_outer;y_ellipse]';comet(x,y)%---------------------------动态图像如下：2．利用两种不同途径求解边值问题dfdxf gdgdxf g f g=+=-+==34430001,,(),()的解.途径一：指令syms f g[f,g]=dsolve('Df=3*f+4*g,Dg=-4*f+3*g','f(0)=0,g(0)=1');disp('f=');disp(f)disp('g=');disp(g)结果（Matlab 7.8版本）f=i/(2*exp(t*(4*i - 3))) - (i*exp(t*(4*i + 3)))/2g=exp(t*(4*i + 3))/2 + 1/(2*exp(t*(4*i - 3)))（Matlab 6.5版本）f=exp(3*t)*sin(4*t)g=exp(3*t)*cos(4*t)>>途径二：%problem2function dy=problem2(t,y)dy = zeros(2,1);dy(1) = 3*y(1)+4*y(2);dy(2) = -4*y(1)+3*y(2);[t,y] = ode45('problem2',[0 2],[0 1]);plot(t,y(:,1),'r',t,y(:,2),'b');图23．假设著名的Lorenz 模型的状态方程表示为⎪⎩⎪⎨⎧-+-=+-=+-=)()()()()()()()()()()()(322133223211t x t x t x t x t x t x t x t x t x t x t x t x σρρβ 其中，设28,10,3/8===σρβ。
MATLAB与控制系统仿真大作业
>> G1=tf(num,den);
>> G1=zpk(G1)
G1 =
5 (s+0.6) ----------------(s+3) (s+2) (s+1)
>> num=[2 1]; den=[1 2.9 1]; G2=tf(num,den); >> G2=zpk(G2)
G2 =
2 (s+0.5)
>> xlabel('x') ylabel('exp(-2.*x)')
三、simulink 建模（每题 10 分，共 20 分）
1、已知单位负反馈的开环传递函数为 G(s) =
2 s2 + 4s
，试利用
simulink
建立系在单位阶跃输入作用下的模型。
要求答案包括：（1）simulink 建模结构图；
（2）在同一个坐标中的阶跃信号和响应曲线图。
2、已知系统的开环传递函数为
G1(s)
=
2s2 + 5s + 6 s2 + 2s + 3
、G2 (s)
=
s2
s+6 + 7s +1
，H
(s)
=
5(s + 2) s +10
求：建立 Simulink 仿真模型，并求出其系统在单位阶跃响应；
要求答案包括：（1）simulink 建模结构图；
x=
-2.9709
0.5491
3.6000
0.0509
2. 已知下列矩阵
2 3 1
−1 3 5
A
=
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期末大作业
1.龟兔赛跑本题旨在可视化龟兔赛跑的过程。
比赛的跑道由周长为P面积为A的矩形构成。
每单位时间，乌龟沿跑道缓慢前进一步，而兔子信心满满，每次以一个固定的概率决定走或不走。
如果选择走，就从2-10步中等概率选择一个步长。
每个单位时间用一个循环表示。
赛跑从矩形跑道左上点（0，0）开始，并沿顺时针方向进行。
不管是乌龟或兔子，谁先到达终点，比赛就告结束。
要求：编写MATLAB程序可视化上述过程。
程序以P，A以及兔子每次休息或前进的概率为输入参量。
程序必须可视化每个时刻龟兔赛跑的进程，并以红色“*”表示乌龟，蓝色的“—”表示兔子。
测试时可取P=460，A=9000。
通过上述例子，可否从理论和实验角度估计兔子休息或前进的概率，是的兔子和乌龟在概率意义下打平手。
2.黄金分割Fibonacci数列F n通过如下递推格式定义F n=F n-1+F n-2，其中F0=F1=1要求：1.计算前51项Fibonacci数，并存入一个向量2.利用上述向量计算比值F n/F n-13.验证该比值收敛到黄金比例(1+5)/23.图像处理此题旨在熟悉图像处理的基本操作，请各位自己选择一张彩色图像pMATLAB以三维数组读取一张彩图。
该彩图上每个像素位置分别存放一个取值0-255的三维向量，其三个分量分别表示该点的红（R）绿（G）蓝（B）强度信息。
要求：编写MATLAB程序，读入原始彩色图像，并且在一个图形窗口界面下显示六张图像。
这六张图分别是原始RGB彩图，及其5各变形：RBG，BRG，BGR，GBR和GRB。
每张子图要求以其对应变形命名。
最后将图像以a.jpg形式保存并黏贴至报告中。
提示：imread, imshow, cat。
MATLAB大作业
贵州大学实验报告姓名学号实验组实验时间指导教师成绩实验项目名称MATLAB大作业实验通过综合考察，让同学们能对学到的知识综合运用，更好地灵活地解决问题目的实验利用matlab提供的一些函数实现数据分析和数值计算。
原理实验运行MATLAB7.0或以上的计算机一台仪器实验内容实验数据1、试编写名为test01.m的MATLAB函数，用以计算下述的值：⎪⎩⎪⎨⎧-<->=ttntttntf的)4/sin()(si对所有)4/sin(其他情况)sin(的)4/sin()(si对所有)4/sin()(ππππ要绘制t关于函数f(t)的图形，其中t的取值范围为ππ66≤≤-t，间距为10/π。
（提示：注意要产生一系列的点，这里可考虑t的输入是向量形式，可以利用find函数找出所需限定值的元素的位置，对其按需要赋值后，再进行绘图；其次，另外一种思路，也可考虑使用循环的形式来实现）①test01.m代码如下：function y=test01(t)if sin(t)>sin(pi/4)y=sin(pi/4);elseif sin(t)<sin(-pi/4)y=sin(-pi/4);elsey=sin(t);end②调用test01函数计算值③绘制t关于函数f(t)的图形，其中t的取值范围为ππ66≤≤-t，间距为10/π代码如下：t=-6*pi:pi/10:6*pi;y=sin(t);y(find(y>sin(pi/4)))=sin(pi/4);%f(t)=sin(pi/4),y>sin(pi/4) y(find(y<sin(-pi/4)))=sin(-pi/4);%f(t)=sin(pi/4),y>sin(pi/ 4)plot(t,y) %绘制图形title('t from -6{\pi} to 6{\pi}') %添加图形标题xlabel('t') %x轴标签ylabel('f(t)') %y轴标签实验结果：2、编写函数，在同一窗口的4个子图中利用plot等语句绘制y=at2图像，其中a=[1 2 5 10],t范围[-2，5]。
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、命令窗口中输入如下命令，在%号后写出每行命令的功能。
（不用写出运算结果） A=1; f0=4; phi=pi/3; w0=2*pi*f0;t=0:0.01:1; %__________________________________________ y=A*cos(w0*t+phi) ;axis off %__________________________________________ PT=plot(t,y); %__________________________________________ set(PT,’LineWidth ’,[5]); %__________________________________________ AX=gca; %__________________________________________ set(AX,’FontSize ’,14); %__________________________________________ axis square %__________________________________________ xlabel(‘时间[s]’); %__________________________________________ T=title(‘余弦波’) %__________________________________________ set(T,’FontSize ’,’16’); %__________________________________________ legend(‘cos ’); %__________________________________________ grid on %__________________________________________二、 程序设计题2、已知两个矩阵：3765213100132610A ---⎡⎤⎢⎥---⎢⎥=⎢⎥-⎢⎥-⎣⎦，1468237021570010B ⎡⎤⎢⎥⎢⎥=⎢⎥--⎢⎥-⎣⎦将矩阵A 中所有等于-1的元素改为3，将矩阵B 中等于0的元素值改为新的矩阵A 中相应位置元素的值，最后将新的矩阵A 的第三、四列的元素全改为-1，编写相应的程序。
MATLAB大作业
1. 某控制系统为单位负反馈系统，其开环传递函数为：()()10.5KG s s s =+。
（1） 要求设计一校正装置，使校正后系统的静态速度误差系数20v K >，超调量%15%σ<，调整时间5s t s <，同时使用Simulink 仿真（示波器）验证。
（2） 绘制校正前后系统的根轨迹。
解：1）查看符合条件的zeta 。
zeta=0:0.001:0.99;sigma=exp(-zeta*pi./sqrt(1-zeta.^2))*100; plot(zeta,sigma) xlabel('\zeta'); ylabel('\sigma');title('\sigma%=e^{-\zeta*\pi/sqrt(1-\zeta^2)}*100%', 'fontsize',16) gridz=spline(sigma,zeta,15)运行结果： z=0.5169σ与ζ的关系曲线如图1-1图1-12）求系统期望主导极点sigma=0.15;zeta=((log(1/sigma))^2/((pi)^2+(log(1/sigma))^2))^(1/2);wn=1/zeta;p=[1 2*zeta*wn wn*wn];s=roots(p)运行结果：s = -1.0000 + 1.6560i，-1.0000 - 1.6560i3）查看根轨迹G0=tf([1],[0.5 1 0]);rlocus(G0);sgrid(.5169,[])根轨迹如图1-2和局部放大根轨迹图1-3图1-2图1-3读出增益为1.88，则校正前系统的稳态误差Kv=1.88/1=1.88。
按要求，偶极子的零点和极值比值应为20/1.88=10.6383,取Zc=0.01，Pc=0.01/11=0.0009,因此校正环节为0009.0)01.0(88.1)(++=s s s Gc4)得出校正后的系统，并进行验证。
matlab期末大作业题目及答案
matlab 期末大作业（30分，每题6分）1. 积分运算（第四数值和五章符号）(1)定积分运算：分别采用数值法（quad ，dblquad ）和符号运算（syms, int ）一重定积分π⎰1. 数值法（quad ）a) 运行代码：b) 运行结果：2. 符号运算（syms ）a) 运行代码：b) 运行结果：二重定积分112200()x y dxdy+⎰⎰1.数值法（dblquad）：a)运行代码：b)运行结果：2.符号运算（syms）：a)运行代码：b)运行结果：(2) 不定积分运算sin dxdy ⎰⎰（(x/a)+b/y） i.运行代码：ii.运行结果：2. 用符号法和数值法求解线性代数方程 （第五章和第二章）⎩⎨⎧=+=+12*22x *213*12x *a11y a a y a （1） 用syms 定义待解符号变量x,y 和符号参数a11,a12,a21,a22，用符号solve 求x,y 通解 1. 运行代码：2. 运行结果：（2） 用subs 带入a11=2,a12=4,a21=6,a22=8,求x 和y 特解，用vpa 输出有效数值4位的结果 1. 运行代码：2. 运行结果：（3） 采用左除(\)和逆乘法求解符号参数赋值后的方程 ⎩⎨⎧=+=+12*8x *63*4x *2y y1. 运行代码：2. 运行结果：3.数值法和符号法求解非线性方程组(第四数值和五章符号 )（1）采用数值法（fsolve ）求解初始估计值为x0 = [-5; -5]的数值解1. 运行代码：2. 运行结果：21x 21x 21e x 2x e x x 2--=+-=-（2）符号法（solve ）的符号结果用eval 或double 转化为数值结果.1. 运行代码：2. 运行结果：4. 解二阶微分方程 (第四数值和五章符号 )⎪⎩⎪⎨⎧===++6)0(',0)0(09322y y y dx dy dx y d（1）数值ode 求特解，用plot (x,y) 画t 在[0,10]范围内（x ，y ）数值曲线 1. 运行代码：2. 运行结果：（2）符号运算dsolve求通解，用ezplot画t在[0,10]范围内（x，y）符号曲线1. 运行代码：2. 运行结果：5. 三维绘图（第六章）已知：x和y都在[-8,8]范围内，采用subplot(3,1,x)绘制三个子图，它们分别是用meshgrid和mesh绘制网格图、用c=contour 绘制等位线和用surf 绘制曲面图1.运行代码：2.运行结果：。
MATLAB大作业
若系统是完全可控的，则可引入状态反馈调节器，且：
u v Kx
这时，闭环系统的状态空间模型为： (A BK)x Bu x y Cx 系统方框图如下：
此时，相稳定裕量Pm=46.7deg,满足题目中430<γ<480的要求。
（5）计算校正后系统阶跃响应曲线及其性能指标。 程序如下： % BodeJiaoZheng4 % 系统校正后性能指标及阶跃响应曲线 global y t; k0=1000;n1=1;d1=conv([1 0],[1 2]); s1=tf(k0*n1,d1); n2=[0.05443 1];d2=[0.009339 1]; s2=tf(n2,d2); sope=s1*s2; sys=feedback(sope,1); step(sys); %绘制阶跃响应曲线 [y,t]=step(sys); % 求出阶跃响应的函数值及其对应时间 [sigma,tp,ts]=ste(y,t) % 调用函数ste()
1 s(s 2)
2）作原系统的波特图与阶跃响应曲线，检查是否满足要求。
在MATLAB命令窗口输入：（SysExpBode.M) >>k0=1000;n1=1;d1=conv([1 0],[1 2]); [mag,phase w]=bode(k0*n1,d1); s1=tf(k0*n1,d1); figure(1);
根据校正后的结构与参数，调用函数ste()，给出上述例题系统的性 能指标及阶跃响应曲线。程序如下：
%计算系统超调量，峰值时间，调节时间。 global y,t s1=tf(20700,[1 20 75 0]); s2=tf([1 3.841],[1 73.12]); sope=s1*s2;sys=feedback(sope,1); step(sys) [y,t]=step(sys); [sigma,tp,ts]=ste(y,t) 运行该程序，可得到系统的阶跃响应曲线，并有系统性能指标： 超调量：sigma=0.2489<30% 峰值时间：tp=0.2209 调节时间：ts=0.3037<0.5s
matlab大作业例子
Matlab作业题目1:(1)程序部分：>> x=0:pi/50:2*pi;>> y=cos(0.5+((3*sin(x))./(1+x.^2))); >> plot(x,y)>> xlabel('x');>> ylabel('y');>> title('x-y');(2)运行结果截图：（1）程序部分：>> clear>> clc>> x=0:pi/100:4*pi;>> y1=sin(x);>> y2=cos(x);>> plot(x,y1,'r-',x,y2,'g:') >> hold on>> y3=x(find(abs(y1-y2)<0.001)); >> plot(y3,sin(y3),'*')（2）运行结果截图：（1）程序部分：>> t=(0:2*pi/100:2*pi)'; >> y1=sin(t)*[1,-1]; >> y2=sin(t).*sin(9*t); >> subplot(2,1,1);>> plot(t,[y1,y2]);>> subplot(2,1,2);>> plot(t,[y1,y2]) （2）运行结果截图：题目4（1）程序部分：>> t=0:pi/50:4*pi;>> y=exp(-t/3);>> y0=exp(-t/3).*sin(3*t); >> plot(t,y,'r-',t,y0,'b:') >> xlabel('\bf\it t')>> ylabel('\bf\it y')>> grid on;（2）运行结果截图：题目5（1）程序部分：>> n=0;>> sum=0;>> while sum<2000n=n+1;sum=sum+n;end>> n（2）运行结果截图：题目6（1）程序部分：for j=1:3n=input('n='); sum=0;for i=1:na=1/(i.^2);sum=sum+a;endPI=sqrt(6*sum) end（2）运行结果截图：题目7（1）程序部分：n0=0;y0=0;while 3*y0<5n0=n0+1;y0=y0+1/(2*n0-1);y=y0-1/(2*n0-1);n=n0-1;endn（2）运行结果截图：题目8（1）程序部分：for i=1:7x=input('put x:');if x<0&x~-3y=x^2+x-6;elseif x>=0&x<5&x~=3&x~=2 y=x*x-5*x+6;elsey=x*x-x-1;endend（2）运行结果截图：题目9①switch语句（1）程序部分：>> score=input('请输入成绩'); switch fix(score/10)case {9,10}disp('A');case {8}disp('B');case {7}disp('C');case {6}disp('D');case {0,5}disp('E');otherwisedisp('error');end（2）运行结果截图：如右图②if语句（1）程序部分：score=input('请输入成绩:') n=score/10;if n>=9&n<10disp Aelseif n>=8&n<9disp Belseif n>=7&n<8disp Celseif n>=6&n<7disp Delseif n>=0&n<6disp Eelsedisp errorend（2）运行结果截图：题目10（1）程序部分：t=input('员工的工作时间:')if t<60m=84*t-700;elseif t>120m=84*(t-120)*1.15+84*120;elsem=84*t;endm2）运行结果截图：题目11（1）程序部分：a=round(20*rand(5,6))n=input('请输入n的值:') tryb=a(n,:);catchb=a(5,:);endblasterr（2）运行结果截图：。
matlab大作业1
计算机实习班级：姓名：学号：题目一：使用while循环及for循环计算1+2+3+…+100分析：无。
解答：MATLAB程序：%使用for-end语句>> clear>> sum=0;>> for i=1:1:100sum=sum+i;end>> sumsum =5050%使用while-end语句>> clear>> sum=0;i=0;>> while i<100i=i+1;sum=sum+i;end>> sumsum =5050小结：了解matlab中的基本语句，与C语言中的相关内容进行对比，从而更好的掌握及运用语句。
题目二：圆轮在平面上滚动时，圆周上任一点的轨迹为摆线，圆内任一点的轨迹为内摆线，圆外任一点的轨迹为外摆线，建立模型，并作出图。
分析：建立模型：设r为轮的半径，R为点的半径，则由分析得方程x=rt-Rsinty=r-Rcost；（t为参数）MATLAB程序：t=0:0.1:10;r=1;R=input('R=');x=r*t-R*sin(t); y=r-R*cos(t); plot(x,y),axis('equal');hold ont=0:0.1:10;r=1;R=input('R=');x=r*t-R*sin(t); y=r-R*cos(t); plot(x,y),axis('equal');hold ont=0:0.1:10;r=1;R=input('R=');x=r*t-R*sin(t); y=r-R*cos(t); plot(x,y),axis('equal');hold on令R=1, 0.5, 2，作出如下图形小结：逐渐将matlab应用到解决实际问题当中，了解matlab强大的绘图功能。
MATLAB 大作业
MATLAB 大作业请各位同学，自己完成matlab 大作业的内容。
禁止相互抄袭，如有雷同，零分计算。
大作业的格式按照实验报告的格式书写，务必标明题号，作业完成后，将生成的报告打印出来提交。
正文的字号以宋体五号字，1.5倍行距的格式打印。
请与18周五前将大作业报告由班级负责人统一收齐交给我，过期不收。
1、 试编写名为test01.m 的MATLAB 函数，用以计算下述的值：⎪⎩⎪⎨⎧-<->=t t n t t t n t f 的对所有其他情况的对所有)4/sin()(si )4/sin()sin()4/sin()(si )4/sin()(ππππ要绘制t 关于函数f (t )的图形，其中t 的取值范围为ππ66≤≤-t ，间距为10/π。
（提示：注意要产生一系列的点，这里可考虑t 的输入是向量形式，可以利用find 函数找出所需限定值的元素的位置，对其按需要赋值后，再进行绘图；其次，另外一种思路，也可考虑使用循环的形式来实现）2、 编写函数，在同一窗口的4个子图中利用plot 等语句绘制y=at 2图像，其中a=[1 2 5 10],t 范围[-2，5]。
3、 求函数32)(3-+=x x x f 在区间[-5,5]上的最大值和最小值。
4、 求解函数⎰102dx e x 的数值积分和符号积分，并比较结果。
5、 求解微分方程3|;1|2)1(002='='=''+==x x y y y x y x 的精确解和解析解，并绘制图形。
假设求解区间为[0,10] 。
6、 说说你对MATLAB 及应用这门课程学习后的体会，另外请说明在所学章节中哪一章的内容你最感兴趣，为什么？哪一章的内容你认为是没有必要学习的，为什么？如果可以选择MATLAB 的学习的内容的话，谈谈你所期望学到的知识类别的前三种。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
MATLAB大作业作业要求：（1）编写程序并上机实现，提交作业文档，包括打印稿（不含源程序）和电子稿（包含源程序），以班为单位交，作业提交截止时间6月24日。
（2）作业文档内容：问题描述、问题求解算法（方案）、MATLAB程序、结果分析、本课程学习体会、列出主要的参考文献。
打印稿不要求MATLAB程序，但电子稿要包含MATLAB 程序。
（3）作业文档字数不限，但要求写实，写出自己的理解、收获和体会，有话则长，无话则短。
不要抄袭复制，可以参考网上、文献资料的内容，但要理解，要变成自己的语言，按自己的思路组织内容。
（4）从给出的问题中至少选择一题(多做不限，但必须独立完成，严禁抄袭)。
（5）大作业占过程考核的20%，从完成情况、工作量、作业文档方面评分。
第一类：绘制图形。
（B级）问题一：斐波那契（Fibonacci）螺旋线，也称黄金螺旋线（Golden spiral），是根据斐波那契数列画出来的螺旋曲线，自然界中存在许多斐波那契螺旋线的图案，是自然界最完美的经典黄金比例。
斐波那契螺旋线，以斐波那契数为边的正方形拼成的长方形，然后在正方形里面画一个90度的扇形，连起来的弧线就是斐波那契螺旋线，如图所示。
问题二：绘制谢尔宾斯基三角形（Sierpinskitriangle）是一种分形，由波兰数学家谢尔宾斯基在1915年提出,它是一种典型的自相似集。
其生成过程为：取一个实心的三角形（通常使用等边三角形），沿三边中点的连线，将它分成四个小三角形，然后去掉中间的那一个小三角形。
接下来对其余三个小三角形重复上述操作，如图所示。
问题三：其他分形曲线或图形。
分形曲线还有很多，教材介绍了科赫曲线，其他还有皮亚诺曲线、分形树、康托（G. Cantor）三分集、Julia集、曼德布罗集合（Mandelbrot set），等等。
这方面的资料很多（如/content/16/0103/14/5315_525141100.shtml），请分析构图原理并用MATLAB实现。
问题四：模拟掷骰子游戏：掷1000次骰子，统计骰子各个点出现的次数，将结果以下表的形式显示，并绘制出直方图。
点数 1 2 3 4 5 6出现次数166 150 164 162 184 174问题五：利用MATLAB软件绘制一朵鲜花，实现一定的仿真效果。
提示：二维/三维绘图，对花瓣、花蕊、叶片、花杆等的形状和颜色进行详细设置。
第二类：插值与拟合。
（B级）问题一：有人对汽车进行了一次实验，具体过程是，在行驶过程中先加速，然后再保持匀速行驶一段时间，接着再加速，然后再保持匀速，如此交替。
注意，整个实验过程中从未（1）分别使用最近点插值、线性插值、三次埃尔米特插值和三次样条插值进行计算[0,110]时间段50个时间点的速度。
（2）绘制插值图形并标注样本点。
问题二：估算矩形平板各个位置的温度。
已知平板长为5m，宽为3m，平板上3×5栅格点上的温度值为44，25，20，24，30；42，21，20，23，38；25，23，19，27，40。
（1）分别使用最近点插值、线性插值和三次样条插值进行计算。
（2）用杆图标注样本点。
（3）绘制平板温度分布图。
问题三：自行车道的设计。
对9条道路上的自行车道宽度以及自行车与过往机动车之间（1）对数据进行线性拟合。
（2）绘制拟合曲线和样本点。
（3）如果自行车与过往机动车之间安全距离的最小距离是1.8m，试计算相应的自行车道宽度的最小值。
问题四：在水资源工程学中，水库的大小与为了蓄水而拦截的河道中的水流速度密切相关。
对于某些河流来说，这种长时间的历史水流记录很难获得。
然而通常容易得到过去若干年间关于降水量的气象资料。
鉴于此，推导出流速与降水量之间的关系式往往特别有用。
只要获得那些年份的降水量数据，就可以利用这个关系式计算出水流速度。
下表是在被水库拦截的某河道中测得的数据。
降水量（cm）88.9 108.5 104.1 139.7 127 94 116.8 99.1流速（m3/s）14.6 16.7 15.3 23.2 19.5 16.1 18.1 16.6（1）对数据进行线性拟合。
（2）绘制拟合曲线和样本点。
（3）如果某年的降水量是120cm，利用拟合直线估算当年的水流速度。
（4）若流域面积为1100km2，估计在其他过程中，如蒸发、深层地下水渗透和消耗用途，损失的降水量占总体降水量的比例。
问题五：假设有已知实测数据如下表所示：x 0.1 0.2 0.3 0.4 0.5 0.6 0.7 0.8 0.9 1.0y 2.3201 2.6470 2.9707 3.2885 3.6008 3.9090 4.2147 4.5191 4.8232 5.1275假设已知该数据可能满足的原型函数为，试求出满足数据的最小二乘解a,b,c,d的值。
提示：曲线拟合并绘图分析第三类：定积分问题。
（B级）问题一：地球密度随着离中心（r=0）距离的变化而变化，不同半径处的密度如表所示，r（km）0 1100 1500 2450 3400 3630 4500 5380 6060 6280 6380 ρ（g/cm3）13 12.4 12 11.2 9.7 5.7 5.2 4.7 3.6 3.4 3 问题二：河道平均流量Q（m3/s）可使用速度和深度的乘积的积分来计算（河道横截面不规则），公式如下。
其中V(x)是离岸x（m）距离处的水速（m/s）,H(x)是离岸x距离处的水深（m）。
根据x 0 1.6 4.1 4.8 6.1 6.8 9V 0 0.08 0.61 0.68 0.55 0.42 0x 0 1.1 2.8 4.6 6 8.1 9H 0 0.21 0.78 1.87 1.44 1.28 0.2第四类：线性方程组求解。
（B级）问题一：多项式插值指的是采用唯一的n-1次多项式对n个数据点进行拟合。
该多项式的一般形式为：p(x)=p1x n-1+p2x n-2+…+p n-1x+p n确定这些系数的一种直接方法是，建立n个线性代数方程，然后求解。
已知一个四次多项式通过x 200 250 300 400 500y 0.746 0.675 0.616 0.525 0.457 （1）建立线性方程组，并求解得到多项式的系数。
（2）计算该线性方程组系数矩阵的条件数，并进行解释。
（3）绘制多项式曲线并求其零点。
问题二：如图所示，5个反应器通过导管连接在一起。
每根导管中化学物的传输率等于流速（Q，单位是m3/s）乘以化学物浓度（c，单位是mg/m3）。
若系统达到稳定状态，流入和流出每个反应器的质量相等。
例如，对于第一个反应器来说，质量守恒可表示为：Q01c01+Q31c3=Q15c1+Q12c1（1）使用LU分解计算平衡方程系数矩阵的逆矩阵。
（2）求各反应器中化学物的稳态浓度。
问题三：静定桁架受力分析。
（1）如图所示，求力和反作用力。
（2）求受力平衡方程系数矩阵的逆矩阵，对于逆矩阵第二行中的零，作何解释。
（3）将节点1的力改为方向向上，计算这种改变对H2和V2的影响。
（4）将节点1的力撤销，而在节点1和2处施加1500N的水平外力，求节点3处垂直反作用力（V3）。
第五类：一元方程求解。
（B级）问题一：在热力学中，下列多项式将干燥空气的零压力比热c p（单位为kJ/（kgK））与温度（单位为K）关联起来了：C p=0.99403+1.671×10-4T+9.7215×10-8T2-9.5838×10-11T3+1.9520×10-14T4（1）绘制在T=0~1200K范围内，c p随温度变化的曲线。
（2）求对应于1.1kJ/（kgK）比热的温度。
问题二：在化学工程中，将水蒸汽（H2O）加热到足够高的温度，使得大部分水发生分解或分离而形成氧气（O2）和氢气（H2）：H2O H2+O2如果假定其中只存在这一种化学反应，那么已经发生分解的H2O所占比列x可以表示为：其中K为该反应的平衡系数，P t为混合物的总压强。
如果P t=4，且K=0.05，那么求满足该式子的x值。
第六类：最优化问题。
（B级）问题一：最大利润问题。
某公司生产两种产品的产量分别为x，ykg，其相应的成本满足以下函数：C(x,y)=x2+2xy+2y2+2000已知产品x的价格为200元/kg，产品y的价格为300元/kg，并假定两种产品全部售完，试求使公司获得最大利润产品产量以及公司的最大利益润。
问题二：作用在螺旋桨上的总阻力可以通过下式估计:摩擦力升力其中，D=阻力，=飞行高度与海平面之间的大气密度比(ratioofairdensity)，W=重量，V=速度。
如图所示，当速度增加时，对阻力的两个部分受到的影响是不同的。
摩擦阻力随速度的增加而增加，但由升力引起的阻力却随速度的增加而下降。
二者的结合导致一个最小的阻力。
（1）如果=0.6、W=16000，求最小阻力及阻力最小时的速度值。
（2）进行敏感性分析以确定当W为12000~20000的过程中，最优值是如何变化的，取=0.6。
螺旋桨上阻力与速度的关系图问题三：如图所示，一个梯子通过支撑角分别与两个面接触，梯子的最大可能长度可以通过计算下面函数取值最小时的值而确定。
对于==2m的情况，绘制L随变化的图形，其中的取值范围为45°~135°。
通过一个墙角连接两个墙面的梯子问题四：对于一架稳定水平航行的喷气机，推力与阻力平衡，升力与重力平衡(如图所示)。
在这种情况下，当阻力与速度的比例最小时，会出现最佳巡航速度。
阻力可以用下式计算:其中是零升力时的阻力系数，是升力系数，AR是展弦比。
在稳定水平飞行的情况下，升力系数可以用下式计算:其中W是喷气机重量（N），是空气密度（kg/）,是速度(m/s), A是机翼平面面积()，然后阻力可以用下式计算:在稳定水平飞行中，喷气机受到的四个主要力使用这些公式，确定在海平面上10千米飞行的670kN喷气机的最佳稳定巡航速度。
在计算中应用以下参数：A=150，AR=6.5，=0.018， kg/。
问题五：某公司经营两种设备，第一种设备每件售价29元，第二种设备每件售价455元。
根据统计，售出第一种设备一件所需的营业时间平均为0.5小时，第二种设备是（2+0.25*X2）,其中X2是第二种设备的售出数量。
已知该公司在这段时间内的总营业时间为800小时，试确定使营业额最大的营业谋划。
提示：即两种设备各准备售出多少件，使得在规定的营业时间里营业额最大。
这是一个有约束的最优化问题求解，考虑用fmincon函数问题六：某车间生产A、B两种产品，已知生产产品A、B需要原料分别为3公斤和4公斤，所需的工时分别为5分钟和3分钟，现在可以应用的原料为120公斤，工时为150分钟，每生产一件A和B分别可获得7元和5元的利润，应当如何安排生产A、B的件数，才能使车间获得最大利润？提示：线性规划问题，考虑用linprog函数问题七：某种作物在全部生产过程中至少需要32公斤氮，磷以24公斤为宜，钾不得超过42公斤。