matlab实验指导答案详解非常详细正确Word格式文档下载.docx

1、（3）_row 不合法，首字符必须是字母（4）for 不合法，for为MATLAB的保留关键字2、求以下变量的值，并在MATLAB中验证。（1）a = 1 : 2 : 5 ;a = 1 3 5（2）b = a a ;b = 1 1 1 3 3 3 5 5 5（3）c = a + b （ 2 , : ） c =4 6 8实验二MATLAB语言矩阵运算掌握基本的矩阵运算及常用的函数。1、下列运算是否合法，为什么？如合法，结果是多少？（1） result1 = a（2） result2 = a * b 不合法，两矩阵的内阶不相等。（3） result3 = a + b（4） result4 = b

2、* d（5） result5 = b ; c *d（6） result6 = a . * b（7） result7 = a . / b（8） result8 = a . * c不合法，两矩阵的维数不相等，点乘相求两矩阵的维数相同。（9） result9 = a . b（10） result10 = a . 2（11） result11 = a 2不合法，a矩阵不是方阵，方阵才能取平方。（12） result12 = 2 . a以上运算的关键不是结果是多少，而是结果怎么得来的，比如result5 = b ; *d，要学会将其分解后去理解，可在命令窗口输入以下语句：ct=b;r5=t * d2、用

3、MATLAB求下面的的方程组。（1）A=7,2,1,-2;9,15,3,-2;-2,-2,11,5;1,3,2,13b=4;7;-1;0x=Abx = 0.4979 0.1445 0.0629 -0.0813（2） A=1,1,1,0;1,2,1,-1;2,-1,0,-3;3,3,5,-6 b=1;8;3;5xyzw=inv（A）*bxyzw= 1.0000 5.0000 -5.0000 -2.00003、已知 （1）求矩阵A的秩（rank）rank（A）（2）求矩阵A的行列式（determinant）det（A）（3）求矩阵A的逆（inverse）inv（A）（4）求矩阵A的特征值及特征向量

4、（eigenvalue and eigenvector）v,d=eig（A）该题目要求大家学会使用MATLAT的帮助系统，一是知道函数名，怎样查函数用法，二是不知道函数名，怎样查函数名及其用法。4、关系运算与逻辑运算 已知a=20,b=-2,c=0,d=1（1） r1 = a b（2） r2 = a b & c d（3） r3 = a = b* （-10）（4） r4 = b | cr1 =1，r2 = 0，r3 =1，r4 = 0三、思考题,求y=？（用format long g查看y的值）方法一：n=-10:1:10format long gx=2 .ny=sum（x）方法二：y=0;fo

5、r t=-10: y=y+2t;endy方法三t=-10while t=10y=y+2t; t=t+1;实验三选择结构程序设计及调试掌握利用if、switch语句编写选择结构程序，学会MATLAB程序编辑、运行及调试方法。1、求分段函数的值。用if语句实现，算出下列表中x对应的y值。X-52581461955clcclearx=-5 % x=input（请输入x的值）if x=0&xy=x2-5*x+6;else y=x2-x-1;2、分别用if和swith语句实现，将百分制成绩转换为成绩等级A、B、C、D、E。其中90100分为A，8090）分为B，7080）分为C，6070）分为D，60分

6、以下为E。对超出百分制范围的成绩，给出错误提示信息。s=input（请输入百分制成绩:）;if s100 g=ERROR;elseif s=90&s1if rem（n,2）=0 n=n/2; n=3*n+1; a=a,n; a2、根据，当n分别取100、1000、10000时，求x的值分别是多少？n=100for k=1:n y=y+1/k2;x=sqrt（6*y）3、编程求满足的最小m值。以下是三种方法及结果验证方法。sum_x=0;n=0;while sum_x1e4break;re3=necho onsum（2.1:12）13）echo off已知y和t的函数关系：求下面表格中与t对应的

7、y值t0.20.40.60.81.0y=1t=1;f=1;for n=1:20 f=f*n; y=y+tn/f实验五函数的编写及调试掌握MATLAB函数的编写及调试方法。1、编写一个函数，计算下面函数的值，给出x的值，调用该函数后，返回y的值。function y=myfun1（x）选择一些数据测试你编写的函数。=0y=sin（x）;elseif xmax_x max_x=x（k）;if x（k）min_x min_x=x（k）;m_x=sum_x/N;rms_x=sqrt（sum_x2/N）;3、编写一个函数，给出一个向量，生成如下范德蒙矩阵。function v=myvander（x）例如

8、：v=myvander（2 3 4 5）得v=生成一些数据测试你写的函数。 v（k,:）=x.（k-1）;编写程序，用如下迭代公式求，a的值分别为：3，17，113。迭代的终止条件为，迭代初值,迭代次数不超过100次。分别对迭代结果和准确值进行比较，并统计迭代次数。a=input（请输入a=x0=1; x1=x0/2+a/（2*x0）;if（abs（x1-x0）=1e-5）break x0=x1;disp（sqrt（a）数值解为:x1迭代次数:sqrt（a）的准确解为:sqrt（a）请输入a=3x1 = 1.7321ans =实验六 MATLAB的绘图1、在同一坐标系下绘制下面三个函数在t0，

9、4的图象。t=0:1e-2:4*pi;y1=t;y2=sqrt（t）;y3=4*pi*exp（-0.1*t）.*sin（t）;plot（t,y1,r,t,y2,g,t,y3,bgridxlabel（t / sylabel（y1 y2 y3legend（tsqrt（t）4*pi*exp（-0.1*t）.*sin（t）2、编写程序，选择合适的步距，绘制下面函数在区间-6，6中的图象。x=-6:6; % 设定自变量x的取值范围leng=length（x）; % 计算向量x的长度for m=1:leng % 计算函数值 y（m）=myfun1（x（m）;plot（x,y）,grid %绘制函数曲线x）

10、,ylabel（y3、用compass函数画下面相量图ua=1;ub=cos（-2*pi/3）+sin（-2*pi/3）*i;uc=cos（2*pi/3）+sin（2*pi/3）*i;compass（ua,ub,uc,ua-ub,ub-uc,uc-ua）4、三维空间曲线绘制z=0:x=cos（z）;y=sin（z）;plot3（x,y,z）5、用mesh或surf函数，绘制下面方程所表示的三维空间曲面，x和y的取值范围设为-3，3。x,y=meshgrid（-3:0.5:3）z=-x.*x/10+y.*y/10;mesh（x,y,z）;figuresurfc（x,y,z）;附:画抛物面thet

11、a=linspace（0,2*pi,60）R=linspace（0,5,20）;X1=R*cos（theta）;Y1=R*sin（theta）;Z1=X1.*X1/10+Y1.*Y1/10;mesh（X1,Y1,Z1）axis square用对分法求解方程在0，1内的解，并验证，在程序中统计出对分次数。提示：先将原方程转化成的形式。对分法的基本思想是：一个一元方程f（x）=0，若f（x1）*f（x2）0，则在x1,x2区间内有实数解。取该区间的中点xm=（x1+x2）/2，判定f（x1）和f（x2）二者中哪一个与f（xm）异号，若f（x1）*f（xm）1e-10 y1=2*exp（-x1）-s

12、in（x1）; y2=2*exp（-x2）-sin（x2）; xm=（x1+x2）/2; ym=2*exp（-xm）-sin（xm）;if y1*ym x2=xm;elseif y2*ym x1=xm; len=x2-x1;format longxm,ym,n实验七 MATLAB数值运算掌握MATLAB常用的数值运算函数。1、求代数方程的5个根，并将其用星号（*）标记在复平面图上。（用roots和plot函数）。p=3 4 7 2 9 12;r=roots（p）plot（r,Real AxisImag Axisr = -0.8612 + 1.4377i -0.8612 - 1.4377i0.6

13、737 + 1.0159i 0.6737 - 1.0159i -0.9583 2、求代数方程p=1 0 0 0 0 -1;axis equal3、求下面函数在0.5,4区间内的过零点。（用fzero函）fplot（x3-2*x2*sin（x）+5*x*cos（x）+1/x,0.5,4）x1=fzero（,0.5,2）x2=fzero（,2,4）x1 = 1.5117x2 = 2.60954、已知R=50欧姆，U=4V，二极管D正向电流与电压的关系为：其中：Ud为二极管正向电压Is为反向饱合电流，取10-12AK为玻尔茨曼常数，1.38*10-23T为绝对温度，取300开尔文（27摄氏度）q为电

14、子电荷1.6*10-19C求此电路中的电流Id和二极管正向电压Ud（要求用fsolve函数求解）R=50;U=4;Is=1e-12;K=1.38e-23;T=300;q=1.6e-19;%Ud+Id*R-U=0;%Id-Is*exp（Ud*q/（K*T）-1）=0;Ud=0:0.7;Id=Is*exp（Ud*q/（K*T）-1）;Ud1=0:4;Id1=（U-Ud1）/R plot（Ud,1000*Id,hold onplot（Ud1,1000*Id1,hold off二极管伏安特性负载线x=fsolve（UdId,1 0.05）function f=UdId（x）Ud=x（1）;Id=x（2

15、）;f=Ud+Id*R-UId-Is*exp（Ud*q/（K*T）-1）;x =0.6707 0.0666 Ud Id5、实验数据处理：已知某压力传感器的测试数据如下表p0.01.12.12.84.25.06.16.98.19.09.9u111317182224293439p为压力值，u为电压值，试用多项式来拟合其特性函数，求出a,b,c,d，并把拟合曲线和各个测试数据点画在同一幅图上。p=0.0 1.1 2.1 2.8 4.2 5.0 6.1 6.9 8.1 9.0 9.9;u=10 11 13 14 17 18 22 24 29 34 39;plot（p,u,r+abcd=polyfit（

16、p,u,3）p1=0:10;u1=polyval（abcd,p1）;plot（p1,u1,axis（-1 11 8 41）;pup-uu（p）=a*p3+b*p2+c*p+d,2）实验八 MATLAB应用1、用Simulink求解下图所示电路0100微秒内的响应。已知R=6*10-4欧，C=1700微法，L=6*10-9享，uc（0）=15kV。模块参数设置:Integrator1的Initial condition:15kV在命令窗口为R,L,C赋值。仿真参数设置如下：Start time:Stop time:100e-6Solver Type:Variable-stepSolver:ode45 Max step size:1e-7Min step size:autoInitial step size:Relative tolerance:1e-3Absolute tolerance:1e-6实验四 MATLAB数值运算p=3 4 7 2 9 12roots（p） 0.6737 + 1.0159ip2=1 0 0 0 0 -1roots（p2） -0.8090 + 0.5878i