MATLAB实验指导书加程序+上机实例Word格式.docx

1、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:0.1:6*pi;y=5*s in （x）;plot（x,y）个总体了解（6）运行MATLAB的演示程序，demo ,以便对MATLAB 有四、思考题1、 以下变量名是否合法？为什么？（1）x2（2）3col（3）_row（4）for2、 求以下变量的值，并在MATLAB中验证。（1）a =1 : 2 : 5 ;（2）b

2、 =a a（3）c =a + b （ 2 ,:）实验二 MATLAB语言矩阵运算、实验目的：掌握基本的矩阵运算及常用的函数 。、实验内容：1471 2 32 4abcd 854 5 61 3361、下列运算是否合法，如合法，结果是多少？ resultl = a（2）result2 = a * b（3）result3 = a + b（4）result4 = b * d（5）result5 = b ; c * d（6）result6 = a . * b（7）result7 = a . / b（8）result8 = a . * c（9）result9 = a . b（10）resultIO = a

3、 . A2（11）resultll = a 人2（12）resultll = 2 . a a2、用MATLAB求下面的的方程组X1915X2（1） c11X313X4xy z2yz w 82xy3w3x3y5z6w 57 2 129 15 3 23、 已知A2 2 1151 3 2 13求矩阵A的秩（rank）（2）求矩阵A的行列式（determinant）（3）求矩阵A的逆（inverse）求矩阵 A的特征值及特征向量 （eigenvalue and eigenvector）4、 关系运算与逻辑运算已知 a=20,b=-2,c=0,d=1（1） r1 = a b（2） r2 = a b &（

4、3） r3 = a = b* （-10）（4）r4 = b | c10y 2n 2 10n 10三、思考题2 9 2102 2 ，求 y= ?（用 format long 查看 y 的值）实验三程序的编辑及调试掌握MATLAB程序编辑、运行及调试方法。、实验内容:1、启动MATLAB后，点击File|New|M-File ，启动 MATLAB的程序编辑及调试器（Editor/Debugger ），编辑以下程序，点击File|Save保存程序，注意文件名最好用英文字符。点击Debug|Run运行程序，在命令窗口查看运行结果，程序如有错误则改正。运行下面的程序，按程序提示输入n=1,2,3,5,7

5、等数来验证这一结论。%classic 3n+1 problem from n umber theory.while 1n=in put（ En ter n,n egative quits: ）;if nif rem（n, 2）=0 n=n/2;elsen=3*n+1;a=a ,n ;2、编程求满足i 12i 10000的最小m值。用对分法求解方程 2e sinx在o, 1内的解，并验证，在程序中统计出对分次数 。提示：先将原方程转化成f（x） 2e x sinx 0的形式。对分法的基本思想是：一个一元方程f（x）=0，若f（x1）*f（x2）0 ，则在x1,x2区间内有实 数解。取该区间的中点

6、 xm=（x1+x2）/2 ，判定f（x1）和f（x2）二者中哪一个与 f（xm）异号，若f（x1）*f（xm）0 ，则解存在的区间缩小为 x1,xm，否则解存在的区间缩小为 xm,x2。重复这样的步骤，直到区间的长度小于一个可以接受的小数 （比如1e-10 ）则认为中点即是原方程的解。实验四 函数的编写及调试掌握MATLAB函数的编写及调试方法。1、编写一个函数，计算下面函数的值，给出标量x的值，调用该函数后，返回y的值。function y=myfun1（x）sin x,x 6,选择一些数据测试你编写的函数2、编写一个函数求向量 x中元素的平均值、最大值、最小值、均方根值。方均根值（Roo

7、tfunction m_x,max_x,min_x,rms_x=myfun2（x）Mean Square）的计算公式为：用下面数据测试你写的函数（1）x=sin（0:0.01:6*pi） （2）x=rand（1,200），得到的x为200个（0, 1）之间均匀分布的随机数 。得v=生成一些数据测试你写的函数统 计 迭 代 次 数。实验五MATLAB的绘图1、 在同一坐标系下绘制下面三个函数在 t 0 , 4 的图象。y i ty 2 ty 3 4 e 0.1t sin（ t）2、 编写程序，选择合适的步距，绘制下面函数在区间-6，6中的图象。sin x, x 0y（x） x, 0x3x 6,

8、x 33、 用compass函数画下面相量图ua = 1 ; ub = cos（-2*pi/3）+si n（-2*pi/3）*i ; uc=cos（2*pi/3）+si n（ 2*pi/3）*i;compass（ua,ub,uc,ua-ub,ub-uc,uc-ua）4、三维空间曲线绘制 z=0:4*pi;x=cos（z）;y=s nz）；plot3（x,y,z）5、用mesh或surf函数，绘制下面方程所表示的三维空间曲面 ，x和y的取值范围设为-3，3。10 10在同一坐标系下，用不同颜色和线型绘制以下两个函数在 t -2 ，2 范围内的图象。实验六MATLAB数值运算掌握MATLAB常用的

9、数值运算函数。1、 求代数方程3x5 4x4 7x3 2x2 9x 12 0的5个根，并将其用星号（*）标记在 复平面图上。（用roots和plot函数）。5 . _2、 求代数方程x 1 0的5个根，并将其用星号（*）标记在复平面图上。（用roots和plot函数）。3、 求下面函数在0.5,4区间内的过零点。（用fzero函）3 2 1f （x） x 2x sin（x） 5xcos（x）0 1 1_+ 50DU 3Eo 已知R=50欧姆，U=4V ,二极管D正向电流与电压的关系为U dqKTIse其中：Ud为二极管正向电压Is为反向饱合电流，取10-12AK为玻尔茨曼常数，1.38*10-

10、2T为绝对温度，取300开尔文（27摄氏度）q为电子电荷1.6*10 -19C求此电路中的电流Id和二极管正向电压 Ud （要求用fsolve函数求解）5、实验数据处理：已知某压力传感器的测试数据如下表p0.01.12.12.84.25.06.16.98.19.09.9u14171822242934393 2p为压力值，U为电压值，试用多项式u（p） ap bp cp d来拟合其特性函数，求出a,b,c,d，并把拟合曲线和各个测试数据点画在同一幅图上实验七MATLAB应用1、以原点为奇对称中心的方波 y（wt），可以用相应频率的基波及其奇次谐波合成n 1,2,3,取的阶数越多，越接近方波，但总

11、消除不了边缘上的尖峰 ，这称为吉布斯效应。设方波频率为 50Hz,时间 t取00.04 秒（f=50Hz,w=2*pi*f,h=1e-5,tf=40e-3,t = 0:h:tf） ，编写程序，画出如下用1次谐波、1,3次谐波、1,3,5,7,9次谐波，1,3,5, 次1谐波合成的近似方波。（产生方波的函数为：square）2、用Simulink求解下图所示电路0100微秒内的响应。已知R=6*10-4欧,C=1700 微法，L=6*10-9 享，u0）=15kV模块参数设置：Integrator1 的 Initial condition: 15kV在命令窗口为R,L,C赋值。仿真参数设置如下：Start time:Stop time:100e-6Solver Type:Variable-stepSolver:ode45Max step size:1e-7Min step size:autoIn itial step size:Relative tolera nce:1e-3Absolute tolera nce:1e-6MATLAB实验程序 实验1第1题x=2*si n（ 85*pi/