matlab课后作业 北交大 钱满义
学生实验作业matlab
实验报告（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实验报告
北京交通⼤学现代控制理论matlab实验报告实验⼀利⽤MATLAB进⾏线性系统的模型转换及联结课程：现代控制理论姓名：王柬⽂学号：11212070班级：⾃动化1103指导教师：张勇标题 (1)⽬录 (2)⼀实验⽬的： (3)⼆实验原理： (3)三实验步骤： (6)四实验要求： (7)五实验内容： (8)1 (8)2 (18)3 (29)4 (36)5 (37)六实验感想： (40)⼀实验⽬的：1、学习系统状态空间模型的建⽴⽅法、了解状态空间模型与传递函数、零极点模型之间相互转换的⽅法；2、通过编程、上机调试，掌握系统状态空间模型与传递函数相互转换的⽅法。
3、通过编程、上机调试，掌握系统模型的联结⽅法。
⼆实验原理：⼀、连续系统（1）状态空间模型x Ax Buy Cx Du=+=+& （1.1）其中：nx R ∈是系统的状态向量，mu R ∈是控制输⼊，py R ∈是测量输出，A 是n n ?维状态矩阵、B 是n m ?维输⼊矩阵、C 是p n ?维输出矩阵、D 是直接转移矩阵。
在MATLAB 中，⽤（A,B,C,D ）矩阵组表⽰。
系统传递函数和状态空间模型之间的关系如式（1.2）所⽰。
1()()G s C sI A B D -=-+ （1.2）（2）传递函数模型11101110()(),()m m m m n n n n b s b s b s b num s H s m n den s a s a s a s a ----++++==≤++++L L 在MATLAB 中，直接⽤分⼦/分母的系数表⽰1010[,,,][,,,]m m n n num b b b den a a a --==L L（3）零极点增益模型1212()()()()()()()m n s z s z s z H s ks p s p s p ---=---L L在MATLAB 中，⽤[z, p, k]⽮量组表⽰，即1212[,,,];[,,,];[];m n z z z z p p p p k k ===L L ⼆、离散系统（1）传递函数模型11101110()m m m m n n n n b z b z b z b H z a z a z a z a ----++++=++++L L （2）零极点增益模型1212()()()()()()()m n z z z z z z H z kz p z p z p ---=---L L（3）状态空间模型(1)()()()()()x k Ax k Bu k y k Cx k Du k +=+=+三、三种模型间的转换表⽰状态空间模型和传递函数的MATLAB 函数。
信号与系统matlab课后作业-北京交通大学
信号与系统MATLAB平时作业学院：电子信息工程学院班级：：学号：教师：钱满义MATLAB 习题M3-1 一个连续时间LTI系统满足的微分方程为y ’’(t)+3y ’(t)+2y(t)=2x ’(t)+x(t)(1)已知x(t)=e -3t u(t)，试求该系统的零状态响应y zs (t)； (2)用lism 求出该系统的零状态响应的数值解。
利用(1)所求得的结果，比较不同的抽样间隔对数值解精度的影响。
解：(1) 由于''()3'()2()2'()(),0h t h t h t t t t δδ++=+≥则2()()()t t h t Ae Be u t --=+ 将()h t 带入原方程式化简得(2)()()'()2'()()A B t A B t t t δδδδ+++=+所以1,3A B =-=2()(3)()t t h t e e u t --=-+又因为3t ()()x t e u t -= 则该系统的零状态响应3t 23t 2t ()()()()(3)()0.5(6+5)()zs t t t y t x t h t e u t e e u t e e e u t ----=*=*-+=-- (2)程序代码 1、ts=0;te=5;dt=0.1;sys=tf([2 1],[1 3 2]);t=ts:dt:te;x=exp(-3*t).*(t>=0);y=lsim(sys,x,t)2、ts=0;te=5;dt=1;sys=tf([2 1],[1 3 2]);t=ts:dt:te;x=exp(-3*t).*(t>=0);y1=-0.5*exp(-3*t).*(exp(2*t)-6*exp(t)+5).*[t>=0];y2=lsim(sys,x,t)plot(t,y1,'r-',t,y2,'b--')xlabel('Time(sec)')legend('实际值','数值解')用lism求出的该系统的零状态响应的数值解在不同的抽样间隔时与（1）中求出的实际值进行比较将两种结果画在同一幅图中有图表 1 抽样间隔为1图表 2 抽样间隔为0.1图表 3 抽样间隔为0.01当抽样间隔dt减小时，数值解的精度越来越高，从图像上也可以看出数值解曲线越来越逼近实际值曲线，直至几乎重合。
MATLAB课后习题集(张亚运)-推荐下载
7、命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史
窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到 M 文件中。
8、 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件 有什么区别？
MATLAB 计算机建模与仿真
习题集
姓名： 张亚运 学号： 1121408028 班级：M11 电气工程及其自动化
对全部高中资料试卷电气设备，在安装过程中以及安装结束后进行高中资料试卷调整试验；通电检查所有设备高中资料电试力卷保相护互装作置用调与试相技互术通关，1系电过，力管根保线据护敷生高设产中技工资术艺料0不高试仅中卷可资配以料置解试技决卷术吊要是顶求指层，机配对组置电在不气进规设行范备继高进电中行保资空护料载高试与中卷带资问负料题荷试2下卷2，高总而中体且资配可料置保试时障卷，各调需类控要管试在路验最习；大题对限到设度位备内。进来在行确管调保路整机敷使组设其高过在中程正资1常料中工试，况卷要下安加与全强过，看度并22工且22作尽22下可22都能22可地护以缩1关正小于常故管工障路作高高；中中对资资于料料继试试电卷卷保破连护坏接进范管行围口整，处核或理对者高定对中值某资，些料审异试核常卷与高弯校中扁对资度图料固纸试定，卷盒编工位写况置复进.杂行保设自护备动层与处防装理腐置，跨高尤接中其地资要线料避弯试免曲卷错半调误径试高标方中高案资等，料，编试要5写、卷求重电保技要气护术设设装交备备置底4高调、动。中试电作管资高气，线料中课并敷3试资件且、设卷料中拒管技试试调绝路术验卷试动敷中方技作设包案术，技含以来术线及避槽系免、统不管启必架动要等方高多案中项；资方对料式整试，套卷为启突解动然决过停高程机中中。语高因文中此电资，气料电课试力件卷高中电中管气资壁设料薄备试、进卷接行保口调护不试装严工置等作调问并试题且技，进术合行，理过要利关求用运电管行力线高保敷中护设资装技料置术试做。卷到盒程护处中装，高置当中高不资中同料资电试料压卷试回技卷路术调交问试叉题技时，术，作是应为指采调发用试电金人机属员一隔，变板需压进要器行在组隔事在开前发处掌生理握内；图部同纸故一资障线料时槽、，内设需，备要强制进电造行回厂外路家部须出电同具源时高高切中中断资资习料料题试试电卷卷源试切，验除线报从缆告而敷与采设相用完关高毕技中，术资要资料进料试行，卷检并主查且要和了保检解护测现装处场置理设。备高中资料试卷布置情况与有关高中资料试卷电气系统接线等情况，然后根据规范与规程规定，制定设备调试高中资料试卷方案。
北京交通大学《信号与系统》 课后matlab作业
Matlab课后作业1.M2-1（1）Matlab程序：t=-5:0.01:5;x=(t>0)-(t>2);plot(t,x);axis([-5,5,-2,2]);仿真结果：（8）Matlab程序：t=-10:0.01:10;pi=3.14;x=sin(pi*t)./(pi*t).*cos(30*t);plot(t,x);仿真结果：M2-2Matlab程序：t=-2:0.001:2;x=(t>-1)-(t>0)+2*tripuls(t-0.5,1,0); plot(t,x);axis([-2,2,-2,2]);仿真结果：M3-3（1）function yt=f(t)yt=t.*(t>0)-t.*(t>=2)+2*(t>=2)-3*(t>3)+(t>5); （2）Matlab程序：t=-10:0.01:11;subplot(3,1,1);plot(t,f(t));title('x(t)');axis([-1,6,-2,3]);subplot(3,1,2);plot(t,f(0.5*t));axis([-1,11,-2,3]);title('x(0.5t)');subplot(3,1,3);plot(t,f(2-0.5*t));title('x(2-0.5t)');axis([-9,5,-2,3]);仿真结果：M2-9(1)Matlab程序：k=-4:7;x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; stem(k,x);仿真结果：（2）Matlab程序：k=-12:21;x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; N=length(x);y=zeros(1,3*N-2);y(1:3:end)=x;stem(k,y);仿真结果：Matlab程序：k=-1:3;x=[0,0,-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; x1=x(1:3:end);stem(k-1,x1);仿真结果：（3）Matlab程序：k=-6:5;x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; stem(k,x);仿真结果：程序>> k=-2:9;>> x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; >> stem(k,x);结果程序>> k=-4:7;>> x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; >> xk=fliplr(x);>> k1=-fliplr(k);>> stem(k1,xk);结果M3-1（1）程序>> ts=0;te=5;dt=0.01; >> sys=tf([2 1],[1 3 2]); >> t=ts:dt:te;>> x=exp(-3*t).*(t>=0); >> y=lsim(sys,x,t);>> plot(t,y);>> xlabel('Time(sec)') >> ylabel('y(t)')结果（2）程序>> ts=0;te=5;dt=0.0001; >>sys=tf([2 1],[1 3 2]); >>t=sys:dt:te;>>x=exp(-3*t).*(t>=0); >>y=lsim(sys,x,t);>>plot(t,y);>>xlabel('Time(sec)') >>ylabel('y(t)')结果M3-4>> x=[0.85,0.53,0.21,0.67,0.84,0.12]; >> k1=-2:3;>> h=[0.68,0.37,0.83,0.52,0.71];>> k2=-1:3;>> y=conv(x,h);>> k=(k1(1)+k2(1)):(k1(end)+k2(end)); >> stem(k,y)结果M6-1(1)>> num=[16 0 0];>> den=[1 5.6569 816 2262.7 160000]; >> [r,p,k]=residue(num,den)得r =0.0992 - 1.5147i0.0992 + 1.5147i-0.0992 + 1.3137i-0.0992 - 1.3137ip =-1.5145 +21.4145i-1.5145 -21.4145i-1.3140 +18.5860i-1.3140 -18.5860ik =[]所以可得 X(s)=j s j j s j j s j 5860.183140.13137.10992.05860.183140.13137.10992.04145.215145.15147.10992.021.4145j -1.5145s j 5147.1-0992.0++--+-++-++++++x(t)=3.0108e-1.5145tcos(21.4145t-1.5054)u(t)+2.635e-1.314tcos(18.586t+1.6462)u(t ) (2)X(s)=)2552^)(5(2^+++s s s s解:>> num=[1 0 0 0];den=conv([1 5],[1 5 25]);[r,p,k]=residue(num,den)[angle,mag]=cart2pol(real(r),imag(r))得r =-5.0000 + 0.0000i-2.5000 - 1.4434i-2.5000 + 1.4434ip =-5.0000 + 0.0000i-2.5000 + 4.3301i-2.5000 - 4.3301ik =1angle =3.1416-2.61802.6180mag =5.00002.88682.8868所以X(s)=3301.45.24434.15.23301.45.24434.15.25s 5.0-1j s j j s j +++-+-+--+++x(t)=δ(t)+5e-5tu(t)+5.7736e-2.5tcos(4.3301t-2.618)u(t)M6-2程序>> t=0:0.1:10;>> y1=(2.5*exp(-t)-1.5*exp(-3*t)).*(t>=0);>> y2=((1/3)+2*exp(-t)-(5/6)*exp(-3*t)).*(t>=0);>> y=((1/3)+(9/2)*exp(-t)-(7/3)*exp(-3*t)).*(t>=0);>> plot(t,y1,'r-',t,y2,'g--',t,y,'b-')>> xlabel('Time');>> legend('零输入响应','零状态响应','完全响应')结果M6-5>> num=[1 2];>> den=[1 2 2 1];>> sys=tf(num,den);>> pzmap(sys)>> num=[1 2];den=[1 2 2 1];[r,p,k]=residue(num,den) [angle,mag]=cart2pol(real(r),imag(r))1.0000 + 0.0000i-0.5000 - 0.8660i-0.5000 + 0.8660ip =-1.0000 + 0.0000i-0.5000 + 0.8660i-0.5000 - 0.8660ik =[]angle =-2.09442.0944mag =1.00001.00001.0000所以H(s)=866.05.0866.05.0866.05.0866.05.01s 1j s j j s j +++-+-+--++系统冲激响应h(t)=e-tu(t)+2e-0.5tcos(0.866t-2.0944)u(t)>> num=[1 2];>> den=conv([1 0],[1 2 2 1]);>> [r,p,k]=residue(num,den)r =-1.0000 + 0.0000i-0.5000 + 0.8660i-0.5000 - 0.8660i2.0000 + 0.0000ip =-1.0000 + 0.0000i-0.5000 + 0.8660i-0.5000 - 0.8660i0.0000 + 0.0000ik =[][angle,mag]=cart2pol(real(r),imag(r))angle =3.14162.0944-2.0944mag =1.00001.00001.00002.0000所以Y(s)=s j s j j s j 2866.05.0866.05.0866.05.0866.05.0-1s 1-+++--+-++++ 系统阶跃响应y(t)=e-tu(t)+2e-0.5tcos(0.866t+2.0944)u(t)因为系统的冲激响应h(t)=e-tu(t)-1.00001e-0.5tcos(0.866t)u(t)+1.73205e-0.5tsin(0.866t)u(t) 所以系统的频率响应H(j ω)=5.0)866.0(866.05.0)866.0(866.05.0)866.0(5.05.0)886.0(5.01j 1j j j j ++--+--++-+++++ωωωωω。
MATLAB课后习题解答
信号与系统MATLAB平时作业第二章（1）代码：t=-4:0.001:4;>> T=2;>> xt=rectpuls(t-1,T);>> plot(t,xt);>> axis([-4,4,-0.5,1.5]) 图形：（2）代码：t=sym('t');>> y=heaviside(t); >> ezplot(y,[-1,1]); >> grid on;>> axis([-1,1,-0.2,1.2]) 图形：(3)代码：> A=10;a=-1;B=5;b=-2;>> t=0:0.001:10;>> xt=A*exp(a*t)-B*exp(b*t); >> plot(t,xt);图形：（4）代码：t=sym('t');y=t*heaviside(t);ezplot(y,[-1,3]);grid on ;axis([-1 3 -0.1 3.1])图形：（5）代码：xt=abs(A*sin(w0*t+phi)); >> plot(t,xt);>> A=2;w0=10*pi;phi=pi/3;图形：（6）代码：> A=1;w0=1;B=1;w1=2*pi;>> t=0:0.001:20;>> xt=A*cos(w0*t)+B*sin(w1*t); >> plot(t,xt);图形：（7）代码：A=4;a=-0.5;w0=2*pi;>> t=0:0.001:10;>> xt=A*exp(a*t).*cos(w0*t); >> plot(t,xt);图形：（8）代码：w0=30;t=-15:0.001:15;xt=sinc(t/pi).*cos(w0*t); plot(t,xt);图形：M2-2（1）xt程序：t=-1:0.001:1;x=4*t.*[t>=0]-4*t.*[t>=1/2]-4*t.*[t>=1/2]+4*[t>=1/2]+4*t.*[t>=1]-4*[t>=1 ]+[t>=-1]-[t>=0];>> plot(t,x);图形：(2)程序：t=-1:0.001:1;x=(4*t.*[t>=0]-4*t.*[t>=1/2]-4*t.*[t>=1/2]+4*[t>=1/2]+4*t.*[t>=1]-4*[t>= 1]+[t>=-1]-[t>=0]).*cos(50*t);>> plot(t,x);图形：M2-3（1）yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5);（2）t=0:0.001:6;yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5);subplot(3,1,1)plot(t,yt)title('x(t)')axis([0,6,-2,2])t=0:0.001:12;yt=((0.5)*t).*(t>=0&t<=4)+2*(t>=4&t<=6)-1*(t>=6&t<=10); subplot(3,1,2)plot(t,yt);axis([0,12,-2,2]);title('(0.5)*t');t=-6:0.001:6;yt=((2-(0.5)*t)).*(t>=0&t<=4)+2*(t>=-2&t<=0)-1*(t>=-6&t<=-2); subplot(3,1,3)plot(t,yt);axis([-6,6,-2,2]);title('2-(0.5)*t');图形：M2-4奇分量程序：t=-2:0.001:2;x1t=(t.*[t>=-1]-t.*[t>=0]+[t>=-1]-[t>=0]+t.*[t<=1]-t.*[t<=0]-[t<=1]+[t<=0]) /2;>> plot(t,x1t);图形：偶分量：程序：t=-2:0.001:2;>>x2t=(t.*[t>=-1]-t.*[t>=0]+[t>=-1]-[t>=0]-t.*[t<=1]+t.*[t<=0]+[t<=1]-[t<=0]) /2;>> plot(t,x2t);图形：M2-5 （1）function [x,k]=impseq(0,-50,50) k=[-50:50];x=[k==0];stem(k,x)图形：（2）k=[-50:50];>> uk=[zeros(1,50),ones(1,51)];>> stem(k,uk)图形：（3）function [f,k]=impseq(0,-50,50)；k=[-50:50];x=10*(1/2).^k.*[k>=0];>> stem(k,x);图形：（4）function [f,k]=impseq(k0,-50,50) k=[-10:10];x1=[(k+2)>=0];x2=[(k-6)>=0];x=x1-x2;stem(k,x);图形：（5）function [f,k]=impseq(k0,-50,50) k=[-50:50];>> x=k.*[k>=0];>> stem(k,x);图形：（6）k=-50:50;xk=(5*0.8.^k).*cos(0.9*pi*k);>> stem(k,xk);图形：M2-6 Ω0=0.1*pi;k=-50:50;>> xk=sin(0.1*pi*k);>> stem(k,xk);图形：Ω0=0.5*pi;k=-50:50;xk=sin(0.5*pi*k);stem(k,xk);图形：Ω0=0.9*pi;k=-50:50;xk=sin(0.9*pi*k); stem(k,xk);图形：Ω0=1.1*pi;k=-50:50;xk=sin(1.1*pi*k); stem(k,xk);图形：Ω0=1.5*pi;k=-50:50;xk=sin(1.5*pi*k); stem(k,xk);图形：Ω0=1.9*pi;k=-50:50;xk=sin(1.9*pi*k); stem(k,xk);图形：结论：随着Ω0的变动，波形也随之而成周期性变换M2-7 （1）程序：t=-2:0.001:2;x1t=cos(6*pi*t);plot(t,x1t);hold onk=-2:2;x1k=cos(0.6*pi*k);stem(k,x1k,'r');hold off图形：（2）程序：t=-1:0.001:1;x1t=cos(14*pi*t); plot(t,x1t);hold onk=-1:1;x1k=cos(1.4*pi*k); stem(k,x1k,'r'); hold off图形：（3）程序：t=-0.5:0.001:0.5;x1t=cos(26*pi*t);plot(t,x1t);hold onk=-0.5:0.5;x1k=cos(2.6*pi*k);stem(k,x1k,'r');hold off图形：M2-8 （1）用square函数：程序：k=0:40;A=1;P=0.6;y=A*square(P*k);stem(k,y)axis([k(1)-1 k(end)+1 -(A+1) (A+1)]); 图形：（2）程序：k=0:40;x2k=sawtooth(0.1*pi*k,0.5);stem(k,x2k);图形：M2-9(1)k=-4:7;xk=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1];stem(k,xk)图形：(2)k=-12:21x=[-3,0,0,-2,0,0,3,0,0,1,0,0,-2,0,0,-3,0,0,-4,0,0,2,0,0,-1,0,0,4,0,0,1,0,0,-1]; subplot(2,1,1)stem(k,x)title('3倍内插)t=-1:2;y=[-2,-2,2,1];subplot(2,1,2)stem(t,y)title('3倍抽取’)axis([-3,4,-4,4])图形：（3）k=-4:7;x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; subplot(2,1,1)stem(k+2,x)title('x[k+2]')subplot(2,1,2)stem(k-4,x)title('x[k-4]')图形：（4）k=-4:7;x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; stem(-fliplr(k),fliplr(x))title('x[-k]')图形：第三章M3-1(1)代码：ts=0;te=5;dt=0.01;>> sys=tf([2 1],[1 3 2]);>> t=ts:dt:te;>> x=exp(-3*t);>> y=lsim(sys,x,t);>> plot(t,y);>> xlabel('Time(sec)');>> ylabel('y(t)');图形：（2）代码：ts=0;te=5;dt=1; sys=tf([2 1],[1 3 2]); t=ts:dt:te;x=exp(-3*t);y=lsim(sys,x,t); >> plot(t,y); xlabel('Time(sec)'); ylabel('y(t)');图形：M3-2 代码：ts=0;te=5;dt=0.01;>> sys=tf([6],[1 3 6]);>> t=ts:dt:te;>> y=impulse(sys,t);>> plot(t,y);>> xlabel('Time(sec)');>> ylabel('h(t)');图形：（3）代码：ts=0;te=5;dt=0.01; sys=tf([6],[1 3 6]); t=ts:dt:te;y=step(sys,t);plot(t,y);xlabel('Time(sec)'); ylabel('u(t)')图形：M3-3 (1)代码：ts=0;te=5;dt=1;sys=tf([1],[1 0.2 1]);t=ts:dt:te;y=step(sys,t);plot(t,y);xlabel('Time(sec)');ylabel('u(t)')图形：（2）代码：ts=0;te=5;dt=1; sys=tf([1],[1 1 1]); t=ts:dt:te;y=step(sys,t);plot(t,y);xlabel('Time(sec)'); ylabel('u(t)')图形：（3）代码：ts=0;te=5;dt=1;sys=tf([1],[1 2 1]);t=ts:dt:te;y=step(sys,t);plot(t,y);xlabel('Time(sec)');ylabel('u(t)')图形：M3-4 代码：x=[0.85,0.53,0.21,0.67,0.84,0.12];kx=-2:3;h=[0.68,0.37,0.83,0.52,0.71];kh=-1:3;y=conv(x,h);k=kx(1)+kh(1):kx(end)+kh(end);stem(k,y,'r');图形：M3-6h(t)=x(t)x(t)=u(t)*u(t)+u(t-1)*u(t-1)+2u(t)*u(t-1)=r(t)-2r(t-1)+r(t-2)y(t)=x(t) h(t)=1/2*t^2.*u(t)-3/2*(t-1)^2.*u(t-1)-1/2*(t-3)^2.*u(t-3)0, t<0t^2/2, 0<=t<1 Y(t) = -t^2+3*t-3/2 1<=t<2( t-3)^2/2 2<=t<30 t>=3(b)T=0.1;k=-1:T:4;f1=1*((k>=0)&(k<=1));f2=tripuls(k-1,2);y=T*conv(f1,f2);tmin=-2;tmax=8;t1=tmin:0.1:tmax;plot(t1,y)grid on图形：T=0.1：图形：T=0.01T=0.001：M3-7 代码：b=[1 0.42 -0.19];a=[0.31 0.68];k=0:20;x=(0.6).^k.*[k>=0];y=filter(b,a,x);stem(k,y);图形：M3-8 代码：b=[1 0.7 -0.45 -0.6];a=[0.8 -0.44 0.36 0.02];k=0:30;h=impz(b,a,k)stem(k,h);h =1.25001.5625-0.2656-1.6305-0.81630.29140.56840.2019-0.1520-0.1887-0.04040.06650.05950.0038-0.0263-0.01770.00200.00970.0049-0.0017-0.0034-0.00120.00090.00110.0002-0.0004-0.0004-0.00000.00020.0001-0.0000 图形：第四章M4-1周期矩形信号：代码：n=-20:20;X=-j*1/2*sin(n/2*pi).*sinc(n/2);subplot(2,1,1);stem(n,abs(X));title('幅度谱')xlabel('nw');subplot(2,1,2);stem(n,angle(X));title('相位谱')图形：三角波信号：代码：>> n=-20:20;X=sinc(n)-0.5*((sinc(n/2)).^2); subplot(2,1,1);stem(n,abs(X));title('幅度谱')xlabel('nw');subplot(2,1,2);stem(n,angle(X));title('相位谱')图形：M4-6 (4)代码：x=[1,2,3,0,0];X=fft(x,5);subplot(2,1,1);m=0:4;stem(m,real(X)); title('X[m]实部') subplot(2,1,2); >> stem(m,imag(X)); title('X[m]虚部');图形：M4-7 (3)代码：k=0:10;x=0.5.^k;subplot(3,1,1);stem(k,x)title('x[k]')X=fft(x,10);subplot(3,1,2);m=0:9;stem(m,real(X));title('X[m]实部')>> subplot(3,1,3);>> stem(m,imag(X));title('X[m]虚部')图形：M5-2代码：t=0:0.05:2.5;T=1;xt1=rectpuls(t-0.5,T);subplot(2,2,1)plot(t,xt1)title('x(t1)')axis([0,2.5,0,2])xt2=tripuls(t-1,2);subplot(2,2,2)plot(t,xt2)title('x(t2)')axis([0,2.5,0,2])xt=xt1+xt2.*cos(50*t);subplot(2,2,[3,4])plot(t,xt)title('x(t)')figure;b=[10000];a=[1,26.131,341.42,2613.1,10000];[H,w]=freqs(b,a,w);subplot(2,1,1)plot(w,abs(H));set(gca,'xtick',[0 1 2 3 4 5 6 7 8]);set(gca,'ytick',[0.985 0.99 0.995 1 1.005]);grid; title('幅度曲线')subplot(2,1,2)plot(w,angle(H));set(gca,'xtick',[0 1 2 3 4 5 6 7 8])set(gca,'ytick',[-2 -1.5 -1 -0.5 0]);grid;title('相位曲线')figure;sys=tf([10000],[1 26.131 341.42 2613.1 10000]); yt1=lsim(sys,xt,t);subplot(2,1,1);plot(t,yt1);title('y(t1)')yt2=lsim(sys,xt.*cos(50*t),t);subplot(2,1,2);plot(t,yt2);title('y(t2)');图形：第六章 M6-1已知连续时间信号的s 域表示式如下，使用residue 求出X(s)的部分分式展开式，并写出x(t)的实数形式表达式。
matlab课后习题与答案
习题二1.如何理解“矩阵是MATLAB最基本的数据对象”？答：因为向量可以看成是仅有一行或一列的矩阵，单个数据（标量）可以看成是仅含一个元素的矩阵，故向量和单个数据都可以作为矩阵的特例来处理。
因此，矩阵是MATLAB最基本、最重要的数据对象。
2.设A和B是两个同维同大小的矩阵，问：(1)A*B和A.*B的值是否相等？答：不相等。
(2)A./B和B.\A的值是否相等？答：相等。
(3)A/B和B\A的值是否相等？答：不相等。
(4)A/B和B\A所代表的数学含义是什么？答：A/B等效于B的逆右乘A矩阵，即A*inv(B)，而B\A等效于B矩阵的逆左乘A矩阵，即inv(B)*A。
3.写出完成下列操作的命令。
(1)将矩阵A第2~5行中第1, 3, 5列元素赋给矩阵B。
答：B=A(2:5,1:2:5); 或B=A(2:5,[1 3 5])(2)删除矩阵A的第7号元素。
答：A(7)=[](3)将矩阵A的每个元素值加30。
答：A=A+30;(4)求矩阵A的大小和维数。
答：size(A);ndims(A);(5)将向量t的0元素用机器零来代替。
答：t(find(t==0))=eps;(6)将含有12个元素的向量x转换成34矩阵。
答：reshape(x,3,4);(7)求一个字符串的ASCII码。
答：abs(‘123’); 或double(‘123’);(8) 求一个ASCII 码所对应的字符。
答：char(49);4. 下列命令执行后，L1、L2、L3、L4的值分别是多少？A=1:9;B=10-A;...L1=A==B;L2=A<=5;L3=A>3&A<7;L4=find(A>3&A<7);答：L1的值为[0, 0, 0, 0, 1, 0, 0, 0, 0]L2的值为[1, 1, 1, 1, 1, 0, 0, 0, 0]L3的值为[0, 0, 0, 1, 1, 1, 0, 0, 0]L4的值为[4, 5, 6]5. 已知23100.7780414565532503269.5454 3.14A -⎡⎤⎢⎥-⎢⎥=⎢⎥⎢⎥-⎣⎦完成下列操作：(1) 取出A 的前3行构成矩阵B ，前两列构成矩阵C ，右下角32⨯子矩阵构成矩阵D ，B 与C 的乘积构成矩阵E 。
MATLAB结课大作业
一、判断系统稳定性的方法，并举例说明。
方法一：用Nyquist稳定判据判断系统的稳定性Nyquist稳定判据：若想使得闭环系统稳定，则开环系统G（S）H（S）的Nyquist曲线逆时针绕临界点（-1，j0）点的圈数R必需等于G（S）H（S）（系统的开环传函）位于S的右半平面开环极点数P。
即：Z=P-R Z=0 稳定；Z≠0 不稳定，Z为闭环正实部根的个数。
方法二：用Bode图判断系统的稳定性函数调用格式为：margin( )或[Gm Pm wcp wcg]=margin(G)对于最小相位系统：当相角裕度P m(γ)>0o 或幅值裕度G m(h) >1时，表示系统稳定当相角裕度P m(γ)<0o 或幅值裕度G m(h) <1时，表示系统不稳定幅值裕度G m(h)、相角裕度P m(γ)越大，系统稳定程度越好。
在使用时，G m(h)、P m(γ)是成对使用的，有时仅使用一个裕度指标P m（γ）。
方法三：用代数稳定判据法判断系统的稳定性(1)系统数学模型为传递函数形式G(S)=tf(num,den):执行语句：roots(G.den{1});注：“{}”表示维数(2)系统数学模型为零极点增益形式G(S)=zpk(z,p,k);执行语句：G.p{1};(3)系统数学模型为状态空间形式G(S)=ss(A,B,C,D);执行语句：eig(G.A);注：eig（）表示计算系统的极点方法四：用根轨迹法判断系统的稳定性若根轨迹在参数取值过程中,部分在左半平面,部分在右半平面,则系统的稳定性与可变参数的取值有关。
函数命令调用格式：[k poles]=rlocfind(G)方法五：用单位阶跃响应曲线判定系统稳定性例：已知系统的开环传函为：5（S+2）G（S）= ----------------------------（S+10）（S³+3S²+2S+5）判断系统的稳定性解：方法一：用Nyquist稳定判据判断系统的稳定性>>G=tf（5*[1 2],conv([1 10],[1 3 2 5])）;>>roots(G.den{1})Ans=-10.0000-2.9042-0.0479+1.3112i-0.0479-1.3112i>>nyquist(G)执行以上程序得到nyquist图可以看出，nyquist曲线包围临界点的圈数R=0。
MATLAB教程2012a习题解答1-7章完整版-张志涌-北航
hut =
(3*cos(t))/10 - 3/(10*exp(3*t)) + sin(t)/10
（2）变换法（复验）
syms s;
HU=laplace(h,t,s)*laplace(u,t,s);
huL=simple(ilaplace(HU,s,t))%拉氏变换及反变换
syms a positive
solve(x^2-a*x+a^2)
ans =
a/2 + (3^(1/2)*a*i)/2
a/2 - (3^(1/2)*a*i)/2
a =
在此，
〖目的〗
理解准确符号数值的创建法。
高精度误差的观察。
〖解答〗
（1）x=7/3
x=7/3;a=x,b=sym(x),c=sym(x,'d'),d=sym('7/3'),
d0_=limit(d,t,0,'left') %求dy/dt|t=0-
dpi_2=limit(d,t,pi/2) %求dy/dt|t=pi/2
y =
abs(sin(t))
d =
sign(sin(t))*cos(t)
d0_ =
-1
dpi_2 =
0
〖目的〗
符号积分的解析解和符号数值解。
符号计算和数值计算的相互校验。
ans =
a
symvar(sym('z*exp(j*th)'),1)
ans =
z
（
（
〖目的〗
体验变量限定假设的影响
〖解答〗
（1）求三阶方程 正实根
matlab课后练习习题及答案详解
matlab课后习题及答案详解第1章MATLAB概论与其余计算机语言对比较，MATLAB语言突出的特色是什么？MATLAB拥有功能强盛、使用方便、输入简捷、库函数丰富、开放性强等特色。
MATLAB系统由那些部分构成？MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。
安装MATLAB时，在选择组件窗口中哪些部分一定勾选，没有勾选的部分此后怎样补安装？在安装MATLAB时，安装内容由选择组件窗口中个复选框能否被勾选来决定，能够依据自己的需要选择安装内容，但基本平台（即MATLAB选项）一定安装。
第一次安装没有选择的内容在补安装时只要依据安装的过程进行，不过在选择组件时只勾选要补装的组件或工具箱即可。
MATLAB操作桌面有几个窗口？怎样使某个窗口离开桌面成为独立窗口？又怎样将离开出去的窗口从头搁置到桌面上？在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是封闭窗口的Close按钮，一个是能够使窗口成为独立窗口的Undock 按钮，点击Undock按钮就能够使该窗口离开桌面成为独立窗口，在独立窗口的view菜单中选择Dock,,菜单项就能够将独立的窗口从头防备的桌面上。
怎样启动M文件编写/调试器？在操作桌面上选择“成立新文件”或“翻开文件”操作时，M文件编写/调试器将被启动。
在命令窗口中键入edit命令时也能够启动M文件编写/调试器。
储存在工作空间中的数组能编写吗？怎样操作？储存在工作空间的数组能够经过数组编写器进行编写：在工作空间阅读器中双击要编写的数组名翻开数组编写器，再选中要改正的数据单元，输入改正内容即可。
命令历史窗口除了能够察看前方键入的命令外，还有什么用途？命令历史窗口除了用于查问从前键入的命令外，还能够直接履行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。
怎样设置目前目录和搜寻路径，在目前目录上的文件和在搜寻路径上的文件有什么差别？目前目录能够在目前目录阅读器窗口左上方的输入栏中设置，搜寻路径能够经过选择操作桌面的file菜单中的SetPath菜单项来达成。
MATLAB课后习题集附标准答案
第2章 MATLAB概论1、与其他计算机语言相比较，MATLAB语言突出的特点是什么？答：起点高、人机界面适合科技人员、强大而简易的作图功能、智能化程度高、功能丰富，可扩展性强.2、MATLAB系统由那些部分组成？答：开发环境、MATLAB数学函数库、MATLAB语言、图形功能、应用程序接口3、安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？答：在安装MATLAB时，安装内容由选择组件窗口中各复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装.第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可.矚慫润厲钐瘗睞枥庑赖。
4、MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？聞創沟燴鐺險爱氇谴净。
答：在MATLAB操作桌面上有五个窗口，在每个窗口的右下角有两个小按钮，一个是关闭窗口的Close 按钮，一个是可以使窗口称为独立的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面称为独立窗口，在独立窗口的view菜单中选择Dock，菜单项就可以将独立的窗口重新防止的桌面上.残骛楼諍锩瀨濟溆塹籟。
5、如何启动M文件编辑/调试器？答：在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动.在命令窗口中键入edit命令时也可以启动M文件编辑/调试器.酽锕极額閉镇桧猪訣锥。
6、存储在工作空间中的数组能编辑吗？如何操作？答：存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可.彈贸摄尔霁毙攬砖卤庑。
7、命令历史窗口除了可以观察前面键入的命令外，还有什么用途？答：命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中.謀荞抟箧飆鐸怼类蒋薔。
终版matlab习题及答案.doc
1/2
>> limit(symsum(1/((2*k)^2-1),k,1,n),n,inf)
ans =
1/2
（2）
>> limit(n*symsum(1/(n^2+k*pi),k,1,n),n,inf)
ans =
1
第十一题
（1）
>> syms a t;x=a*(cos(t)+t*sin(t));y=a*(sin(t)-t*cos(t));
>> P=y*x^3+exp(y);Q=x*y^3+x*exp(y)-2*y;
>> ds=[diff(x,t);diff(y,t)];I=int([P Q]*ds,t,0,pi)
I =
-2/15*c*(-2*c^4+15*b^4)/b^4/a
第十二题
>> syms a b c d e;A=vander([a b c d e])
ans =
-(-2*sin(t)-t*cos(t))/sin(t)*cos(t)
>> f=diff(y,t,2)/diff(x,t,2);subs(f,t,sym(pi)/3)
ans =
3/8-1/24*pi*3^(1/2)
第九题
>> syms x y t
>> s=int(exp(-t^2),t,0,x*y);
2.用MATLAB语句输入矩阵 和
3．假设已知矩阵 ，试给出相应的MATLAB命令，将其全部偶数行提取出来，赋给 矩阵，用 命令生成 矩阵，用上述命令检验一下结果是不是正确。
4．用数值方法可以求出 ，试不采用循环的形式求出和式的数值解。由于数值方法是采用double形式进行计算的，难以保证有效位数字，所以结果不一定精确。试采用运算的方法求该和式的精确值。
MATLAB课后习题资料
第一部分 MATLAB 运算基础1. 先求下列表达式的值，然后显示MATLAB 工作空间的使用情况并保存全部变量。
(1) 0122sin851z e=+ (2) 221ln(1)2z x x =++，其中2120.455i x +⎡⎤=⎢⎥-⎣⎦ (3) 0.30.330.3sin(0.3)ln , 3.0, 2.9,,2.9,3.022a a e e az a a --+=++=--(4) 2242011122123t t z t t t t t ⎧≤<⎪=-≤<⎨⎪-+≤<⎩，其中t =0:0.5:2.52. 已知：1234413134787,2033657327A B --⎡⎤⎡⎤⎢⎥⎢⎥==⎢⎥⎢⎥⎢⎥⎢⎥-⎣⎦⎣⎦求下列表达式的值：(1) A+6*B 和A-B+I （其中I 为单位矩阵） (2) A*B 和A.*B (3) A^3和A.^3 (4) A/B 及B\A(5) [A,B]和[A([1,3],:);B^2]3. 设有矩阵A 和B123453166789101769,111213141502341617181920970212223242541311A B ⎡⎤⎡⎤⎢⎥⎢⎥-⎢⎥⎢⎥⎢⎥⎢⎥==-⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎢⎥⎣⎦⎣⎦(1) 求它们的乘积C 。
(2) 将矩阵C 的右下角3×2子矩阵赋给D 。
(3) 查看MATLAB 工作空间的使用情况。
4. 完成下列操作：(1) 求[100,999]之间能被21整除的数的个数。
(2) 建立一个字符串向量，删除其中的大写字母。
第二部分 MATLAB 矩阵分析与处理1. 设有分块矩阵33322322E R A O S ⨯⨯⨯⨯⎡⎤=⎢⎥⎣⎦，其中E 、R 、O 、S 分别为单位矩阵、随机矩阵、零矩阵和对角阵，试通过数值计算验证22E R RS A OS +⎡⎤=⎢⎥⎣⎦。
2. 产生5阶希尔伯特矩阵H 和5阶帕斯卡矩阵P ，且求其行列式的值Hh 和Hp 以及它们的条件数Th 和Tp ，判断哪个矩阵性能更好。
MATLAB教程2012a习题解答1-7章完整版-张志涌-北航
●MATLAB R2012a●课后习题答案全解●第一章基础准备及入门习题1及解答⏹ 1.数字1.5e2，1.5e3 中的哪个与1500相同吗？〖解答〗1.5e3⏹ 2.请指出如下5个变量名中，哪些是合法的？abcd-2 xyz_3 3chan a变量ABCDefgh〖解答〗2、5是合法的。
⏹ 3.在MATLAB环境中，比1大的最小数是多少？〖解答〗1+eps⏹ 4.设a = -8 , 运行以下三条指令，问运行结果相同吗？为什么？w1=a^(2/3)w2=(a^2)^(1/3)w3=(a^(1/3))^2〖解答〗（1）不同。
具体如下w1=a^(2/3) %仅求出主根w2=(a^2)^(1/3) %求出(-8)^2的主根 w3=(a^(1/3))^2%求出(-8)主根后再平方w1 = -2.0000 + 3.4641iw2 = 4.0000 w3 =-2.0000 + 3.4641i（2）复数的多方根的，下面是求取全部方根的两种方法： （A ）根据复数方根定义a=-8;n=2;m=3;ma=abs(a);aa=angle(a); for k=1:m%m 决定循环次数 sa(k)=(aa+2*pi*(k-1))*n/m;%计算各根的相角 endresult=(ma^(2/3)).*exp(j*sa) %计算各根result =-2.0000 + 3.4641i 4.0000 - 0.0000i -2.0000 - 3.4641i（B ）利用多项式023=-a r 求根p=[1,0,0,-a^2]; r=roots(p) r =-2.0000 + 3.4641i -2.0000 - 3.4641i 4.0000⏹ 5.指令clear, clf, clc 各有什么用处？〖解答〗 clear 清除工作空间中所有的变量。
clf 清除当前图形。
clc 清除命令窗口中所有显示。
⏹ 6.以下两种说法对吗？（1）“MATLAB 进行数值的表达精度与其指令窗中的数据显示精度相同。
matlab课后习题解答第二章
matlab课后习题解答第⼆章第2章符号运算习题2及解答1 说出以下四条指令产⽣的结果各属于哪种数据类型，是“双精度”对象，还是“符号”符号对象3/7+; sym(3/7+; sym('3/7+'); vpa(sym(3/7+)〖⽬的〗不能从显⽰形式判断数据类型，⽽必须依靠class指令。
〖解答〗c1=3/7+c2=sym(3/7+c3=sym('3/7+')c4=vpa(sym(3/7+)Cs1=class(c1)Cs2=class(c2)Cs3=class(c3)Cs4=class(c4)c1 =c2 =37/70c3 =c4 =Cs1 =doubleCs2 =symCs3 =symCs4 =sym2 在不加专门指定的情况下，以下符号表达式中的哪⼀个变量被认为是⾃由符号变量.sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)')〖⽬的〗理解⾃由符号变量的确认规则。
〖解答〗symvar(sym('sin(w*t)'),1)ans =wsymvar(sym('a*exp(-X)'),1)ans = asymvar(sym('z*exp(j*th)'),1) ans = z5求符号矩阵=333231232221131211a a a a a a a a a A 的⾏列式值和逆，所得结果应采⽤“⼦表达式置换”简洁化。
〖⽬的〗理解subexpr 指令。
〖解答〗A=sym('[a11 a12 a13;a21 a22 a23;a31 a32 a33]')DA=det(A) IA=inv(A);[IAs,d]=subexpr(IA,d) A =[ a11, a12, a13] [ a21, a22, a23] [ a31, a32, a33] DA =a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31 IAs =[ d*(a22*a33 - a23*a32), -d*(a12*a33 - a13*a32), d*(a12*a23 - a13*a22)] [ -d*(a21*a33 - a23*a31), d*(a11*a33 - a13*a31), -d* (a11*a23 - a13*a21)] [ d*(a21*a32 - a22*a31), -d*(a11*a32 - a12*a31), d*(a11*a22 - a12*a21)] d =1/(a11*a22*a33 - a11*a23*a32 - a12*a21*a33 + a12*a23*a31 + a13*a21*a32 - a13*a22*a31)8（1）通过符号计算求t t y sin )(=的导数dtdy。
- 1、下载文档前请自行甄别文档内容的完整性，平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益，请联系客服反馈,我们会尽快为您处理(人工客服工作时间：9:00-18:30)。
第二单元M2-1（1）>> t=-5:0.01:5; x=(t>0)-(t>2); plot(t,x);axis([-5,5,-2,2]);-5-4-3-2-1012345-2-1.5-1-0.500.511.52（2）t=-5:0.01:5; >> x=(t>0); >> plot(t,x);>> axis([-5,5,-2,2]);-5-4-3-2-1012345-2-1.5-1-0.500.511.52(3)t=-1:0.01:10;>> x=10*exp(-t)-5*exp(-2*t); >> plot(t,x);-20246810-10-55（4）t=-1:0.01:10; >> x=t.*(t>0); >> plot(t,x);-20246810012345678910(5)t=-10:0.001:10;>> pi=3.14;x=2*abs(sin(10*pi*t+pi/3)); >> plot(t,x);-10-8-6-4-2024681000.20.40.60.811.21.41.61.82(6)t=-10:0.001:10;>> x=cos(t)+sin(2*pi*t); >> plot(t,x);-10-8-6-4-20246810-2-1.5-1-0.500.511.52(7)>> t=-10:0.001:10; pi=3.14;>> x=4.*exp(-0.5.*t).*cos(2.*pi*t); plot(t,x);-10-8-6-4-20246810-600-400-200200400600(8)t=-10:0.001:10; >> pi=3.14;x=sin(pi*t)./(pi*t).*cos(30*t); >> plot(t,x);-10-8-6-4-20246810-1-0.8-0.6-0.4-0.200.20.40.60.81M2-2t=-2:0.001:2;>> x=(t>-1)-(t>0)+2*tripuls(t-0.5,1,0); >> plot(t,x)>> axis([-2,2,-2,2])-2-1.5-1-0.500.51 1.52-2-1.5-1-0.500.511.52t=-2:0.001:2; >> pi=3.14;>> x=(t>-1)-(t>0)+2*tripuls(t-0.5,1,0); >> y=x.*cos(50*pi*t); >> plot(t,y)>> axis([-2,2,-2,2])-2-1.5-1-0.500.51 1.52-2-1.5-1-0.500.511.52M2-3function yt=x2_1(t)yt=t.*(t>0)-t.*(t>=2)+2*(t>=2)-3*(t>3)+(t>5);t=-10:0.01:11; >> subplot(3,1,1) >> plot(t,x2_1(t)) >> title('x(t)')>> axis([-1,6,-2,3]) >> subplot(3,1,2) >> plot(t,x2_1(0.5*t)) >> axis([-1,11,-2,3]) >> title('x(0.5t)') >> subplot(3,1,3)>> plot(t,x2_1(2-0.5*t)) >> title('x(2-0.5t)') >> axis([-9,5,-2,3])-1123456-202x(t)246810-22x(0.5t)-8-6-4-224-22x(2-0.5t)四 偶分量 t=-2:0.01:2;>> x=0.5*t.*(t>-1)-t.*(t>0)+0.5*t.*(t>1)+0.5*(t>-1)-0.5*(t>1); >> plot(t,x);>> axis([-2,2,-1,1])-2-1.5-1-0.500.51 1.52-1-0.8-0.6-0.4-0.200.20.40.60.81奇分量 t=-2:0.01:2;>> x=0.5*t.*(t>-1)-(t>0)-0.5*t.*(t>1)+0.5*(t>-1)+0.5*(t>1); >> plot(t,x);>> axis([-2,2,-1,1])-2-1.5-1-0.500.51 1.52-1-0.8-0.6-0.4-0.200.20.40.60.81五（1） x[k]=delta[k] k=-50:50;>> delta=[zeros(1,50),1,zeros(1,50)]; >> stem(k,delta)（2）x[k]=u[k]>> k=-50:50;uk=[zeros(1,50),ones(1,51)]; stem(k,uk)(3) k=-50:50;>> xk=10*(0.5.^k).*[zeros(1,50),ones(1,51)]; >> stem(k,xk)(4) k=-50:50;>> xk=[zeros(1,48),ones(1,7),zeros(1,46)]; >> stem(k,xk)(5) k=-50:50;uk=k.*[zeros(1,50),ones(1,51)]; stem(k,uk)(6)k=-50:50;>> xk=5.*(0.8.^k).*cos(0.9*pi.*k).*[ones(1,101)];>> stem(k,xk)x 105-50-40-30-20-1001020304050M2-6W=0.1pik=-50:50;>> w=0.1*pi;>> xk=sin(w.*k);>> stem(k,xk);W=0.5pik=-50:50;>> w=0.5*pi; >> xk=sin(w.*k); >> stem(k,xk);W=0.9pik=-50:50;>> w=0.9*pi; >> xk=sin(w.*k); >> stem(k,xk);W=1.1pik=-50:50;>> w=1.1*pi; >> xk=sin(w.*k); >> stem(k,xk);W=1.5pik=-50:50;>> w=1.5*pi; >> xk=sin(w.*k); >> stem(k,xk);W=1.9pik=-50:50;>> w=1.9*pi; >> xk=sin(w.*k); >> stem(k,xk);M2-7>> t=-1*pi:pi/100:3*pi; >> subplot(3,2,1)>> x1=cos(6*pi*t); >> plot(t,x1)>> subplot(3,2,2)>> x2=cos(14*pi*t); >> plot(t,x2)>> subplot(3,2,3)>> x3=cos(26*pi*t); >> plot(t,x3)>> subplot(3,2,4)>> k=0:30;>> x1k=cos(3*pi/5*k); >> stem(k,x1k)>> subplot(3,2,5)>> k=0:30;>> x2k=cos(7*pi/5*k); >> stem(k,x2k)>> subplot(3,2,6)>> k=0:30;>> x3k=cos(13*pi/5*k); >> stem(k,x3k)-50510-101-50510-101-50510-101M2-8 k=0:40; >> pi=3.14; >> w=2*pi/10;>> ft=-square(w.*(k+1),50); >> stem(k,ft)k=0:40;>> pi=3.14;>> w=2*pi/20;>> xk=sawtooth(w.*k,0.5); >> stem(k,xk)M2-9>> k=-4:7;x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; >> stem(k,x)>> l=-4:7;>> k=-12:21;>> x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; >> N=length(x);>> y=zeros(1,3*N-2); >> y(1:3:end)=x; >> stem(k,y);-4-3-2-101234>> x=[0 0 -3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; >> x1=x(1:3:end); >> k=-1:3;>> stem(k-1,x1)>> k=-6:5;>> x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; >> stem(k,x)>> k=-2:9;>> x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1];>> stem(k,x)>> k=-4:7;x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; >> xk=fliplr(x);>> k1=-fliplr(k);>>stem(k1,xk)-4-3-2-101234第三单元1）>> ts=0;te=5;dt=0.01; >> sys=tf([2 1],[1 3 2]); >> t=ts:dt:te;>> x=exp(-3*t).*(t>=0); >> y=lsim(sys,x,t); >> plot(t,y);>> xlabel('Time(sec)') >> ylabel('y(t)')012345-0.0500.050.10.150.20.250.3Time(sec)y (t )>> ts=0;te=5;dt=0.0001; >> sys=tf([2 1],[1 3 2]); >> t=ts:dt:te;>> x=exp(-3*t).*(t>=0); >> y=lsim(sys,x,t); >> plot(t,y);>> xlabel('Time(sec)') >> ylabel('y(t)')012345-0.0500.050.10.150.20.250.3Time(sec)y (t )4）x=[0.85,0.53,0.21,0.67,0.84,0.12]; k1=-2:3;h=[0.68,0.37,0.83,0.52,0.71]; k2=-1:3; y=conv(x,h);k=(k1(1)+k2(1)):(k1(end)+k2(end)); stem(k,y);-3-2-101234568）>> k=0:30;>> a=[1 0.7 -0.45 -0.6];>> b=[0.8 -0.44 0.36 0.02];>> h=impz(b,a,k);>> stem(k,h)第五单元（1）>> t=0:0.05:2.5;>> x1t=(t>=0)-(t>=1); >> x2t=tripuls(t-1,2,0); >> xt=x1t+x2t.*cos(50*t); >> plot(t,xt)00.51 1.52 2.5-1-0.50.511.522)>> w=linspace(0,2*pi,200); >> b=[10^4];>> a=[1 26.131 3.4142*10^2 2.6131*10^3 10^4]; >> H=freqs(b,a,w); >> subplot(2,1,1); >> plot(w,abs(H));>> set(gca,'xtick',[0 1 2 3 4 5]);>> set(gca,'ytick',[0 0.4 0.707 1]);grid; >> xlabel('\omega');>> ylabel('|H(j\omega)|') >> subplot(2,1,2); >> plot(w,angle(H));>> set(gca,'xtick',[0 1 2 3 4 5]);grid; >> xlabel('\omega');>> ylabel('\phi(\omega)');0123451ω|H (j ω)|12345-2-1.5-1-0.50ωφ(ω)3) >> ts=0;te=5;dt=0.01;>> sys=tf([10^4],[1 26.131 3.4142*10^2 2.6131*10^3 10^4]); >> t=ts:dt:te; >> t=0:0.05:2.5; x1t=(t>=0)-(t>=1); x2t=tripuls(t-1,2,0); x=x1t+x2t.*cos(50*t); >> y=lsim(sys,x,t); >> plot(t,y);00.51 1.52 2.5-0.20.20.40.60.811.2>> ts=0;te=5;dt=0.01;>> sys=tf([10^4],[1 26.131 3.4142*10^2 2.6131*10^3 10^4]); >> t=ts:dt:te; >> t=0:0.05:2.5;>> x1t=(t>=0)-(t>=1); >> x2t=tripuls(t-1,2,0); >> x3t=cos(50*t); >> x=(x1t+x2t).*x3t; >> y=lsim(sys,x,t); >> plot(t,y);0.511.522.5-0.1-0.08-0.06-0.04-0.0200.020.040.060.080.1结论：输入信号乘以cos 函数，幅度减半，左右各移wo 单位。