山东建筑大学数学实验期末作业matlab

matlab⼤作业实验报告,《Matlab程序设计》期末实验报告-⼤作业2015.doc《MATLAB程序设计》实验报告学院： 学号： 姓名：⼀、题⽬：1、(10分)已知矩阵，⽤Matlab代码实现以下要求：(1)将矩阵赋给变量A，并在屏幕上显⽰A；(2)将A按列进列逆序重排，重排后的矩阵赋给变量B，并在屏幕上显⽰B；(3)⽤reshape命令将A重排为⼀个2⾏6列矩阵并赋给变量C；(4)将A重排为⼀个列向量，将其赋给变量D，求D的平均值，在屏幕上显⽰D和它的平均值；(5)⽤命令查看变量A的维数，并显⽰运⾏结果。

2、(10分)写代码实现以下要求：构造菜单项‘Plot’，菜单项Plot有两个⼦菜单项Plot sin(选择此项后执⾏画出曲线，线型为虚线，线条颜⾊为红⾊)和Plot cos(选择此项后执⾏画出曲线 ，线型为实线，线条宽度为2)。

3、(20分)已知，实现下列操作：(1)在同⼀个图形窗⼝，同⼀坐标系下⽤不同的颜⾊和线型绘制三条曲线，并添加图例来区分三条曲线(5分)。

(2)⽤subplot命令，以⼦图的⽅式绘制三条曲线，图形排列⽅式为三⾏⼀列(5分)。

(3) 分别⽤直⽅图(bar)、棒状图(stem)和填充图(fill)绘制三条曲线，以⼦图⽅式绘制，排列⽅式为3⾏3列，共9幅⼦图(10分)。

4、(10分)⽤surf命令绘制曲⾯图形，⽤shading interp命令进⾏插值着⾊处理并添加垂直颜⾊棒。

5、(15分)⾃2011年9⽉1⽇起，我国实⾏新的个⼈所得税征收办法，起征点为3500元，请⽤If-else if-else-end结构实现⼈⼯输⼊⽉收⼊后能计算出个⼈所得税的缴纳额并显⽰⽉收⼊10000元时应缴纳的税款。

级数应纳税所得额x(元)税率备注1x<=15003%x指⽉收⼊扣除起征点3500元之后的余额；215008000045%同上6. (10分)⽤while-end循环结构计算级数和的值，输⼊n值，能计算出f的值，并显⽰结果。

实验报告（MATLAB课后作业练习题）学院电子信息学院班级学号姓名任课教师目录实验作业1 (3)第一题、一阶电路 (3)实验作业2 (7)第一题Waterfall Scope（瀑布显示图） (7)Chirp Signal扫频信号源 (7)Uniform Random Number信号源下 (8)Band-Limited White Noise信号源 (8)第二题：设计一个编程开关仿真系统框图 (9)仿真实验作业3 (10)第一题 (10)第二题 (13)仿真实验作业4 (14)第一题 (14)第二题 (16)仿真实验作业5 (19)仿真实验作业6 (21)仿真实验作业7 (23)仿真实验作业8 (26)实验作业1第一题、一阶电路(1)、电路图如下，R=1.4欧，L=2亨，C=0.32法，初始状态：电感电流为零，电容电压为0.5V ，t=0时刻接入1V 的电压，求0<t<15s 时，i (t)，v o (t)的值，并且用Simulink 仿真画出R=1.4、R=5和 R=9的电流与电容电压的关系曲线。

还可以进一步修改信号源参数，使用三角波、正弦波等作为激励信号，观察输出信号的情况。

function xdot=funcforexl23(t,x,~,R,L,C)xdot=zeros(2,1); %矩阵初始化 xdot(1)=-R/L*x(1)-1/L*x(2)+1/L* f(t); %方程1 xdot(2)=1/C * x(1); %方程2 function in=f(t) %输入信号 in=(t>0)*1; %阶跃信号%filename ex123.mL=2; %电感值 C=0.32; %电容值for R=[1.4 5 9] %仿真电阻值分别为1.5, 3, 5欧姆的情况[t,x]=ode45('funcforexl23',[0,15],[0;0.5],[],R,L,C); %也可采用ode23, ode15s 等求解figure(1);plot(t,x(:,1));hold on ; xlabel('time sec'); text(2,0.07,'\leftarrow i_L(t)');grid;figure(2);plot(t,x(:,2));hold on ;xlabel('time sec'); text(2.1,0.75,'\leftarrow u_C(t)');grid; End输入输出的传递函数：11)()()(2++==RCs LCs s F s U s H c① R=1.4时：1448.064.01)(2++=s s s H ±Vs=1Vt=0R L C +-)(t i )(t v o② R=5时：16.164.01)(2++=s s s H③ R=9时：188.264.01)(2++=s s s H连续系统的传递函数如下：借助多项式乘法函数conv 来处理：两个向量分别用num 和den 表示。

《MATLAB语言及应用》期末大作业题目1．数组的创建和访问（20分，每小题2分）：1)利用randn函数生成均值为1，方差为4的5*5矩阵A；2)将矩阵A按列拉长得到矩阵B；3)提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的矩阵C；4)寻找矩阵A中大于0的元素；]5)求矩阵A的转置矩阵D；6)对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E；7)删除矩阵A的第2列和第4列得到矩阵F；8)求矩阵A的特政值和特征向量；9)求矩阵A的每一列的和值；10)求矩阵A的每一列的平均值；程序代码：clear;clc;A=1+sqrt(4)*randn(5) %生成均值为1，方差为4的5*5矩阵A；B=A(:) %将矩阵A按列拉长得到矩阵B;C=A([2 3],[2 4]) %提取矩阵A的第2行、第3行、第2列和第4列元素组成2*2的矩阵C；n=find(A>0) %寻找矩阵A中大于0的元素；x=A(n)D=A' %求矩阵A的转置矩阵D;E1=flipud(A); %对矩阵A进行上下对称交换后进行左右对称交换得到矩阵E；E=fliplr(E1)F=A(:,[1 3 5]) %删除矩阵A的第2列和第4列得到矩阵F；[Av,Ad]=eig(A) %求矩阵A的特征值和特征向量；S=sum(A,1) %求矩阵A的每一列的和值；Avg=S/5 %求矩阵A的每一列的平均值；运行结果：A =2.3333 2.1171 0.8568 2.1971 -0.7526-1.7853 0.4453 -3.8292 1.2944 0.4690 -1.6011 -1.5874 -0.3887 0.7971 0.3448 -0.2100 -0.7769 -1.7828 -4.2700 -1.3165 -1.9771 -0.9730 1.6593 1.0561 2.1601B =2.3333-1.7853-1.6011-0.2100-1.97712.11710.4453-1.5874-0.7769-0.97300.8568-3.8292-0.3887-1.78281.65932.19711.29440.7971-4.27001.0561-0.75260.46900.3448-1.31652.1601C =0.4453 1.2944-1.5874 0.7971n =167111516171820222325x =2.33332.11710.44530.85681.65932.19711.29440.79711.05610.46900.34482.1601D =2.3333 -1.7853 -1.6011 -0.2100 -1.97712.1171 0.4453 -1.5874 -0.7769 -0.97300.8568 -3.8292 -0.3887 -1.7828 1.65932.1971 1.2944 0.7971 -4.2700 1.0561-0.7526 0.4690 0.3448 -1.3165 2.1601E =2.1601 1.0561 1.6593 -0.9730 -1.9771-1.3165 -4.2700 -1.7828 -0.7769 -0.21000.3448 0.7971 -0.3887 -1.5874 -1.60110.4690 1.2944 -3.8292 0.4453 -1.7853-0.7526 2.1971 0.8568 2.1171 2.3333F =2.3333 0.8568 -0.7526-1.7853 -3.8292 0.4690-1.6011 -0.3887 0.3448-0.2100 -1.7828 -1.3165-1.9771 1.6593 2.1601Av =Columns 1 through 40.1004 + 0.2832i 0.1004 - 0.2832i0.6302 -0.5216-0.5969 -0.5969 -0.4811 0.0856 -0.4405 + 0.0006i -0.4405 - 0.0006i -0.3078 0.21200.2732 - 0.4899i 0.2732 + 0.4899i0.0244 -0.1780-0.0617 + 0.2024i -0.0617 - 0.2024i 0.5254 0.8025 Column 50.3903-0.49590.02180.1929-0.7511Ad =Columns 1 through 4-2.6239 + 1.7544i 0 0 00 -2.6239 - 1.7544i 0 00 0 -0.2434 00 0 0 3.54550 0 0 0Column 50 0 0 02.2257S =-3.2403 -0.7749 -3.4846 1.0747 0.9049Avg =-0.6481 -0.1550 -0.6969 0.2149 0.18102．符号计算（10分，每小题5分）：1) 求方程组20,0uy vz w y z w ++=++=关于,y z 的解；2) 利用dsolve 求解偏微分方程,dx dyy x dt dt==-的解；程序代码：clc[u,v,w] = solve('u*y^2 + v*z + w = 0','y + z + w = 0','u,v,w')[x y]=dsolve('Dx=y','Dy=-x')运行结果：u =(-v*z+y+z)/y^2v =vw = -y-zx =-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:4a =；程序代码：clccleara=0.5:0.5:4.5;t=-2*pi:0.1:2*pi;N=length(a); for i=1:1:Nx=a(i)*cos(t);y=sqrt(25-a(i).^2).*sin(t);plot(x,y)hold onend运行结果：-5-4-3-2-1123452) 利用plotyy 指令在同一张图上绘制sin y x =和10x y =在[0,4]x ∈上的曲线；程序代码：clcclearx=0:0.01:4;y1=sin(x); y2=10.^x;plotyy(x,y1,x,y2) %用双y 轴绘制二维图形运行结果：3) 用曲面图表示函数22z x y =+；程序代码：clcclear[X,Y] = meshgrid(-2:0.05:2); %产生xy 平面上的网格数据Z = X.^2 + Y.^2;surf(X,Y,Z) %绘着色曲面图hold off运行结果：4) 用stem 函数绘制对函数cos 4y t π=的采样序列；程序代码：clcclearfs=25;Ts=1/fs;n=1:1:200;yn=cos(pi*n*Ts/4);stem(n,yn) %绘离散数据的火柴杆图运行结果：4. 设采样频率为Fs = 1000 Hz ，已知原始信号为)150π2sin(2)80π2sin(t t x ⨯+⨯=，由于某一原因，原始信号被白噪声污染，实际获得的信号为))((ˆt size randn x x +=，要求设计出一个FIR 滤波器恢复出原始信号。

期末综合实验报告1、（1）设矩阵A 、B 满足关系式AB=2B+A ，且， 求B.(2)生成一个3×7的随机数矩阵，并将其各列排序、求各列的最大值、各列的最小值、各列的平均值和各列的元素之和。

解：（1）Matlab 源代码 A=[3 2 5;0 1 4;2 3 1];B=(A-2*eye(3))\A,B=inv(A-2*eye(3))*A 解得 B =-0.4667 2.2667 1.7333 1.0667 -0.4667 -0.5333 0.2667 0.1333 0.8667 B =,132410523⎥⎥⎥⎦⎤⎢⎢⎢⎣⎡=A-0.4667 2.2667 1.73331.0667 -0.4667 -0.53330.2667 0.1333 0.8667(2) Matlab源代码A=rand(3,7),Asort=sort(A),Amax=max(A),Amin=min(A),Amean=mean(A,1),Asum=sum(A)解得A =0.8147 0.9134 0.2785 0.9649 0.9572 0.1419 0.79220.9058 0.6324 0.5469 0.1576 0.4854 0.4218 0.95950.1270 0.0975 0.9575 0.9706 0.8003 0.9157 0.6557Asort =0.1270 0.0975 0.2785 0.1576 0.4854 0.1419 0.65570.8147 0.6324 0.5469 0.9649 0.8003 0.4218 0.79220.9058 0.9134 0.9575 0.9706 0.9572 0.9157 0.9595Amax =0.9058 0.9134 0.9575 0.9706 0.9572 0.9157 0.9595Amin =0.1270 0.0975 0.2785 0.1576 0.4854 0.1419 0.6557Amean =0.6158 0.5478 0.5943 0.6977 0.7476 0.4931 0.8025Asum =1.8475 1.6433 1.78292.0931 2.2428 1.47942.40742、设（1）在第一个子图中画出该函数的三维图形（surf ） ; （2）在第二个子图中绘制该函数的等值线和梯度场。

MATLAB数学实验答案（全）第⼀次练习教学要求：熟练掌握Matlab 软件的基本命令和操作，会作⼆维、三维⼏何图形，能够⽤Matlab 软件解决微积分、线性代数与解析⼏何中的计算问题。

补充命令vpa(x,n) 显⽰x 的n 位有效数字，教材102页fplot(‘f(x)’,[a,b]) 函数作图命令，画出f(x)在区间[a,b]上的图形在下⾯的题⽬中m 为你的学号的后3位（1-9班）或4位（10班以上） 1.1 计算30sin limx mx mx x →-与3sin lim x mx mxx →∞-syms xlimit((902*x-sin(902*x))/x^3) ans =366935404/3limit((902*x-sin(902*x))/x^3,inf)//inf 的意思 ans = 0 1.2 cos1000xmxy e =，求''y syms xdiff(exp(x)*cos(902*x/1000),2)//diff 及其后的2的意思 ans =(46599*cos((451*x)/500)*exp(x))/250000 - (451*sin((451*x)/500)*exp(x))/250 1.3 计算221100x y edxdy +??dblquad(@(x,y) exp(x.^2+y.^2),0,1,0,1)//双重积分 ans = 2.13941.4 计算4224x dx m x +? syms xint(x^4/(902^2+4*x^2))//不定积分 ans =(91733851*atan(x/451))/4 - (203401*x)/4 + x^3/12 1.5 (10)cos ,x y e mx y =求//⾼阶导数syms xdiff(exp(x)*cos(902*x),10) ans =-356485076957717053044344387763*cos(902*x)*exp(x)-3952323024277642494822005884*sin(902*x)*exp(x)1.6 0x =的泰勒展式(最⾼次幂为4).syms xtaylor(sqrt(902/1000+x),5,x)//泰勒展式 ans =-(9765625*451^(1/2)*500^(1/2)*x^4)/82743933602 +(15625*451^(1/2)*500^(1/2)*x^3)/91733851-(125*451^(1/2)*500^(1/2)*x^2)/406802 + (451^(1/2)*500^(1/2)*x)/902 +(451^(1/2)*500^(1/2))/500 1.7 Fibonacci 数列{}n x 的定义是121,1x x ==12,(3,4,)n n n x x x n --=+=⽤循环语句编程给出该数列的前20项（要求将结果⽤向量的形式给出）。

MATLAB数学实验第二版课后练习题含答案课后练习题MATLAB数学实验第二版的课后练习题如下：第一章课后练习题1.编写MATLAB程序，计算并输出下列公式的结果：y = \\frac{1}{\\sqrt{2\\pi\\sigma^2}} e^{-\\frac{(x-\\mu)^2}{2\\sigma^2}}其中，x, $\\mu$, $\\sigma$ 分别由用户输入。

要求输出结果精确至小数点后两位。

答案如下：x=input('请输入 x 的值:');mu=input('请输入 mu 的值:');sigma=input('请输入 sigma 的值:');y=1/sqrt(2*pi*sigma^2) *exp(-(x-mu)^2/ (2*sigma^2));fprintf('y = %.2f\', y);2.编写MATLAB程序，求解下列方程的解：4x + y = 11\\\\x + 2y = 7答案如下：A= [4,1;1,2];B= [11;7];X=inv(A) *B;fprintf('x = %.2f, y = %.2f\', X(1), X(2));第二章课后练习题1.编写MATLAB程序，计算下列多项式的值：P(x) = x^4 - 2x^3 + 3x^2 - x + 1其中，x 由用户输入。

要求输出结果精确至小数点后两位。

答案如下：x=input('请输入 x 的值:');y=x^4-2*x^3+3*x^2-x+1;fprintf('P(%.2f) = %.2f\', x, y);2.编写MATLAB程序，绘制下列函数的图像：f(x) = \\begin{cases} x + 1, & x < 0 \\\\ x^2, & 0 \\leq x < 1 \\\\ 2x - 1, & x \\geq 1 \\end{cases}答案如下：x=-2:0.01:2;y1=x+1;y2=x.^2.* ((x>=0) & (x<1));y3=2*x-1;plot(x,y1,x,y2,x,y3);legend('y1 = x + 1','y2 = x^2','y3 = 2x - 1');总结本文提供了《MATLAB数学实验第二版》的部分课后练习题及其答案。

重庆大学学生实验报告实验课程名称数学实验开课实验室学院年级专业班学生姓名学号开课时间至学年第学期数理学院制开课学院、实验室：实验时间：年月日2）比较函数33)(,)(,)(xxhxxxgxxf=+==在x→0时函数的性态。

观察到什么现象？从观察到的现象，反映了什么结论。

程序：x=-1:0.0001:1;y1=x;y2=x.^3;y3=y1+y2;plot(x,y1,x,y2,x,y3)结果：结论：当x→0时，f(x)与g(x)很接近，而h(x)与前两个函数都不接近。

3）比较函数33)(,)(,)(xxhxxxgxxf=+==在x→∞时函数的性态。

程序如下所示：x=linspace(-100000,100000,30); y1=x;y2=x+x.^3;y3=x.^3;subplot(2,2,1),plot(x,y1),title('f(x)=x'),xlabel('x');ylabel('f(x)'); grid; subplot(2,2,2),plot(x,y2),title('g(x)=x+x^3'),xlabel('x');ylabel('g(x)'); grid;subplot(2,2,3),plot(x,y3),title('h(x)=x^3'),xlabel('x');ylabel('h(x)');grid ；结果：4）在日常生活中我们有这样的经验：与幂函数相比，指数函数是急脾气，对数函数是慢性子。

这就是说，当x→∞时，再小的指数函数也比幂函数变化快，再大的对数函数也比幂函数变化慢。

当x→∞时,比较10xy=与xy1.1=的大小.当x→∞时,比较001.0xy=与xy lg1000=的大小.程序如下所示：x=linspace(5000,8000,500);y1=x.^10;y2=1.1.^x;Subplot(1,2,1),plot(x,y1),xlabel('x');ylabel('y)'); grid; title('y=x^1^0');Subplot(1,2,2),plot(x,y2),xlabel('x');ylabel('y)'); grid; title('y=1.1^x');结果：从上图可以看出来指数函数变化快程序如下所示：x=linspace(5000,8000,500);y1=x.^0.001;y2=1000.*log(x);Subplot(1,2,1),plot(x,y1),xlabel('x');ylabel('y)'); grid; title('y=x^0.001');Subplot(1,2,2),plot(x,y2),xlabel('x');ylabel('y)'); grid; title('y=1000.*log(x)');结果：分析：由以上函数图形可知对数函数变化比幂函数慢。

数学实验练习2.1画出下列常见曲线的图形。

（其中a=1,b=2,c=3）1、立方抛物线3xy=解：x=-5:0.1:0;y=(-x).^(1/3);y=-y;x=0:0.1:5;y=[y,x.^(1/3)];x=[-5:0.1:0,0:0.1:5];plot(x,y)2、高斯曲线2x e=y-解：fplot('exp(-x.^2)',[-5,5])3、笛卡儿曲线)3(13,1333222axy y x t at y t at x =++=+=解：ezplot('x.^3+y.^3-3*x*y',[-5,5])xyx.3+y.3-3 x y = 0或t=-5:0.1:5; x=3*t./(1+t.^2); y=3*t.^2./(1+t.^2); plot(x,y)4、蔓叶线)(1,1322322xa x y t at y t at x -=+=+=解：ezplot('y.^2-x.^3/(1-x)',[-5,5])xyy.2-x.3/(1-x) = 0或t=-5:0.1:5; x=t.^2./(1+t.^2); y=t.^3./(1+t.^2); plot(x,y)5、摆线)cos 1(),sin (t b y t t a x -=-= 解：t=0:0.1:2*pi;x=t-sin(t); y=2*(1-cos(t)); plot(x,y)6、星形线)(sin ,cos 32323233a y x t a y t a x =+== 解：t=0:0.1:2*pi; x=cos(t).^3; y=sin(t).^3;plot(x,y)或ezplot('x.^(2/3)+y.^(2/3)-1',[-1,1])xyx.2/3+y.2/3-1 = 07、螺旋线ct z t b y t a x ===,sin ,cos 解：t=0:0.1:2*pi; x=cos(t); y=2*sin(t); z=3*t; plot3(x,y,z) grid on8、阿基米德螺线θa r = 解：x =0:0.1:2*pi; r=x; polar(x,r)902701809、对数螺线θa e r = 解：x =0:0.1:2*pi; r=exp(x); polar(x,r)90270180010、双纽线))()((2cos 22222222y x a y x a r -=+=θ 解：x=0:0.1:2*pi; r=sqrt(cos(2*x)); polar(x,r)90270或ezplot('(x.^2+y.^2).^2-(x.^2-y.^2)',[-1,1]) grid onxy(x.2+y.2).2-(x.2-y.2) = 011、双纽线)2)((2sin 222222xy a y x a r =+=θ 解：x=0:0.1:2*pi; r=sqrt(sin(2*x)); polar(x,r)90270或ezplot('(x.^2+y.^2).^2-2*x*y',[-1,1]) grid onxy(x.2+y.2).2-2 x y = 012、心形线)cos 1(θ+=a r 解：x =0:0.1:2*pi; r=1+cos(x); polar(x,r)90270练习2.21、求出下列极限值。

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与数学实验的考试试题一、单项选择题（每题2分，共10分）1. MATLAB的全称是什么？A. Matrix LaboratoryB. Microprocessor Application ToolC. Microsoft Advanced Technology ToolD. Microprocessor Application Technology2. 在MATLAB中，以下哪个命令用于绘制函数f(x)=x^2在闭区间[0,1]上的图像？A. plot(0:1, 0:1)B. plot(0:0.01:1, 0:0.01:1)C. plot(0:1, 0:1:1)D. plot(0:0.01:1, 0.^2)3. 以下哪个MATLAB命令用于求解线性方程组？A. solveB. linsolveC. equationD. linear4. 在MATLAB中，用于生成一个3x3单位矩阵的命令是什么？A. eye(3)B. unit(3)C. identity(3)D. I(3)5. 如果变量x和y在MATLAB中分别表示为x = [1 2 3; 4 5 6] 和 y= [1; 2; 3]，那么表达式x * y的结果是什么？A. [5; 15; 29]B. [14; 32; 50]C. [7; 15; 23]D. [3; 6; 9]二、简答题（每题5分，共20分）1. 简述MATLAB在数学实验中的作用和重要性。

2. 解释MATLAB中向量和矩阵的区别，并给出创建它们的基本命令。

3. 在MATLAB中，如何使用for循环生成一个从1到100的奇数向量？4. 描述在MATLAB中使用函数文件的过程，包括如何定义和调用函数。

三、编程题（每题10分，共30分）1. 编写一个MATLAB函数，该函数接受一个向量作为输入，并返回向量中所有元素的和。

```matlabfunction S = sumVector(V)% 请在此处编写代码end```2. 编写一个MATLAB脚本，该脚本生成一个5x5的随机矩阵，并计算其行列式。

《Matlab 程序设计》 期末作业学院： 专业： 班级：学号： 姓名： 成绩：一 、 简答题（每题10分，共50分）1、假设有两个矩阵A 和B ，请用这两个矩阵说明Matlab 数值计算中的“左除”与“右除”有什么区别？2、脚本M 文件与函数M 文件的主要区别是什么？3、设矩阵A=[1 4 3 9；5 30 12 42；2 18 32 15]，试编程找出A 中所有大于或等于5的元素，并把这些元素组成一个列向量，并计算出这些元素的总和（注意附上标注说明）4、利用Matlab 编写程序求取335lim 21x x x →∞++和120()t t dt -⎰，并对求取结果进行分子和分母的提取（注意附上标注说明）。

5、利用Matlab 编写程序求微分方程x e y y y x 2cos 52=+'-''的通解。

二、 综合题（每题20分，共40分）1、在同一张图上分别画出两组坐标)10,,2,1(),4,(),,(3222 =+i i i i i i 的散点图和折线图。

要求只能利用plot 命令进行绘制，X 坐标范围[0,105] ，Y 坐标范围[0,1450]。

2、已知学生的名字和百分制分数。

编写相应的程序要求输入所有学生的百分制分数后，一次性判断每个学生的成绩分别是“满分”、“优秀”、“良好”、“及格”还是“不及格”，并要求最后用单元数组按行分别输出每个学生的名字、分数和成绩判断结果。

学生姓名 Jack Marry Peter Rose Tom成 绩 72 83 56 94 100判断标准：满分(100)，优秀(90-99)，良好(80-89)，及格(60-79)，不及格(<60)三、 论述题（共10分）浅谈一下Fortran 、matlab 、maple 和mathematic 这四种计算语言的特点和区别。

MATLAB程序设计期末大作业姓名：班级：学号：指导教师：题目1给定如图1所示的单位负反馈系统。

图1在系统中分别引入不同的非线性环节（饱和、死区、与滞环），观察系统的阶跃响应，并且分析比较不同的非线性环节对系统性能的影响。

解：1、利用MATLAB中的simulink工具箱，对题设控制系统进行建模，如图1 所示。

则没有任何非线性环节的系统，其阶跃响应曲线如图2 所示。

图22、在系统中加入饱和非线性环节，系统框图如图3所示。

其中，饱和非线性环节的输出上限为0.1，输出下限为-0.1；阶跃信号幅值为1。

图3利用simulink进行仿真，得到的阶跃响应曲线如图4所示。

图4为了比较当饱和非线性环节输出上下限变化时系统阶跃响应的不同，可以利用simulink中的to workspace模块，将多次仿真的结果记录到工作空间的不同数组中，并且绘制到同一幅图像上。

此时，系统框图如图5所示。

图5将4种情况下系统的阶跃响应曲线绘制在同一幅图像中，代码如下：>> plot(tout,out2);>> plot(tout,out2);>> hold on;>> grid on;>> gtext('0.1');>> plot(tout,out1);>> plot(tout,out3);>> gtext('0.2');>> plot(tout,out4);>> gtext('0.5');运行程序，结果如图6所示。

图6从图6中可以看出：当饱和非线性环节的输出范围较窄时，系统的阶跃响应速度较慢，上升时间长；同时，超调量较小，振荡不明显;随着输出范围的扩大，系统的响应速度加快，上升时间大大减小，同时伴有显著的振荡。

这是因为饱和环节会对信号起到限幅作用；不难想象，限制作用越强，系统的输出越不容易超调，响应也会越慢，这从图6也可以看出这一趋势。

实验4 MATLAB 程序设计一、实验目的：1．理解命令式M 文件和函数式M 文件的关系及各自的特点；2．掌握MATLAB 中关系运算和逻辑运算的正确使用方法；3．掌握MATLAB 中两种循环结构及选择结构的正确使用方法；4．能够编制简单的程序解决实际问题。

二、实验内容：1．满足条件01=a ，12=a ，…，n n n a a a +=++12（n=0，1，…)的数列称为Fibonacci 数列，（1）编写一个MATLAB 函数式M 文件fibo.m 来计算Fibonacci 数列的任意项的数值；（2）计算数列中第50项的数值并给出程序运行所用的时间；（3）编写一个MATLAB 命令式M 文件调用函数式M 文件fibo.m 来计算Fibonacci 数列的前n 项和的值；给出n=50时的结果。

（1）计算Fibonacci 数列的程序function x=fibo(n)if n<1|fix(n)~=ndisp('N must be a positive integer.');elseif n==1;x=0;elseif n==2x=[0,1];elsex=zeros(1,n);x(1)=0;x(2)=1;for k=3:nx(k)=x(k-1)+x(k-2);endend（2）计算数列中第50项的数值x= fibo(50);ticF=x(50)toc（3）Fibonacci 数列的前50项和的值x= fibo(50);FF=sum(x)2．利用下面的公式计算π的值14134171513114---++-+-=N N π，选择不同的N 值，比较结果的精确度。

function PI=hw2(n)if n<1|fix(n)~=ndisp('N must be a positive integer.');elsePI=0;for k=1:nx=1/(4*k-3)-1/(4*k-1);PI=PI+x;endendPI=4*PI;n N的最小正整数n。

精选文档MATLAB期末考试一试卷及其参照答案一、填空题1、MATLAB常用操作界面包含、工作空间窗口、、、内存数组编写器、M 文件编写 / 调试器、帮助导航/阅读器、图形窗口等。

2、MATLAB中Inf或inf表示、NaN或nan表示、nargout表示。

3、MATLAB中逗号主要用作；用作输入量与输入量之间的分开符；。

4、工作空间阅读器主要用于内存变量的、和。

5、MATLAB实现将全下标变换为单下标的指令为、据单下标换算出全下标的指令为。

6、二维数组的表记有、、“逻辑1”表记。

二、简答题1、简述 MATLAB历史指令窗的主要作用。

2、简述空数组的功用。

3、简述 MATLAB函数的基本结构。

4、简述绘制二维图形的一般步骤。

三、阅读程序并回答以下问题（每题 4 分，共28 分）1、写出以下指令运转结果。

A=zeros(2,4);A(:)=1:8;s=[2 3 5];A(s)Sa=[10 20 30]'A(s)=Sa2、写出以下指令运转结果。

A=reshape(1:16,2,8)reshape(A,4,4)s=[1 3 6 8 9 11 14 16];A(s)=03、写出以下指令运转结果。

A=[1,2;3,4];B=[-1,-2;2;1];S=3;A.*BA*BS.*AS*B4、下边的函数主要达成什么功能？function f=factor(n)if n<=1f=1;elsef=factor(n-1)*n;end5、写出以下指令运转结果。

ch=‘ ABc123d4e56Fg9’;subch=ch(1:5)k=find(ch>= ‘ a’ &ch<=‘ z’ );ch(k)=ch(k)- ( ‘-a‘’A’ );char(ch)6、写出以下指令运转结果。

A(1,1)={'this is cell'};A{1,2}={[1 2 3;4 5 6]};A{2,1}=[1+2*i];A{2,2}=A{1,2}{1}+(A{1,2}{1}(1,1)+A{1,2}{1}(2,2));celldisp(A)7、下边的程序达成功能是什么？t=0:pi/50:4*pi;y0=exp(-t/3);y=exp(-t/3).*sin(3*t);plot(t,y,'-r',t,y0,':b',t,-y0,':b')xlabel( ‘it t’ );ylabel( ‘it y’ );grid on;四、编程题（ 32 分）1、 1、在同一图上分别用红色实线和绿色虚线绘制y1=sin(x) 和 y2=cos(x)在区间 [0， 4*pi] 的曲线，并用星号* 标出两条曲线的交点以及成立图例。

matlab实验报告和Simulink方法及仿真实验山东建筑大学姓名:XX学号201009班级: 095班山东建筑大学机电工程学院2013.07.01~2013.07.07第二部分Simulink方法及仿真实验一、开环系统1、如图所示，当输入信号分别是方波信号和正弦信号时，分析输出信号的波形。

积分环节的特点是输出量为输入量对时间的累积，输出幅值呈线性增长。

（a）输入信号是方波信号输入信号波形图输出信号波形图分析：当输入信号时方波信号时，输出信号等于对方波信号时间的累积。

当方波信号在高电位时，输出为线性增长；当方波信号在零电位时，输出信号不增加，并且是保持不变的，因此是一条直线，在时间的积累下，输出信号就是如图所示的波形。

2、如图所示，开环系统是一阶惯性环节的开环系统。

当输入信号分别是方波信号和正弦信号时，分析输出信号的波形。

一阶惯性环节逇输出需要延迟一段时间才能够接近所要求的输出量，但是它也是从输入开始的时候就有了输出。

（a）输入信号是方波信号输入信号波形图输出信号波形图分析：当输入信号是方波信号时，输出量一开始就有了，且从0开始逐渐增大，当达到1/2周期时，输出量达到了最大值；当方波信号为零电位时，输出量有不断减小，在一个周期时变为0，如此周期变化下去，因此可以得到输出信号的波形图如图所示。

3如图所示，开环系统是微分环节的开环系统。

当输入信号分别是方波信号和正弦信号时，分析输出信号的波形。

微分环节的输出为输入的倒数，反映了输入的变化趋势，因此，可以对输入的变化趋势经行预测。

微分环节可以对系统提前施加校正作用，以提高系统的灵敏度。

（a）输入信号是方波信号输入信号波形图输出信号波形图分析：当输入信号是方波信号时，在0~T/2和T/2~T内，方波信号处于高电位和点电位不发生变化，因此输出信号应该为0，在T/2时刻，输入信号从高电位变为地点为，变化率趋紧负无穷大，因此输出信号应该是接近负无穷大；同理，在T时刻，输出信号应该是接近正无穷大，但是由于系统的原因，都不会无穷大，而是以比较大的值，并且会有一定的时间延迟。

MATLAB 实验报告实验一信号在MATLAB中的表示1.画出下列信号的波形(1)y1(t)=e−2|t|源程序如下：t = -3:0.00001:3;ft = exp(-2*abs(t));plot(t,ft),grid on%axis([0,0.01,-1.2,1.2])title('y_1 (t)=e^(-2|t|)')运行结果如图：(2)y 2(n )=(0.9)n (−10≤n ≤10)源程序如下：n=-10:10;a=0.9;yn=a.^n;stem(n,yn);grid onaxis([-10,10,-1,3])title('y_2 (n)=〖(0.9)〗^n (-10≤n≤10)')运行结果如图：-3-2-1012300.10.20.30.40.50.60.70.80.91y 1 (t)=e (-2|t|)（3） y3(n)=e iπn3(−10≤n≤10)的实部源程序如下：n=-10:10;yn=exp(i*pi*(n/3));stem(n,real(yn));grid onaxis([-10,10,-1,1])title('（3）y_3 (n)=e^(iπn/3) (-10≤n≤10) 的实部') title('（3）y_3 (n)=e^(iπ n/3) (-10≤n≤10) 的实部') 运行结果如图：2.使用MATLAB命令，产生频率为1KHZ、占空比为50%的周期方波信号。

源程序如下：t = 0:0.000001:0.3;ft = square(2*1000*t,50);plot(t,ft),grid onaxis([0,0.01,-1.2,1.2])title('周期方波信号')运行结果如图：实验2 信号运算的MATLAB 实现1.画出下列信号的波形图（1）y1（t ）=3t^2+t+5源代码：t=-10:0.01:10;y1 = 3*(t.^2)+ t +5;plot(t,y1),grid ontitle('y1（t ）=3t^2+t+5')波形图：00.0010.0020.0030.0040.0050.0060.0070.0080.0090.01-1-0.8-0.6-0.4-0.20.20.40.60.81周期方波信号y1（t）=3t2+t+5（2）y2(t)=t*sin(5*pi*t)*u(t)源代码如下：t=--1:0.0001:5;y2 = t.*sin(5*pi*t).*heaviside(t);plot(t,y2),grid ontitle('y2(t)=t*sin(5*pi*t)*u(t)') 波形如图：y2(t)=t*sin(5*pi*t)*u(t)（3）f(t)=exp(-t).*sin(10*pi*t)+exp(-t/2).*sin(9*pi*t)源代码：t=-10:0.0001:-6;ft =exp(-t).*sin(10*pi*t)+exp(-t/2).*sin(9*pi*t);plot(t,ft),grid ontitle('ft =exp(-t).*sin(10*pi*t)+exp(-t/2).*sin(9*pi*t)') 波形如图：4ft =exp(-t).*sin(10*pi*t)+exp(-t/2).*sin(9*pi*t)2.已知f1（t）=sin,f2(t)=sin，使用MATLAB命令绘出f1（t）-f2（t）和f1（t）/f2（t）的波形图，其中f==1HZ.源程序如下：f = 1;t = 0:0.01:3/f;f1 = sin(2*pi*f*t);f2 = sin(2*pi*6*f*t);subplot(211),plot(t,f1-f2),grid on ,title('f1-f2');subplot(212),plot(t,f1./f2),grid on ,title('f1/f2');波形如图：f1-f215f1/f2实验三1.已知f(t)=u(t)-u(t-1),h(t)=u(t-1)-u(t-3),求时域卷积积分y(t)=f(t)*h(t)源程序如下：p = 0.01;t1=-1:p:1;f=ones(size(t1));t2=1:p:3;h=ones(size(t2));y=conv(f,h);t0=t1(1)+t2(1);t3=length(f)+length(h)-2;t=t0:p:(t3*p+t0);subplot(221)plot(t1,f)title('f(t)')xlabel('t');ylabel('f(t)');subplot(222)plot(t2,h)title('h(t)')xlabel('t');ylabel('h(t)');subplot(223)plot(t,y)g=get(gca,'position');g(3)=2.5*g(3);set(gca,'position',g)title('y(t)=f(t)*h(t)')xlabel('t');ylabel('y(t)'); 波形如图：2.某LTI 离散系统，设激励x(n)= (0.5).^n2,系统单位冲激响应h(n)=u(n+9)-u(n-6),试确定系统的零状态响应y(n)=x(n)*h(n)。

2015-2016数学实验练习题一、选择题1.清除Matlab工作空间（wordspace）变量的命令是（B ）A. clcB. clearC. clfD.delete2. 清除当前屏幕上显示的所有内容，但不清除工作空间中的数据的命令是（ A ）A. clcB. clearC. clfD.delete3. 用来清除图形的命令（ C ）A. clcB. clearC. clfD.delete4. 在MATLAB程序中，使命令行不显示运算结果的符号是（ A ）A. ;B. %C. #D. &5. 在MATLAB程序中，可以将某行表示为注释行的符号是（ B ）A. ;B. %C. #D. &6.在循环结构中跳出循环，执行循环后面代码的命令为 ( B )A. returnB. breakC. continueD. Keyboard7.在循环结构中跳出循环，但继续下次循环的命令为（ C ）A. returnB. breakC. continueD. Keyboard8. MATLAB中用于声明全局变量的关键字是( C )A. infB. symsC. globalD. function9. 用户可以通过下面哪项获得指令的使用说明（ A ）A. helpB. loadC. demoD. lookfor10．在MATLAB命令窗口中键入命令S=zoros(3)；可生成一个三行三列的零矩阵，如果省略了变量名S，MATLAB表现计算结果将用下面的哪一变量名做缺省变量名（ A ）A. ans;B. pi;C. NaN;D. Eps.11. 9/0的结果是（ B ）A. NAN；B. Inf；C. eps；D. 012．在MATLAB中程序或语句的执行结果都可以用不同格式显示，将数据结果显示为分数形式，用下面哪一条命令语句（ D ）A. format long；B. format long e；C. format bank；D. fromat rat13. 下列MATLAB命令中是构造1行3列的(-1,1)均匀分布随机矩阵的命令的是（D）A. randn(1,3)；B. rand(1,3)；C. ones(3)；D. 以上都不对14. 产生四维元素都为1矩阵的语句为（ A ）A. ones(4)B. eye(4)C. zeros(4)D. rand(4)15. 用round 函数对数组[2.48 6.39 3.93 8.52]取整，结果为 ( C )A. [2 6 3 8]B. [2 6 4 8]C. [2 6 4 9]D. [3 7 4 9]16. y=dsolve(‘Dy=1/(1+x^2)-2*y^2’,’y(0)=0’,’x ’); ezplot(y)的功能是（ A ）A. 求微分方程特解并绘图；B. 解代数方程；C. 求定积分；D.求微分方程通解.17. MATLAB 命令roots([1,0,0,-1])的功能是 （ D ）A. 产生向量[1,0,0,1]；B. 求方程310x +=的根；C. 求多项式31x -的值；D. 求方程310x -=的根。