MATLAB应用实验指导书1234结果Word文档下载推荐.docx

1、（6）执行命令x=0:6*pi;y=5*sin（x）;plot（x,y）（7）运行MATLAB的演示程序，demo，以便对MATLAB有一个总体了解。3、矩阵运算 （1）下列运算是否合法，为什么如合法，结果是多少1. result1 = a=1 4;2 5;3 62. result2 = a * b （维数不同）3. result3 = a + b=3 6 2;5 8 114. result4 = b * d=31 22 22;40 49 135. result5 = b ; c * d=31 22 22;40 49 13;-5 -8 76. result6 = a . * b=2 8 -3;

2、4 15 307. result7 = a . / b= ; 8. result8 = a . * c=（维数不同）9. result9 = a . b= ;10. result10 = a . 2=1 4 9;16 25 36 11. result11 = a 2=12. result11 = 2 . a=2 4 8;16 32 64（2）用MATLAB求下面的的方程组。=1, 5 -5 -2X=inv（A）*b（3）已知（1）求矩阵A的秩（rank）=4（2）求矩阵A的行列式（determinant）=12568（3）求矩阵A的逆（inverse）= （4）求矩阵A的特征值及特征向量（ei

3、genvalue and eigenvector）= + - 4、关系运算与逻辑运算 已知a=20,b=-2,c=0,d=1（1） r1 = a b=1（2） r2 = a b & c d=0（3） r3 = a = b* （-10） =1（4） r4 = b | c=0四、思考题1、以下变量名是否合法为什么（1）x2（2）3col（3）_row（4）for定义变量变量名、函数名对字母大小写敏感。MAY、may表示不同变量。变量名第一个字母必须是英文字母，且不能超过65个字符。变量名中不得包含空格、标点但可包含下连符，如my_var 是合法的变量名。2、求以下变量的值，并在MATLAB中验证。

4、（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数值运算与作图掌握MATLAB常用的数值运算函数。1、求代数方程的5个根，并将其用星号（*）标记在复平面图上。（用roots和plot函数）。x =A=3,4,7,2,9,12;x=roots（A）plot（x,*）;grid;2、求代数方程A=1,0,0,0,0,-1;3、求下面函数在,4区间内的过零点。（用fzero函）x = fzero（fun,x0） %查找fun函数在x0附近的零点%估计

5、零点fplot（x3+1/x,4）;hold on;2*x2*sin（x）-5*x*cos（x）hold off;m,n=ginput（2）m =n =y1=fzero（x3-2*x2*sin（x）+5*x*cos（x）+1/x,y1 =y2=fzero（y2 =%建立函数function y=f（x）y=x3-2*x2*sin（x）+5*x*cos（x）+1/x;%调用函数 y1=fzero（fz y2=fzero（4、实验数据处理（选做）已知某压力传感器的测试数据如下表pu1011131417182224293439p为压力值，u为电压值，试用多项式来拟合其特性函数，求出a,b,c,d，并

6、把拟合曲线和各个测试数据点画在同一幅图上。调用方法：polyfit（x,y,n）。用多项式求过已知点的表达式，其中x为源数据点对应的横坐标，可为行向量、矩阵，y为源数据点对应的纵坐标，可为行向量、矩阵，n为你要拟合的阶数用法：linspace（x1,x2,N）功能：linspace是Matlab中的均分计算指令，用于产生x1,x2之间的N点行线性的矢量。其中x1、x2、N分别为起始值、终止值、元素个数。若默认N，默认点数为100。y = polyval（p,x）返回n次多项式p在x处的值。输入变量p是一个长度为n+1的向量，其元素为按降幂排列的多项式系数。y=p1*xn+p2*x（n-1）+.

7、+pn*x+p（n+1） p=,;u=10,11,13,14,17,18,22,24,29,34,39;x=polyfit（p,u,3） %得多项式系数t=linspace（0,10,100）;y=polyval（x,t）; %求多项式得值plot（p,u,t,y,r） %画拟和曲线5、三维空间曲线绘制z=0:4*pi;x=cos（z）;y=sin（z）;plot3（x,y,z） z=0:plot3（x,y,z,rptitle（三维空间曲线text（0,0,0,originxlabel（X）,ylabel（Y）,zlabel（Z6、用mesh或surf函数，绘制下面方程所表示的三维空间曲面，x

8、和y的取值范围设为-3，3。x=-3:3;x,y=meshgrid（x）;z=-x.2/10+y.2/10;mesh（x,y,z）;立体网状图（2）surf（x,y,z）;立体曲面图实验三 MATLAB程序设计掌握MATLAB程序编辑、运行及调试方法。1、熟悉MATLAB的文件启动MATLAB后，点击File|New|M-File，启动MATLAB的程序编辑及调试器（Editor/Debugger），编辑以下程序，点击File|Save保存程序，注意文件名最好用英文字符。点击Debug|Run运行程序，在命令窗口查看运行结果，程序如有错误则改正。注：数论中一个有趣的题目：任意一个正整数，若为偶

9、数，则用2除之，若为奇数，则与3相乘再加上1。重复此过程，最终得到的结果为1。如：2 13 10 5 16 8 4 2 16 3 10 5 16 8 4 2 1运行下面的程序，按程序提示输入n=1,2,3,5,7等数来验证这一结论。%classic 3n+1 problem from number theory.while 1 n=input（Enter n,negative quits: if n1 if rem（n,2）=0 n=n/2; else n=3*n+1; a=a,n; aenda = 1a = 2 1a = 3 10 5 16 8 4 2 1a = 5 16 8 4 2 1a

10、= 7 22 11 34 17 52 26 13 40 20 10 5 16 8 4 2 12.循环语句编程（1）根据，求的近似值。当n=100、1000、10000时，结果是多少 ， ， while 1 k=input（Enter k,negative quits: x=0; for m=1:k x=x+1/m2 x =sqrt（6*x） end （2）用for循环和while循环语句求 1！+2！+10!的值. 4037913k=10;x=0;n=1; n=m*n*1 ; x=x+n; xk=3;m=1; while m=90 disp（ elseif c=80& c90=70&80=60

11、&704. 函数编写（选做）编写一个函数，计算下面函数的值，给出标量x的值，调用该函数后，返回y的值。function y=myfun1（x）选择一些数据测试你编写的函数。if x0&x=3 y=x;3 y=-x+6;运行结果: y=myfun1（-pi/2）y = -1 y=myfun1（0）y = 0 y=myfun1（2） 2 y=myfun1（4）实验四 simulink应用熟悉simulink仿真环境，掌握simulink仿真方法。用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