《系统仿真报告》由会员分享,可在线阅读,更多相关《系统仿真报告(20页珍藏版)》请在金锄头文库上搜索。
1、实验一 MATLAB软件环境与应用基础一、实验目的:1、熟悉MATLAB软件基本操作;2、掌握MATLAB运算基础;3、掌握MATLAB基本绘图方法。二、实验内容:1、用逻辑表达式求下列分段函数的值。12,0 t V 1y = t2 -1, 1 t V 2,其中 t=0:0.5:2.5。12 - 2t +1,2 =0)&(t=1)&(t=2)&(t=5),For n=1:length(i)m(n)=A(i(n),j(n)endm程序运行结果如下:m=758693、1 行 100 列的 Fibonacc 数组 a,a(1)=a(2)=1,a(i)=a(i-1)+a(i-2),用 for 循环指令
2、来寻求该数 组中第一个大于10000的元素,并指出其位置i。3、n=100;a=ones(1,n);for i=3:na(i)=a(i-1)+a(i-2);if a(i)10000a(i),break;end;end,i程序运行结果如下:ans=10946i=214、根据y =1 + + + + -,求 3 52n -1(1) y=3my=y-f(n)mn=n-1breakendendmymn程序运行结果如下:my=2.9944mn=56实验二MATLAB绘图与数值计算一、实验目的:1. 掌握数值插值与曲线拟合的方法及其应用。2. 掌握求数值导数、数值积分、代数方程数值求解、常微分方程数值求解
3、的方法。3. 掌握定义符号对象、求符号函数极限及导数、求符号函数积分的方法。二、实验内容:1、某气象观测站测得某日烦6: 0018:00之间每隔2h的室内外温度(C)见表4.1.表4.1室内外温度观测结果(C)时间681012141618室内温度t118.020.022.025.030.028.024.0室外温度t215.019.024.028.034.032.030.0试用3次样条插值分别求出该日室内外6:3017:30之间每隔2h各点的近似温度(C)。(1) h=6:2:18;t1=18 20 22 25 30 28 24;t2=15 19 24 28 34 32 30;h1=6.5:2:
4、17.5;t11=interp1(h, t1, hl, spline)t22=interp1(h, t2, h1, spline)t11 =18.502020.498622.519326.377530.205126.8178t22 =15.655320.335524.908929.638334.256830.95942、已知lgx在1,101区间10个整数采样点的函数值见表4.2.试求lgx的5次拟合多项式 p(X),并绘制出lgx和p(x)在1,101区间的函数曲线。表4.2 lgx在10个采样点的函数值x1112131415161718191101lgx01.04141.32221.491
5、41.61281.70761.78531.85131.90851.95902.0043(2) x=1 11 21 31 41 51 61 71 81 91 101;y=0 1.0414 1.3222 1.4914 1.6128 1.7076 1.7853 1.8513 1.9085 1.9590 2.0043p=polyfit(x,y,5)x1=1:0.5:101;y1=log10(x1);p1=polyval(p, x1)plot(x1, y1, :o, x1, p1, -*)0.0000-0.00000.0001-0.00580.1537-0.1326ln( x + ey)3、求。y项(3
6、) clearfxy=sym(log(x+exp(y)/sqrt(xA2+yA2)result=limit(limit(fxy,x,1),y,0) fxy =log(x+exp(y)/sqrt(xA2+yA2) result = log(2)4、计算I = Jj f (x)dxdy = jj 2(2- x- y)dxdy,其中D为直线y = x2所围部分。(4) clear syms x yf=(2-x-y)/2;y1=x;y2=xA2;X=solve(x-xA2=0)fdy=int(f,y,xA2,x)I=int(fdy,x,X(1),X(2)X =01fdy =x-5/4*xA2-1/2*
7、x*(x-xA2)+1/4*xA4I =11/1205、求下列变上限积分对变量x的导数:jxa+xdx。(5) clearsyms a x t y1 y2y1=sqrt(a+t)y2=int(y1,t,x,xA2)diff(y2,x)y1 =(a+t)A(1/2)y2 =2/3*(a+xA2)A(3/2)-2/3*(a+x)A(3/2)ans =2*(a+xA2)A(1/2)*x-(a+x)A(1/2)6、求解高阶微分方程y 10y + 9y = m。其中矿(0)=, y(0) = 0(6) clearyl=dsolve(D2y-10*Dy+9*y=exp(2*x),Dy(0)=33/7,y(
8、0)=6/7)yi =exp(t)*(3/8-l/8*exp(2*x)+exp(9*t)*(27/56+l/72*exp(2*x)+l/9*exp(2*x)7、设方程刀 + 2y + z 2Jxyz = ,确定了函数z = z(x,y),求和。 ox oy(7) clearsyms x y zf=x+2*y-2*sqrt(x*y*z);fx=diff(f,x);fy=diff(f,y);fz=diff(f,z);zx=-fx/fzzx=-fx/fzzy=-fy/fzzx =-(-l+l/(x*y*z)A(l/2)*y*z)*(x*y*z)A(l/2)/x/yzx =-(-l+l/(x*y*z)
9、A(l/2)*y*z)*(x*y*z)A(l/2)/x/yzy =-(-2+l/(x*y*z)A(l/2)*x*z)*(x*y*z)A(l/2)/x/yxy + 乎=08、解方程组,y-z = l o%2 -5x + 6 = 0(8) clearx,y ,z=solve(x*yA2+zA2=0,y-z= l,xA2-5 *x+6)x =2233y =l/3+l/3*i*2 人(1/2)1/3-1/3*遥2人(1/2)l/4+l/4*i*3 人(1/2)1/4-1/4*遥3人(1/2)-2/3+1/3*i*2(1/2)-2/3-1/3*i*2W/2)-3/4+1/4*i*3A(1/2)-3/4-
10、1/4*i*3(1/2)(9)x=6:1:21;y=4 6.4 8.0 8.4 9.28 9.5 9.7 9.86 10 10.2 10.32 10.42 10.5 10.55 10.58 10.6;x1=6:0.5:21;h=interp1(x,y,x1,spline);plot(x,y,:o,x1,h,-r)A=1,2,4,5,7,6; r=roots(A)r =0.4477 + 1.3908i0.4477 - 1.3908i-0.8749 + 1.2992i-0.8749 - 1.2992i-1.1457实验三Simulink仿真应用一、实验目的:1 .熟悉Simulink环境。2. 熟
11、悉建立Simulink仿真框图并进行系统仿真。3. 熟悉编写S函数和M函数对系统进行仿真。二、实验内容:1、Van der Pol方程所描述系统的仿真。设置初始条件尤(0) = x (0) = 0.25,仿真时间t=20s及Van der Pol方程12J + (2 -1) j + J = 0。要求:(1) 绘出Simulink仿真模块图并进行仿真得到系统的时间响应和相平面图。(2) 编写M函数和M文件得到系统的时间响应和相平面图。(3) 编写S函数并绘出S函数仿真框图,仿真得到系统的时间响应和相平面图。TLAB曰回区IEile Edit Dehug Desktop Window HelpD
12、曲尊苗 | Current Directory: F:WATv .但Shortcuts 回 How to Add 回 Whats NewCurrent D(?) .;1m1 *日回区All Files zlm.mdl晶 zlml.mdlCurrent DirectiComband Hopen( subpl % 13-1 open( subpl subpl subpl % 13-1 open( subpl subpl subpl subplFile Edi t Vi ew Simulation Format Tools Helpojtl tlie Help menu.GainlReady125%
13、is matching: end-tag .?linkccslinkinfo. xmlim st ep size. The simulation step size will be limited to be less ta % 13-11-27 下午2:44 % open( F:MATzlm. mdl) subplot (121 )3 plot (t3 zl313 z2)3 % 13-11-27 下午2:47 % open( F:MATzlm. mdl)% 13-11-27 下午2:53 %3 % 13-11-27 下午3:11 %D HaFile Edit V昼苜AHHLIo o -1 -1 o o