MATLAB结课作业(2016.12)
Matlab作业习题与答案详解(附程序)

Matlab作业习题与答案详解(附程序)clear all;clc;close all;x=-10:0.01:20;y=4*sin(x)./x;ymin=min(y)二、蒙特卡罗算法的数值计算当前的油位高度是2.3米,见图1。
模拟油流进储油罐的过程(图维数任选),请计算罐内油量。
三维的效果图参见图2。
储油罐由两部分组成,中间是圆柱体,两端是球罐体。
(本题简化自2011年UCMCM A题《储油罐的变位识别与罐容表标定》,细节参见原题原题附件2 cumcm2010A.doc。
)图1图2主程序:clc;clear all;close all;center1=[-3.375,0,1.5]; %左球罐中心center2=[3.375,0,1.5]; %右球罐中心n=10000; %每次的撒点数delta=0.02; %层高h=3;en=h/delta;Show; %画出油罐for i=0:en-1x=(rand(1,n)-0.5)*10; %随机生成n个点y=(rand(1,n)-0.5)*h;z=(rand(1,n)*delta+i*delta);Z=[x;y;z];[dis1dis2]=juli(center1,center2,Z); %算出各点对应的距离index=find(((x>-4&x<4)&dis2<1.5)|(x<-4|x>4)&dis1<1.625); %找出在罐内的点plot3(x(index),y(index),z(index),'.k'); %画出在罐内的点drawnowend子程序1:function [dis1 dis2]=juli(a,b,q)d11=q(1,:)-a(1);d12=q(2,:)-a(2);d13=q(3,:)-a(3);d1=sqrt(d11.^2+d12.^2+d13.^2);d21=q(1,:)-b(1);d22=q(2,:)-b(2);d23=q(3,:)-b(3);d2=sqrt(d21.^2+d22.^2+d23.^2);d1(d1>d2)=d2(d1>d2);dis1=d1;dis2=sqrt(d12.^2+d13.^2);子程序2:function tu=Show%===圆柱部分==figure('color','w')h=3;y='3/2*cos(s)';z='3/2*sin(s)+1.5';x='t';ezmesh(x,y,z,[0 pi*2 -4 4]);axis equalhidden offhold on%===左罐部分==z='cos(s)*(2.6406-t.^2).^(1/2)+1.5'; y='sin(s)*(2.6406-t.^2).^(1/2)';x='t-3.375';ezmesh(x,y,z,[0 pi*2 -1.625 -0.625]);axis equalhidden offhold on%===右罐部分===z='cos(s)*(2.6406-t.^2).^(1/2)+1.5';y='sin(s)*(2.6406-t.^2).^(1/2)';x='t+3.375';ezmesh(x,y,z,[0 pi*2 0.625 1.625]);colormap(gray)axis equalhold onaxis offhidden off三、元胞自动机的简单应用1.简单交通流模拟请模拟一个4车道的交通流,车辆密度为0.3,其余规则自行定义。
MATLAB 结业解答

MATLAB 期末结课大作业(必做题) 要求:抄题、写出操作命令、运行结果,并根据需要,贴上运行图。
1、已知函数x x e x f x sin cos )(-=,作出函数的大致图像。
>> fplot('exp(x)*cos(x)-sin(x)',[-10,10]) >> fplot('exp(x)*cos(x)-sin(x)',[-5,5])2、求下列极限:(1)30sin lim x x x x -→ (2) xx x ⎪⎭⎫ ⎝⎛+∞→11lim (1)>> clear %清除>> syms x; %声明符号变量x>> s1=(x-sin(x))/x^3; %建立符号函数s1>> limit (s1,x,0) % 对是关于变量x 在[x ,0)上求值 ans =1/6(2)>> s2=(1+1/x)^x; %建立符号函数s2>> limit (s2,x,inf) % 对是关于变量x 在[x ,+∞)上求值 ans =exp(1)3、求下列函数的导数:(1)x e y x sin = (2) x e x x y 22sin +=(1)>> clear %清除>> syms x; %创建符号变量x>> y=exp(x)*sin(x); %创建符号函数y>> diff(y) %对y 值进行求导ans =exp(x)*sin(x)+exp(x)*cos(x)(2)>> clear %清除>> syms x; %创建符号变量x>> y=sin(x)+x^2*exp(2*x); %创建符号函数>> diff(y) %对y 进行求导ans =cos(x)+2*x*exp(2*x)+2*x^2*exp(2*x)4、求.)1(532⎰-dx x x 和.sin ⎰xdx e x(1)>> clear>> syms x %创建符号变量x>> s1=x^2*(1-x^3)^5; %创建符号函数s1>> int (s1) %对s1进行求导数ans =-1/18*x^18+1/3*x^15-5/6*x^12+10/9*x^9-5/6*x^6+1/3*x^3(2)>> s2=exp(x)*sin(x); %创建符号函数s2>> int(s2) %对s2进行求导数ans =-1/2*exp(x)*cos(x)+1/2*exp(x)*sin(x)5、求.)(102⎰-dx x x 和.1102⎰-dx x x (1)>> clear %清除>> syms x %定义符号变量x>> s1=x-x^2; %定义符号函数s1>> int(s1,0,1) %对s1关于变量x 在[0,1)上求定积分ans =1/6(2)>> clear %清除>> syms x %定义符号变量x>> s2=x*sqrt(1-x^2);%定义符号函数s1>> int(s2,0,1) %对s2关于变量x 在[0,1)上求定积分ans =1/36、已知二元函数),(cos )sin(2xy xy z +=试求y x z x z y z x z ∂∂∂∂∂∂∂∂∂222,,,。
matlab结业作业

《Matlab/Simulink 与控制系统仿真》结业作业一、必做题Matlab学习心得体会:本学期我们新开了一门MATLAB/Simulink与控制系统仿真的课程,虽说是这个学期才开始学习,但是我们却对MATLAB这个软件的基本使用一点都不陌生。
而本学期的专门学习,更让我获益匪浅。
第一部分是利用MATLAB函数来对控制系统进行分析。
①例如要求取一个系统的阶跃响应时,我就可以调用MATLAB中的step函数②MATLAB在稳定性分析应用时,通过直接求根函数roots来判定系统的稳定性③在利用MATLAB分析系统的根轨迹时,使用rlocus函数来求取根轨迹图,使用pzmap函数求取系统的传递函数零极点④计算给定一组根的根轨迹增益函数 [k,poles]=rlocufind(sys);[k,poles]=rlocfind(sys,p)⑤MATLAB在稳定性分析应用时,也常通过使用用于计算系统稳定裕度的margin函数来完成,它可以从频率响应数据中计算出幅值裕度、相角裕度以及对应的频率⑥对于使用MATLAB来对系统进行频域响应的分析,又往往提供了绘制函数nyguist和bode来绘制系统的奈奎斯特曲线图和伯德图⑦其他:grid%添加栅格 hold on%保持,继续在原图绘画 title%添加标题legend(x,y)%添加文字标注 gyext('TEXT')%在任意可选位置添加文字信息 xlabel(x,y);ylabel(x,y)%添加坐标轴分割图形窗口 subplot(1,2,1)%分割为一行两列,并即将绘制第一个第二部分是利用MATLAB的Simulink模块来对控制系统进行建模仿真。
Simulink是一个对动态系统(包括连续系统、离散系统和混合系统)进行建模、仿真和综合分析的集成软件包,是MATLAB的一个附加组件,其特点是模块化操作、易学易用,而且能够使用MATLAB提供的丰富的仿真资源。
MATLAB作业2_2016_ 作业本

MATLAB 作业二1、 试求出如下极限。
222222001cos()lim ()x y x y x y x y e +→→-++解:>>syms x y;>>f=(1-cos(x^2+y^2))/(x^2+y^2)/exp(x^2+y^2);>> L=limit(limit(f,y,0),x,0)L =2、 试求出下面函数的导数。
22atan ln()y x y x=+ 解:>> syms x y;>>f=atan(y/x)-log(x^2+y^2);>> pretty(-simple(diff(f,y)/diff(f,x)))x - 2 y-------y + 2 x3、假设1cos u -=,试验证22u u x y y x ∂∂=∂∂∂∂。
证明:二者相等亦可以由二者之差为零来证明,故由下面的语句直接证明。
>> syms x y; u=acos(x/y);diff(diff(u,x),y)-diff(diff(u,y),x)ans =4、 假设20(,)xy t f x y e dt -=⎰,试求222222x f f f y x x y y ∂∂∂-+∂∂∂∂。
解: 由下面的命令可以得出所需结果。
>> syms x y tf=int(exp(-t^2),t,0,x*y);x/y*diff(f,x,2)-2*diff(diff(f,x),y)+diff(f,y,2)simple(ans)ans =-2*exp(-x^2*y^2)*(-x^2*y^2+1+x^3*y)5、 假设已知函数矩阵323(,,)sin y x e z f x y z x y z ⎡⎤+=⎢⎥+⎣⎦,试求出其Jacobi 矩阵。
解:Jacobi 矩阵可以由下面的语句直接得出。
>> syms x y zF=[3*x+exp(y)*z; x^3+y^2*sin(z)];jacobian(F,[x,y,z])ans =[ 3, exp(y)*z, exp(y)][ 3*x^2, 2*y*sin(z), y^2*cos(z)]6、 试求解下面的不定积分问题。
matlab2016年课后习题及答案详解【精选】

2012年12月7日星期五第2章 MATLAB 矩阵运算基础2.1 在MATLAB 中如何建立矩阵,并将其赋予变量a ?⎥⎦⎤⎢⎣⎡194375>> a=[5 7 3;4 9 1]2.2 有几种建立矩阵的方法?各有什么优点?可以用四种方法建立矩阵:①直接输入法,如a=[2 5 7 3],优点是输入方法方便简捷;②通过M 文件建立矩阵,该方法适用于建立尺寸较大的矩阵,并且易于修改;③由函数建立,如y=sin(x),可以由MATLAB 的内部函数建立一些特殊矩阵;④通过数据文件建立,该方法可以调用由其他软件产生数据。
2.3 在进行算术运算时,数组运算和矩阵运算各有什么要求?进行数组运算的两个数组必须有相同的尺寸。
进行矩阵运算的两个矩阵必须满足矩阵运算规则,如矩阵a 与b 相乘(a*b )时必须满足a 的列数等于b 的行数。
2.4 数组运算和矩阵运算的运算符有什么区别?在加、减运算时数组运算与矩阵运算的运算符相同,乘、除和乘方运算时,在矩阵运算的运算符前加一个点即为数组运算,如a*b 为矩阵乘,a.*b 为数组乘。
2.5 计算矩阵与之和。
⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡897473535⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡638976242>> a=[5 3 5;3 7 4;7 9 8];>> b=[2 4 2;6 7 9;8 3 6];>> a+bans =777914131512142.6 求的共轭转置。
⎥⎦⎤⎢⎣⎡+-+-+-+-++=i 44i 93i 49i 67i 23i 57i 41i 72i 53i 84x >> x=[4+8i 3+5i 2-7i 1+4i 7-5i;3+2i 7-6i 9+4i 3-9i 4+4i];>> x’4.0000 - 8.0000i 3.0000 - 2.0000i 3.0000 -5.0000i 7.0000 +6.0000i 2.0000 +7.0000i 9.0000 - 4.0000i 1.0000 - 4.0000i 3.0000 + 9.0000i 7.0000 + 5.0000i 4.0000 - 4.0000i 2.7 计算与的数组乘积。
matlab结课作业

学校代码: 10128《MatLab在工科数学中的应用》结课作业学生姓名:学号:学院:专业:班级:指导教师:一、矩阵运算方面的题(参照第一章练习题3)生成矩阵A,通过修改矩阵A的第2行第3列元素为7,把矩阵扩充为4行5列,其中第4行第3列元素为6,其余列为0,删除第2列元素的操作生成矩阵B, 按行的逆顺序取A的1,2,3,5列,且第四行为自然数1到4,生成新的矩阵C;计算B*C,B.* C, B+C,2*B, B和C 的行列式,B\C, B/C,并找出A中大于3的元素,且将其替换为1。
A = 1 2 3 4 56 4 9 8 21 4 3 8 5上机实验步骤输入>>A=[1,2,3,4,5;6,4,9,8,2;1,4,3,8,5];B=A;>>B(2,3)=7;B(4,3)=6;B(:,2)=[]运行结果如下:输入>>C=A(3:-1:1,[1,2,3,5]) 运行结果如下:输入>>C(4,:)=1:4运行结果如下:输入>>B*C运行结果如下:输入>> B.*C运行结果如下:输入>> B+C运行结果如下:输入 >> 2*B运行结果如下:输入>> B\C 运行结果如下:输入>> det(C)运行结果如下:输入>> det(B) 运行结果如下:输入 >> B/C运行结果如下:二、绘图方面的题,包括二维和三维图形(参照第二章练习题)(1)、二维图形的绘制用ezplot 绘制函数5e xy −3sin (x +2y )=0在[-5,5]上的图形。
上机操作: 在MatLab 命令窗口输入:ezplot('5*exp(x*y)-3*sin(x*2*y)',[-5,5,-12,12]) % 运行结果如下:(2)、三维图形的绘制画出z =√x 2+y2)√x 2+y2所表示的三维曲面,x,y 的取值范围是[-10,10]。
MATLAB课程作业练习题

内蒙古工业大学电力学院《MATLAB在电力系统中的应用》练习题汇总编者:顾秀芳单位:电力学院电力工程系时间:2012年8月实验一:《MATLAB操作环境及基本操作》实验指导书1、实验目的:了解MATLAB的运行环境,界面功能;熟悉MATLAB编程的基本方法;掌握MATLAB的基本画图功能;熟悉Simulink模块的基本操作及子系统的封装方法。
2、实验内容:操作MATLAB界面的主要菜单并了解其中的功能;练习矩阵的运算;编写或读懂并运行MATLAB程序;练习并实践MATLAB的画图功能的实现方法。
练习Simulink模块的基本操作;熟悉并练习子系统的封装方法;3、实验要求:通过预习,了解MATLAB的界面的组成,各界面实现的主要功能;了解MATLAB实现运算、逻辑判断的基本语法、函数的含义等;掌握MATLAB基本绘图功能,并实现简单的绘图;熟悉Simulink的及子系统封装的主要功能及基本操作。
4、实验步骤:(1)熟悉MATLAB界面的基本功能;(2)熟悉编程基础,并读通指定程序的含义,运行结果,并分析;(3)熟悉MATLAB的基本绘图语句,读通指定并读通指定程序的含义,运行结果并分析;实现给定的绘图任务;(4)熟悉Simulink模块的操作,并完成给定模型的搭建;(5)熟悉Simulink子系统模型搭建的步骤,并完成给定子系统的搭建及参数设置。
6 、实验题目:1)MATLAB环境的了解包括:窗口菜单的主要功能;各窗口实现的功能,包括:命令窗口,workspace/history,command history三个窗口的应用。
2)MATLAB编程基础内容包括:矩阵的形成及运算(矩阵相乘、点乘);逻辑运算;函数的应用。
练习内容一:矩阵的形成及操作;1:形成矩阵2:从矩阵中提取子矩阵:3:矩阵A中元素A41的值;4:修改矩阵A的元素,使A41=3;5:矩阵A中最后2行和最后3行交汇的子矩阵的值;练习内容二:矩阵的运算1:矩阵定义为:求出下列矩阵运算的结果,并解释部分操作无法实现的原因:(1) R=A+B;(2) R=A*D;(3) R=A.*D;(4) R=A*C;(5) R=A.*C;(6) R=A\B;(7) R=A.\B;(8) R=A.^B.1.网络如图,图中标注为支路导纳yij.试通过程序完成Y导纳矩阵的形成,j含对地支路练习内容三:编程语句的理解1:运行书P20中例题,比较标量循环和矩阵运算的差异;比较程序:1.一个窗口画一个图t=0:0.1:2*pi;y1=sin(t);plot(t,y1,'r:')axis([0,2*pi,-2,2]);xlabel('huduzhi');ylabel('hanshuzhi');title('正弦函数和余弦函数');grid ongtext({'y1=sin(t)'});2.一个窗口绘两个图(1)一条plot语句画t=0:0.1:2*pi;y1=sin(t); y2=cos(t);plot(t,y1,'r:', t,y2,'b-.')axis([0,2*pi,-2,2]);xlabel('t');ylabel('y1,y2');title('正弦函数和余弦函数');grid ongtext({'y1=sin(t)'; 'y2=cos(t)'});(1)分多条plot语句画t=0:0.1:2*pi;y1=sin(t); y2=cos(t); y3=tan(t) plot(t,y1,'r-')holdplot(t,y2,'b--')plot(t,y3,'k-.')axis([0,2*pi,-2,2]);xlabel('t');ylabel('y1,y2,y3');grid ongtext({'y1=sin(t)'; 'y2=sin(t)';'y3=sin(t)'}); legend('y1','y2','y3')3、一个窗口分割成三个小窗口t=0:0.1:2*pi;y1=sin(t);subplot(3,1,1)plot(t,y1,'r:')axis([0,2*pi,-2,2]);xlabel('t');ylabel('y1');title('正弦函数');grid ongtext('y1=sin(t)');subplot(3,1,2)plot(t,cos(t),'r:')axis([0,2*pi,-2,2]);xlabel('t');ylabel('cost');title('余弦函数');grid ongtext('y2=sin(t)');subplot(3,1,3)plot(t,tan(t),'r:')axis([0,2*pi,-2,2]);xlabel('t');ylabel('tant');title('正切函数');grid ongtext('y3=sin(t)');作业:t=0:0.1:2*pi;%y1=sin(t); y2=cos(t); y3=tan(t);y4=cot(t) subplot(2,2,1)plot(t,y1,'r-')axis([0,2*pi,-2,2]);xlabel('t');ylabel('y1'); subplot(2,2,2)plot(t,y2,'b--*')axis([0,2*pi,-2,2]);xlabel('t');ylabel('y2'); subplot(2,2,3)plot(t,y3,'k:')axis([0,2*pi,-2,2]);xlabel('t');ylabel('y3'); subplot(2,2,4)plot(t,y4,'m-.',t,y1.*y4,'m-.+') axis([0,2*pi,-2,2]);xlabel('t');ylabel('y4');grid ongtext({'y4=cot(t)'; 'y1*y4'});legend('y4','y1*y4')2、编制程序或其它途径实现下列绘图功能结合上面矩阵的运算及绘图命令读书中P149\152\154程序第2章:了解Simulink的基本操作;熟悉simulink及SPS库中模块,练习搭建:1)P45图2-22a)2)P102图5-3,课后完成P5-29, P111图5-13,3)搭建下下题表达的模型3-5、考虑如图所示感应电机的等效电路,输入的交流电压源为220V、50Hz,其它参数为。
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综合大作业(附详细答案)

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)

一、MATLAB技术结课作业:根据自己专业相关的研究内容,或自己感兴趣的内容,自拟题目,用MATLAB编程语言编写一个程序,实现一定的软件功能。
要求提交三类文件:(1)程序相关文件:.m文件、.mat文件、.fig文件等;(2)一个软件说明书,存成.doc文件,软件说明书的格式和示例见附件1和附件2;(3)编程和调试的记录文件:*.txt(使用diary on 语句记录)。
所有文件打包后,命名为“姓名+学号”,发送至许老师qq邮箱:505231691@二、注意事项:(1)不能使用其他编程语言,如C语言等;(2)不允许相互抄袭;(3)软件说明书是说明自己所编写软件的功能,而不是让说明MATLAB软件;(4)程序不能过于简单,可以通过增加交互输入来增加复杂度。
三、如果找不到合适的题目,可以完成以下题目之一:1. 字符统计,功能:(1)通过人机交互输入英文文献(txt文档)的地址和名称;(2)统计该英文文献中字符a~z出现的频率;(3)以a~z为横坐标,频率为纵坐标,画出频率曲线图,并标出每个频率值。
2. 函数描述,功能:(1)通过人机交互输入一元代数函数f(x)(由多项式系数向量表示);(2)求出该函数的所有关键点(包括与x轴、y轴的交点、局部极值点);(3)绘制该函数的二维图形,并标出上述关键点。
3. 线性方程组求解,功能:(1)通过人机交互输入方程组AX=B 的系数矩阵和增广矩阵;(2)求解方程组并得出解及其意义(包括‘唯一解’, ‘特解’, ‘最小二乘解’ )。
4. 数据拟合和插值,功能:(1)通过人机交互输入n 和数据点坐标向量X,Y ;(2)在图形窗口绘制n (n<20)个子图,给出上述数据的n 次拟合和插值结果(的最小整数)。
5. 完数统计,(若一个数等于它的各个真因子之和,则称该数为完数,如:6=1+2+3,所以6是完数),功能:(1)通过人机交互输入正整数n;(2)列出[1,n]之间的全部完数。
MATLAB结课大作业

一、判断系统稳定性的方法,并举例说明。
方法一:用Nyquist稳定判据判断系统的稳定性Nyquist稳定判据:若想使得闭环系统稳定,则开环系统G(S)H(S)的Nyquist曲线逆时针绕临界点(-1,j0)点的圈数R必需等于G(S)H(S)(系统的开环传函)位于S的右半平面开环极点数P。
即:Z=P-R Z=0 稳定;Z≠0 不稳定,Z为闭环正实部根的个数。
方法二:用Bode图判断系统的稳定性函数调用格式为:margin( )或[Gm Pm wcp wcg]=margin(G)对于最小相位系统:当相角裕度P m(γ)>0o 或幅值裕度G m(h) >1时,表示系统稳定当相角裕度P m(γ)<0o 或幅值裕度G m(h) <1时,表示系统不稳定幅值裕度G m(h)、相角裕度P m(γ)越大,系统稳定程度越好。
在使用时,G m(h)、P m(γ)是成对使用的,有时仅使用一个裕度指标P m(γ)。
方法三:用代数稳定判据法判断系统的稳定性(1)系统数学模型为传递函数形式G(S)=tf(num,den):执行语句:roots(G.den{1});注:“{}”表示维数(2)系统数学模型为零极点增益形式G(S)=zpk(z,p,k);执行语句:G.p{1};(3)系统数学模型为状态空间形式G(S)=ss(A,B,C,D);执行语句:eig(G.A);注:eig()表示计算系统的极点方法四:用根轨迹法判断系统的稳定性若根轨迹在参数取值过程中,部分在左半平面,部分在右半平面,则系统的稳定性与可变参数的取值有关。
函数命令调用格式:[k poles]=rlocfind(G)方法五:用单位阶跃响应曲线判定系统稳定性例:已知系统的开环传函为:5(S+2)G(S)= ----------------------------(S+10)(S³+3S²+2S+5)判断系统的稳定性解:方法一:用Nyquist稳定判据判断系统的稳定性>>G=tf(5*[1 2],conv([1 10],[1 3 2 5]));>>roots(G.den{1})Ans=-10.0000-2.9042-0.0479+1.3112i-0.0479-1.3112i>>nyquist(G)执行以上程序得到nyquist图可以看出,nyquist曲线包围临界点的圈数R=0。
(完整版)武汉大学matlab期末课程作业

(完整版)武汉⼤学matlab期末课程作业“MATLAB及其应⽤”课程作业院系:姓名:学号:联系⽅式:1. 请指出如下5个变量名中,哪些是合法的? abcd-2 xyz_3 3chan a 变量 ABCDefgh答:xyz_3,ABCDefgh 是合法的。
2. 在命令窗中,运⾏命令a=sqrt(2)。
然后请回答以下问题:计算结果a吗?该计算结果只是5近似吗?请在命令窗中,显⽰出具有最多位有MATLAB 数值结果显⽰的默认设置。
该计算结果只是5近似。
3. 命令clear, clf, clc 各有什么⽤处?答:clear 可以清除matlab ⼯作空间中保持的变量。
clf 可以清除图形窗。
clc 清除命令窗中显⽰内容。
4. 想要在MATLAB 中产⽣⼆维数组=987654321S ,下⾯哪些命令能实现⽬的? S=[1,2,3;4,5,6;7,8;9]S=[1 2 3;4 5 6;7 8 9]S=[1,2,3;4,5,6;7,8,9] %整个命令在中⽂状态下输⼊答:第⼆条S=[1 2 3;4 5 6;7 8 9]能实现⽬的。
中⽂状态下逗号不是有效字符。
1.说出以下四条命令产⽣的结果各属于哪种数据类型,是“双精度”对象,还是“符号”对象?3/7+0.1, sym(3/7+0.1), vpa(sym(3/7+0.1),4), vpa(sym(3/7+0.1))答:3/7+0.1结果是双精度。
sym(3/7+0.1)结果是符号。
vpa(sym(3/7+0.1),4)结果是符号。
vpa(sym(3/7+0.1))结果是符号。
过程如图:2.已知a1=sin(sym(pi/4)+exp(sym(0.7)+sym(pi/3)))产⽣精准符号数字,请回答:以下产⽣的各种符号数哪些是精准的?若不精准,误差⼜是多少?能说出产⽣误差的原因吗?a2=sin(sym(pi/4)+exp(sym(0.7))*exp(sym(pi/3)))a3=sin(sym('pi/4')+exp(sym('0.7'))*exp(sym('pi/3')))a4=sin(sym('pi/4')+exp(sym('0.7+pi/3')))a5=sin(sym(pi/4)+exp(sym(0.7+pi/3)))a6=sin(sym(pi/4)+sym(exp(0.7+pi/3)))a7=sin(sym(pi/4+exp(0.7+pi/3)))a8=sym(sin(pi/4+exp(0.7+pi/3)))(提⽰:可⽤vpa观察误差;注意数位的设置)。
Matlab程序设计(2016大作业)

Matlab程序设计课程大作业题目名称:_________________________________ 班级:_________________________________ 姓名:_________________________________ 学号:_________________________________ 课程教师:温海骏学期:2015-2016学年第2学期完成时间:MATLAB 优化应用§1 线性规划模型一、线性规划问题:问题1:生产计划问题假设某厂计划生产甲、乙两种产品,现库存主要材料有A 类3600公斤,B 类2000公斤,C 类3000公斤。
每件甲产品需用材料A 类9公斤,B 类4公斤,C 类3公斤。
每件乙产品,需用材料A 类4公斤,B 类5公斤,C 类10公斤。
甲单位产品的利润70元,乙单位产品的利润120元。
问如何安排生产,才能使该厂所获的利润最大。
问题2:投资问题某公司有一批资金用于4个工程项目的投资,其投资各项目时所得的净收益(投入资金百分比)如下表:工程项目收益表由于某种原因,决定用于项目A 的投资不大于其他各项投资之和而用于项目B 和C 的投资要大于项目D 的投资。
试确定该公司收益最大的投资分配方案。
问题3:运输问题有A 、B 、C 三个食品加工厂,负责供给甲、乙、丙、丁四个市场。
三个厂每天生产食品箱数上限如下表:四个市场每天的需求量如下表:从各厂运到各市场的运输费(元/每箱)由下表给出:求在基本满足供需平衡的约束条件下使总运输费用最小。
§2 多目标规划模型多目标规划定义为在一组约束下,多个不同的目标函数进行优化设计。
数学模型:12min ()()().()0,1,2,,m j f x f x f x st g x j k⎡⎤⎣⎦≤=其中x=(x 1 ,x 2 , … ,x n )为一个n 维向量;f i (x)为目标函数,i=1,2, … ,m; g j (x)为系统约束, j=1, 2, … ,k 。
matlab课后练习习题及答案详解

matlab课后习题及答案详解第1章MATLAB概论与其余计算机语言对比较,MATLAB语言突出的特色是什么?MATLAB拥有功能强盛、使用方便、输入简捷、库函数丰富、开放性强等特色。
MATLAB系统由那些部分构成?MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。
安装MATLAB时,在选择组件窗口中哪些部分一定勾选,没有勾选的部分此后怎样补安装?在安装MATLAB时,安装内容由选择组件窗口中个复选框能否被勾选来决定,能够依据自己的需要选择安装内容,但基本平台(即MATLAB选项)一定安装。
第一次安装没有选择的内容在补安装时只要依据安装的过程进行,不过在选择组件时只勾选要补装的组件或工具箱即可。
MATLAB操作桌面有几个窗口?怎样使某个窗口离开桌面成为独立窗口?又怎样将离开出去的窗口从头搁置到桌面上?在MATLAB操作桌面上有五个窗口,在每个窗口的右上角有两个小按钮,一个是封闭窗口的Close按钮,一个是能够使窗口成为独立窗口的Undock 按钮,点击Undock按钮就能够使该窗口离开桌面成为独立窗口,在独立窗口的view菜单中选择Dock,,菜单项就能够将独立的窗口从头防备的桌面上。
怎样启动M文件编写/调试器?在操作桌面上选择“成立新文件”或“翻开文件”操作时,M文件编写/调试器将被启动。
在命令窗口中键入edit命令时也能够启动M文件编写/调试器。
储存在工作空间中的数组能编写吗?怎样操作?储存在工作空间的数组能够经过数组编写器进行编写:在工作空间阅读器中双击要编写的数组名翻开数组编写器,再选中要改正的数据单元,输入改正内容即可。
命令历史窗口除了能够察看前方键入的命令外,还有什么用途?命令历史窗口除了用于查问从前键入的命令外,还能够直接履行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。
怎样设置目前目录和搜寻路径,在目前目录上的文件和在搜寻路径上的文件有什么差别?目前目录能够在目前目录阅读器窗口左上方的输入栏中设置,搜寻路径能够经过选择操作桌面的file菜单中的SetPath菜单项来达成。
Matlab大作业题目2016

控制系统仿真大作业例题:中央空调控制系统的MA TLAB仿真运用PID算法,利用MATLAB仿真技术对经验法建立的中央空调控制系统模型进行仿真,通过Matlab仿真分析。
1、基于Matlab的变压器差动保护闭环仿真研究(潘高锋)应用Matlab建立微机保护仿真系统,并对不同原理的变压器差动保护进行仿真和比较。
仿真系统采用积木式结构,根据微机保护的实现原理构建模块,实现保护的闭环仿真,对保护的动作过程进行分析。
2、基于MATLAB/SIMULINK的交流电机调速系统建模与仿真(郝士杰)根据直接转矩控制原理,利用MATLAB/SIMULINK软件构造一个交流电机调速系统,,实现高效的调速系统设计。
利用MATLAB/SIMULINK仿真验证该方法的有效性。
3、基于Bode图的电液比例速度控制的设计(张坤玉)对电液比例速度控制系统进行建模分析,并利用Bode图对铣耳机组纵向进给装置电液比例速度控制系统的动态分析和计算机仿真。
设计校正装置使得系统满足① 频域指标:系统的幅角裕量r=72。
;开环穿越频率186rad/s;闭环频率-3dB处的频率是100rad/s;幅值裕量Kg>8dB。
② 时域指标:超调量<3%;上升时间tr=0.05s;调节时间t=0.25s.③ 静态指标:系统静态误差为零。
4、Matlab在动态电路分析中的应用(杨继陈)用Matlab计算动态电路,可得到解析解和波形图。
一阶电路先计算3要素,后合成解析结果;RLC串联和并联的二阶电路采用自编的通用函数计算,自编函数采用了Matlab求解微分方程的符号运算方法;一般的二阶电路和高阶电路采用拉氏变换列写电路方程,再用拉氏反变换得到解析结果.通过实例分析,分析Matlab在动态电路分析方面的优越性.5、大型天线指向控制系统的MA TLAB分析与仿真(杜涵潇)建立大型天线指向控制系统结构图,若要求系统斜披响应的稳态误差小于1%,阶跃响应的超调量小于5%,调节时间小于2s(△=2%),要求:设计合适的校正网络,并绘制校正后系统的单位阶跃响应曲线;当输入为0时,分析扰动对系统输出的影响。
matlab作业练习及答案(有程序和截图)

1. 给出如下方程组:16323313915111081697612534141543s u p w s u p w s u p w s u p w +++=+++=+++=+++=求s ,u ,p ,w 的值,并求出系数行列式。
[答案:s =-0.1258,u =-8.7133,p =11.2875,w =-0.0500。
行列式=7680。
]>> A=[16 32 33 13;5 11 10 8;9 7 6 12;34 14 15 1];>> B=[91 16 5 43]';>> C=A\B2. 求矩阵H1()H X X X X -''=其中1731565419289121110X ⎡⎤⎢⎥⎢⎥=⎢⎥⎢⎥⎣⎦〔答案:H =〔0.7294,0.9041,0.4477,0.9188〕'〕>> X=[17 31 5;6 5 4;19 28 9;12 11 10]; >> H=X/(X'*X)*X'>> X=[17 31 5;6 5 4;19 28 9;12 11 10]; >> H=X*inv((X'*X))*X'3.用MA TLAB语言实现下面的分段函数,()/,||,h x Dy f x h Dx x Dh x D>⎧⎪==≤⎨⎪-<-⎩。
y=h*(x>D)+h/D*x.x>=-D&x<=D)-h*(x<-D)4.636263 0S=2124822 ii==++++++∑①试不采用循环的形式,用数值方法求出上式的解。
②由于数值方法采用double形式进行计算,难以保证有效数字,试采用符号运算的方法求该式的精确解。
并给出保留16位有效数字的结果。
>> sum(2.^[1:63])>> sum(sym(2).^[1:63])5.编写一个矩阵相加函数mat_add(),使其具体的调用格式为,要求该函数能接受任意多个矩阵进行加法运算。
matlab结课作业

《MatLab 在工科数学中的应用》结课作业学校代码: 10128学生姓名: 学号: 学 院: 专 业: 班 级: 指导教师:一、矩阵运算方面的题(参照第一章练习题3)生成矩阵A ,通过修改矩阵A 的第1行第3列元素为2,把矩阵扩充为4行5列,其中第4行第3列元素为7,其余列为0,删除第2列元素的操作生成矩阵B , 按行的逆顺序取A 的1,2,4,5列,且第四行为自然数1到4,生成新的矩阵C ;计算B *C ,B .* C , B +C ,2*B , B 和C 的行列式,B \C , B /C ,并找出A 中大于3的元素,且将其替换为1。
010213468497322A ⎛⎫ ⎪= ⎪ ⎪⎝⎭上机实验步骤 输入>>A=[0,1,0,2,1; 3,4,6,8,4; 9,7,3,2,2];B=A;>>B(1,3)=2;B(4,3)=7;B(:,2)=[] 运行结果如下:输入>>C=A(3:-1:1,[1,2,4,5]) 运行结果如下:输入>>C(4,:)=1:4 运行结果如下: 输入 >>B*C运行结果如下: 输入 >> B.*C运行结果如下:输入>> B+C运行结果如下:输入>> 2*B运行结果如下:输入>> det(C)运行结果如下:输入>> B/C运行结果如下:输入>> B\C运行结果如下:输入>> det(B)运行结果如下:二、绘图方面的题,包括二维和三维图形(参照第二章练习题)(1)、二维图形的绘制用ezplot 绘制函数2(sin ),[0,6]2(1cos )x t t t y t π=-⎧∈⎨=-⎩的图形。
(2)、三维图形的绘制在同一平面中的两个窗口分别用polar 绘制心形线3(1cos )ρθ=-和马鞍面z=x*x-y*y的图形,并以不同的角度观察马鞍面。
MATLAB讲座12作业及答案(附程序)

Matlab讲座1. 抛体运动轨迹的计算: 已知一质量为m的小球(体积不计), 以45度的倾角, 25的初速度斜上抛, 运动过程不考虑空气阻力,求其运动轨迹。
扩展:若考虑阻力,阻力大小与速度正比,此时运动轨迹怎样?Matlab程序:clear all;clc;close all;M=200;X=0;Y=5;theta=pi/4;V=25;Vx=V*cos(theta);Vy=V*sin(theta);k=0.01;dt=0.03;g=9.8;n=1000;axis([0 70 0 25])hold onplot(X,Y,'r*','linewidth',1.5);title('有阻力时物体斜上抛的运动轨迹图')xlabel('横向距离')ylabel('纵向距离')for i=1:nF=k*V^2;Fx=F*Vx/V;Fy=F*Vy/V;Ax=Fx/M;Ay=Fy/M+g;X=(Vx^2-(Vx-Ax*dt)^2)/(2*Ax)+X;Vx=Vx-Ax*dt;Y=(Vy^2-(Vy-Ay*dt)^2)/(2*Ay)+Y;Vy=Vy-Ay*dt;V=sqrt(Vx^2+Vy^2);plot(X,Y,'r*','linewidth',1.5);hold onpause(0.1);if Y<0break;endend2 光的折射现象: 已知,有n种介质, 折射率分别为n^2; 每种介质厚度为1m, 光线求以20度入射的轨迹。
Matlab程序:clear all;clc;close all;figure('color','k')n=5;thin=1;in_theta=20/180*pi;axis([0 5 0 35])hold on% in_theta1=pi/2-asin(sin(in_theta)/(1)^2); % yfirst=tan(in_theta1);yfirst=0;ylast=0;for i=1:nx=i-1:0.01:i;y=tan(in_theta)*x-(yfirst-ylast);% ylast=tan(in_theta)*i;ylast=y(end);plot(x,y,'r','linewidth',2)hold onin_theta=pi/2-asin(sin(in_theta)/(i)^2); yfirst=tan(in_theta)*i;pause(0.01)endfor i=1:nA=[i i];B=[0 35];plot(A,B,'linewidth',2)hold onend3.1 定义两个函数,一个为画圆函数,输入为圆心和半径以及颜色,输出为相应图形;一个画矩形函数,输入为矩形的中心位置和长宽以及颜色,输出为相应图形。
matlab实验报告

MATLAB程序设计——结课作业姓名:学号:班级:已知一组实验数据点x y{EMBED Equation.DSMT4 \* MERGEFORMAT |(,)i iy 00.20.40.60.81 1.2 1.4 1.6 1.7 1.8x02468.411.114.42355130416(1)以函数对实验数据进行拟合,试确定a,b值(a,b值的确定至少采用除两点法以外的一种方法),描绘出实验数据点和拟合曲线在一起的图形,并计算误差。
答:方法一:两点法(two points.m):将实验数据(14.4,1.2)(55,1.6)带入原函数求解a、b:输入程序:>> S=solve('y+14.4*x=14.4/1.2','y+55*x=55/1.6');>> S.xans =0.55110837438423645320197044334975>> S.yans =4.0640394088669950738916256157635可得:a=4.064,b=0.552.则所得函数为:调用two points.m文件对函数进行绘图,并分析误差。
绘图如下:分析结果:设对于数据点,拟合曲线在处的值为,则1)最大误差为= 1.76262)平均误差为=0.79673)均方根误差为=0.98844)误差平方和为=8.7639结论:本次拟合结果误差较小,可以采纳。
方法二:最小二乘法(min2x.m)由已知函数可变换成,将x、y取倒数可得到新数据如下:x’1/x 0.5 0.25 1/6 1/8.4 1/1.1 1/14.4 1/23 1/55 1/130 1/416 y’1/y 1/0.2 1/0.4 1/0.6 1/0.8 1 1/1.2 1/1.4 1/1.6 1/1.7 1/1.8根据新公式,利用最小二乘法可求得a、b的值调用min2x.m文件求解得到a、b的值,并绘制图像、分析误差:a=0.3337,b=8.9944.分析结果:1)最大误差为= 1.01432)平均误差为=0.26723)均方根误差为=0.43144)误差平方和为=1.8610结论:本次拟合结果误差较小,可以采纳。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
计算应用软件
MATLAB程序设计及绘图
班级:
学号:
姓名:
2016.12
题目1:分别利用For循环和While循环求下值:
63
263 0
21222
i
K==++++
∑ 。
For循环代码:
k=0;
for i=0:1:63
k=k+2^i;
end
k
While循环代码:
k=0;
i=0;
while (i<64)
k=k+2^i;
i=i+1;
end
k
结果:k = 1.8447e+019
题目2:已知方程组:
347124
57423
859
652108
x y z w
x y z w
x z w
x y z w
+--=
⎧
⎪-++=-
⎪
⎨
+-=
⎪
⎪-+-+=-
⎩
,求方程组的解。
代码:
A=[3,4,-7,-12;5,-7,4,2;1,0,8,-5;-6,5,-2,10]; B=[4;-3;9;-8];
X=A\B
结果:X =
-1.4841
-0.6816
0.5337
-1.2429
题目3:某试验试件在受火后,其弹性模量发生变化,数据如下:受火温度T=[20, 100, 200, 300, 400, 500, 600, 700, 800, 900];
受火后弹性模E=[21.95,1.85,1.76,1.70,1.55,1.50,1.42,1.35,1.3,1.32]*1011;
用合适的多项式拟合弹性模量随受火温度变化的规律。
代码:
T=[20,100,200,300,400,500,600,700,800,900];
E=[21.95,1.85,1.76,1.70,1.55,1.50,1.42,1.35,1.3,1.28]*10^11;
plot(T,E,'LineWidth',4)
结果:
题目4:(1)编写一个函数式文件,使其产生如下函数段:
22()50.52646
x
x f x x
x x x <=⎧⎪
=-<<=⎨⎪->⎩
(2)编写命令式文件,绘制x 在[0,10]范围内f (x ).f (x +2)的曲线。
(1)代码: function f=Untitled(x) if x<=2 f=2*x; elseif x>6 f=x-4; else
f=5-0.5*x; end
(2)代码: i=0;
for x=0:0.1:10 i=i+1;
y1(1,i)=Untitled(x); end i=0
for x=0:0.1:10 i=i+1;
y2(1,i)=Untitled(x+2); end y3=y1.*y2; x=0:0.1:10;
plot(x,y3,'LineWidth',5)
题目5:设f(x)=x5-4x4+3x2-2x+6;
(1)取x=[-2,8]之间函数的值(取100个点),画出曲线,看它有几个零点。
(提示:用polyval函数)
(2)用roots函数求此多项式的根。
(1)代码:p=[1,-4,3,-2,6];
x=-2:0.1:8;
y=polyval(p,x);
plot(x,y,'LineWidth',5)
(2)代码:
r=roots(p)
n=length(r);
for i=1:n
if r(i)>=-2 && r(i)<=8 if (isreal(r(i)))
r(i)
end
end
end
结果:
ans =
3
ans =
1.6956
题目6:在x =[-10,10],y =[-10,10]范围内画出函数z =的三维图形。
代码:
x=-10:0.01:10; y=-10:0.01:10;
[X,Y]=meshgrid(x,y);
Z=sin(sqrt(X.^2+Y .^2))./sqrt(X.^2+Y .^2); mesh(X,Y ,Z) 结果:
题目7:由100*rand(30)作为全班成绩S (i ),用饼图分别表示60i S <、6070i S ≤<、
7080i S ≤<、8090i S ≤<、90i S ≤所占比例。
代码: S=100*rand(30);
k1=0;k2=0;k3=0;k4=0;k5=0; for i=1:900 if S(i)<60 k1=k1+1;
elseif S(i)>=60&S(i)<70
k2=k2+1;
elseif S(i)>=70&S(i)<80
k3=k3+1;
elseif S(i)>=80&S(i)<90
k4=k4+1;
else
k5=k5+1;
end
end
x=[k1 k2 k3 k4 k5];
y=[1 1 1 1 1];
pie(x,y,{'Si<60','60<=Si<70','70<=Si<80','80<=Si<90','90<=Si'}); 结果:
题目8:列举“MATLAB程序设计课程”在本专业的一个相关应用实例。
求解任意长度为L的简支梁在任意大小的集中F的作用下的弯矩图和剪力图,最大剪力值:F1max和最大弯矩值Mmax:。
代码:
clc
clear;
disp('输入集中力到A点的距离a,梁的长度l,集中力F')
a=input('a (/m)值:');
L=input('L (/m)值:');
F=input('F (/KN)值:');
M=[];
n=0;
for x=0:0.1:L;
if x>=0 && x<=a
n=n+1;
F1(n)=(L-a)/L*F;
M(n)=(L-a)/F*x;
else
n=n+1;
F1(n)=-a/L*F;
M(n)=a/L*F*(L-x);
end
end
x=0:0.1:L;
y2=0*x;
plot(x,M,'linewidth',3)
title('弯矩图')
xlabel('x /m')
ylabel('M /KN*m') figure(2)
plot(x,F1,'linewidth',3) hold on
plot(x,y2,'linewidth',8) legend('剪力F1','梁y=0') title('剪力图')
xlabel('x /m')
ylabel('F1 /KN')
F1max=F*a/L
Mmax=F*a*(L-a)/L
结果:
计算应用软件. 工程力学专业
F1max =
5
Mmax =
250
11 / 11。