(1)
t=-2:0.001:4;
T=2;
xt=rectpuls(t-1,T); plot(t,xt)
axis([-2,4,-0.5,1.5]) 图象为：
(2)
t=sym('t');
y=Heaviside(t); ezplot(y,[-1,1]);
grid on
axis([-1 1 -0.1 1.1]) 图象为：
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])
图象为：
A=2;w0=10*pi;phi=pi/6;
t=0:0.001:0.5;
xt=abs(A*sin(w0*t+phi)); plot(t,xt)
图象为：
(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)
图象为：
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=cos(w0*t).*sinc(t/pi); plot(t,xt)
axis([-15,15,-1.1,1.1]) 图象为：
（1）function yt=x2_3(t)
yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5); （2） function yt=x2_3(t)
yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5);
t=0:0.001:6;
subplot(3,1,1)
plot(t,x2_3(t))
title('x(t)')
axis([0,6,-2,3])
subplot(3,1,2)
plot(t,x2_3(0.5*t))
title('x(0.5t)')
axis([0,11,-2,3])
subplot(3,1,3)
plot(t,x2_3(2-0.5*t))
title('x(2-0.5t)')
axis([-6,5,-2,3])
图像为：
M2-9
（1）
k=-4:7;
xk=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1];
stem(k,xk,'file')
（2）
k=-12:21;
x=[-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,'file')
title('3倍插')
t=-1:2;
y=[-2,-2,2,1];
subplot(2,1,2)
stem(t,y,'file')
title('3倍抽取')
axis([-3,4,-4,4])
k=-4:7;
x=[-3,-2,3,1,-2,-3,-4,2,-1,4,1,-1]; subplot(2,1,1)
stem(k+2,x,'file')
title('x[k+2]')
subplot(2,1,2)
stem(k-4,x,'file')
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),'file') title('x[-k]')
（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)
y =
0.6649
-0.0239
-0.0630
-0.0314
-0.0127
从（1）（2）对比我们当抽样间隔越小时数值精度越高。
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,'file');
M3-8
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,'file')
M4-1
（1）对于周期矩形信号的傅里叶级数cn =-1/2j*sin(n/2*pi)*sinc(n/2) n=-15:15;
X=-j*1/2*sin(n/2*pi).*sinc(n/2);
subplot(2,1,1);
stem(n,abs(X),'file');
title('幅度谱')
xlabel('nw');
subplot(2,1,2);
stem(n,angle(X),'file');
title('相位谱')
（2）对于三角波信号的频谱是：Cn=-++
n=-15:15;
X=sinc(n)-0.5*((sinc(n/2)).^2);
subplot(2,1,1);
stem(n,abs(X),'file');
title('幅度谱')
xlabel('nw');
subplot(2,1,2);
stem(n,angle(X),'file');
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),'file'); title('X[m]实部') subplot(2,1,2);
stem(m,imag(X),'file'); title('X[m]虚部')
M4-7（3）
k=0:10;
x=0.5.^k;
subplot(3,1,1);
stem(k,x,'file')
title('x[k]')
X=fft(x,10);
subplot(3,1,2);
m=0:9;
stem(m,real(X),'file');
title('X[m]实部')
subplot(3,1,3);
stem(m,imag(X),'file');
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,1000 0];
w=linspace(0,2*pi,200);
[H,w]=freqs(b,a,w);
subplot(2,1,1)
plot(w,abs(H));
title('幅度曲线')
subplot(2,1,2)
plot(w,angle(H));
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
（1）
num=[41.6667];
den=[1 3.7444 25.7604 41.6667];
[r,p,k]=residue(num,den)
r =-0.9361 - 0.1895i -0.9361 + 0.1895i 1.8722
p =-0.9361 + 4.6237i -0.9361 - 4.6237i -1.8722
k =[]
（2）
num=[16 0 0];
den=[1 5.6569 816 2262.7 160000];
[r,p,k]=residue(num,den)
r =0.0992 - 1.5147i 0.0992 + 1.5147i -0.0992 + 1.3137i -0.0992 - 1.3137i p =-1.5145 +21.4145i -1.5145 -21.4145i -1.3140 +18.5860i -1.3140 -18.5860i k =[]
（3）
num=[1 0 0 0];
den=conv([1 5],[1 5 25]);
[r,p,k]=residue(num,den)
r =-2.5000 - 1.4434i -2.5000 + 1.4434i -5.0000
p =-2.5000 + 4.3301i -2.5000 - 4.3301i -5.0000
k =1
（4）
num=[833.3025];
den=conv([1 4.1123 28.867],[1 9.9279 28.867]);
[r,p,k]=residue(num,den)
r =-2.4819 + 1.0281i -2.4819 - 1.0281i 2.4819 - 5.9928i 2.4819 + 5.9928i p =-2.0562 + 4.9638i -2.0562 - 4.9638i -4.9640 + 2.0558i -4.9640 - 2.0558i k =[]
M6-2
sys=tf(b,a);
x=1*(t>0)+0*(t<=0);
y1=lsim(sys,x,t);
plot(t,y1);
title('零状态响应')
figure;
[A B C D]=tf2ss(b,a);
sys=ss(A,B,C,D);
x=1*(t>0)+0*(t<=0);z
i=[1 2];
y=lsim(sys,x,t,zi);
plot(t,y);
title('完全响应')
figure;
y2=y-y1;
plot(t,y2);
title('零输入响应')
M6-5
a=[1 2 2 1]; b=[1 2];
sys=tf(b,a) ;
pzmap(sys)
冲击响应
a=[1 2 2 1]; b=[1 2];
sys=tf(b,a) ;
impulse(sys )
阶跃响应
a=[1 2 2 1]; b=[1 2];
sys=tf(b,a); step(sys)
频率响应
a=[1 2 2 1];
b=[1 2];
w=linspace(-10,10
,20000);
H=freqs(b,a,w);
plot(w,abs(H))
M7-1
（1）
num=[2 16 44 56 32];
den=[3 3 -15 18 -12];
[r,p,k]=residuez(num,den)
r =-0.0177 9.4914 -3.0702 + 2.3398i -3.0702 - 2.3398i
p =-3.2361 1.2361 0.5000 + 0.8660i 0.5000 - 0.8660i
k =-2.6667
（2）
num=[4 -8.86 -17.98 26.74 -8.04];
den=[1 -2 10 6 65];
[r,p,k]=residuez(num,den)
r =1.0849 + 1.3745i 1.0849 - 1.3745i 0.9769 - 1.2503i 0.9769 + 1.2503i
p =2.0000 + 3.0000i 2.0000 - 3.0000i -1.0000 + 2.0000i -1.0000 - 2.0000i k =-0.1237
M7-2
k=0:10;
a=[2 -1 3];
b=[2 -1];
y1=filter(b,a,x);
stem(k,y1)
title('零状态响应')
figure;
h=impz(b,a,k);
x=0.5.^k;
y2=conv(x,h);
N=length(y2);
stem(0:N-1,y2)
title('零输入响应')
完全响应
y=y1+y2
M7-3
(1)
num=[2 16 44 56 32];
den=[3 3 -15 18 -12];
zplane(num,den)
由图可知极点不是都在单位圆，所以系统不是稳定系统。
（2）
num=[4 -8.68 -17.98 26.74 -8.04];
den=[1 -2 10 6 65];
zplane(num,den)
由图可知极点不是都在单位圆，所以系统不是稳定系统。
习题2 1． 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度”对象，还是“符号” 对象？ 3/7+0.1, sym(3/7+0.1), vpa(sym(3/7+0.1)) a=class(3/7+0.1)%双精度 b=class(sym(3/7+0.1))%符号 c=class(vpa(sym(3/7+0.1),4))%符号 d=class(vpa(sym(3/7+0.1)))%符号 2． 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认为是独立自由变量。 sym('sin(w*t)') , sym('a*exp(-X)' ) , sym('z*exp(j*th)') a=sym('sin(w*t)'); symvar(a) b=sym('a*exp(-X)'); symvar(b) c=sym('z*exp(j*th)'); symvar(c) 3． 求以下两个方程的解： （提示：关于符号变量的假设要注意） （1）试写出求三阶方程05.443 =-x 正实根的程序。注意：只要正实根，不要出现其他根。 x=sym('x','positive'); f=x^3-44.5; x=solve(f,x) （2）试求二阶方程02 2=+-a ax x 在0>a 时的根。 a=sym('a','positive'); syms x; f=x^2-a*x+a^a; x=solve(f,x) 4． 观察一个数（在此用@记述）在以下四条不同指令作用下的异同： a = @ , b = sym( @ ), c = sym( @ ,' d ' ), d = sym( '@ ' ) 在此，@ 分别代表具体数值 7/3 , pi/3 , pi*3^(1/3) ；而异同通过vpa(abs(a-d)) , vpa(abs(b-d)) , vpa(abs(c-d))等来观察。 a=7/3 b=sym(7/3) c=sym(7/3,'d') d=sym('7/3') vpa(abs(a-d)) vpa(abs(b-d)) vpa(abs(c-d)) a=pi/3 b=sym(pi/3) c=sym(pi/3,'d')
第1章 MATLAB概论 1.1与其他计算机语言相比较，MATLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MATLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MATLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 在安装MATLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MATLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MATLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 1.8 如何设置当前目录和搜索路径，在当前目录上的文件和在搜索路径上的文件有什么区别？ 当前目录可以在当前目录浏览器窗口左上方的输入栏中设置，搜索路径可以通过选择操作桌面的file 菜单中的Set Path菜单项来完成。在没有特别说明的情况下，只有当前目录和搜索路径上的函数和文件能够被MATLAB运行和调用，如果在当前目录上有与搜索路径上相同文件名的文件时则优先执行当前目录上的文件，如果没有特别说明，数据文件将存储在当前目录上。 1.9 在MATLAB中有几种获得帮助的途径？
《信号与系统MATLAB实践》第一次上机作业 实验一、熟悉MATLAB基本操作 三、基本序列运算 1.数组的加减乘除和乘方运算 A=[1 2 3]; B=[4 5 6]; C=A+B; D=A-B; E=A.*B; F=A./B; G=A.^B; subplot(2,4,1);stem(A) subplot(2,4,2);stem(B) subplot(2,4,3);stem(C) subplot(2,4,4);stem(D) subplot(2,4,5);stem(E) subplot(2,4,6);stem(F) subplot(2,4,7);stem(G) 2.绘制函数波形 （1）t=0:0.001:10
x=3-exp(-t); plot(t,x) ylabel('f(t)'); xlabel('t'); title('(1)'); （2）t=0:0.001:10 x=5*exp(-t)+3*exp(-2*t); plot(t,x) ylabel('f(t)'); xlabel('t'); title('(2)');
（3）t=0:0.001:3 x=exp(-t).*sin(2*pi*t); plot(t,x) ylabel('f(t)'); xlabel('t'); title('(3)'); （4）t=0:0.001:3 x=sin(3*t)./(3*t);
plot(t,x) ylabel('f(t)'); xlabel('t'); title('(4)'); （5）k=1:1:6 x=(-2).^(-k); stem(k) xlabel('k'); ylabel('f(k)'); title('(5)');
第二章 1 求下列表达式的值。 （1） w=sqrt(2)*(1+0.34245*10^(-6)) (2) a=3.5; b=5; c=-9.8; x=(2*pi*a+(b+c)/(pi+a*b*c)-exp(2))/tan(b+c)+a (3) a=3.32; b=-7.9; y=2*pi*a^(2)*[(1-pi/4)*b-(0.8333-pi/4)*a] (4) t=[2,1-3*i;5,-0.65]; z=1/2*exp(2*t)*log(t+sqrt(1+t^(2))) 2 求下列表达式 A=[-1,5,-4;0,7,8;3,61,7]; B=[8,3,-1;2,5,3;-3,2,0]; (1) A+6*B A^2-B+eye (2) A*B A.*B B.*A (3) A/B B\A (4) [A,B] [A([1,3],:);B^2] 3 根据已知，完成下列操作 （1） A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14]; K=find(A>10&A<25); A(K) （2） A=[23,10,-0.778,0;41,-45,65,5;32,5,0,32;6,-9.54,54,3.14]; B=A(1:3,:) C=A(:,1:2) D=A(2:4,3:4) E=B*C （3） E 2 用if语句 score=input('请输入成绩:'); if score>=90&&score<=100 disp('A'); elseif score>=80&&score<=89 disp('B'); elseif score>=70&&score<=79 disp('C'); elseif score>=60&&score<=69; disp('D'); elseif score<60&&score>=0; disp('E'); else disp('出错'); end 用switch语句 score=input('请输入成绩:'); switch fix(score/10) case {9,10} disp('A'); case {8} disp('B'); case {7} disp('C'); case {6} disp('D'); case {0,1,2,3,4,5} disp('E'); otherwise disp('出错'); end 第四章1题 1) X=0:10; Y=x-x.^3/6; P lot(x,y) 2) t=0:0.01:2*pi; x=8.*cos(t); y=4*sqrt(2).*sin(t); plot(x,y) 2题 第2章符号运算 习题2及解答 1 说出以下四条指令产生的结果各属于哪种数据类型，是“双精度” 对象，还是“符号”符号对象？ 3/7+0.1; sym(3/7+0.1); sym('3/7+0.1'); vpa(sym(3/7+0.1)) 〖目的〗 ●不能从显示形式判断数据类型，而必须依靠class指令。 〖解答〗 c1=3/7+0.1 c2=sym(3/7+0.1) c3=sym('3/7+0.1') c4=vpa(sym(3/7+0.1)) Cs1=class(c1) Cs2=class(c2) Cs3=class(c3) Cs4=class(c4) c1 = 0.5286 c2 = 37/70 c3 = 0.52857142857142857142857142857143 c4 = 0.52857142857142857142857142857143 Cs1 = double Cs2 = sym Cs3 = sym Cs4 = sym 2 在不加专门指定的情况下，以下符号表达式中的哪一个变量被认 为是自由符号变量. sym('sin(w*t)'),sym('a*exp(-X)'),sym('z*exp(j*th)') 〖目的〗 ●理解自由符号变量的确认规则。 〖解答〗 symvar(sym('sin(w*t)'),1) ans = w symvar(sym('a*exp(-X)'),1) ans = a symvar(sym('z*exp(j*th)'),1) ans = z 3 求以下两个方程的解 （1）试写出求三阶方程05.443 =-x 正实根的程序。注意：只要正实根，不要出现其他根。 （2）试求二阶方程022=+-a ax x 在0>a 时的根。 〖目的〗 ● 体验变量限定假设的影响 〖解答〗 （1）求三阶方程05.443 =-x 正实根 reset(symengine) %确保下面操作不受前面指令运作的影响 syms x positive solve(x^3-44.5) ans = (2^(2/3)*89^(1/3))/2 （2）求五阶方程02 2 =+-a ax x 的实根 syms a positive %注意：关于x 的假设没有去除 solve(x^2-a*x+a^2) Warning: Explicit solution could not be found. > In solve at 83 ans = [ empty sym ] syms x clear 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 4 观察一个数（在此用@记述）在以下四条不同指令作用下的异同。 a =@, b = sym( @ ), c = sym( @ ,' d ' ), d = sym( '@ ' ) 在此，@ 分别代表具体数值 7/3 , pi/3 , pi*3^(1/3) ；而异同通过vpa(abs(a-d)) , vpa(abs(b-d)) , vpa(abs(c-d))等来观察。 〖目的〗 ● 理解准确符号数值的创建法。 ● 高精度误差的观察。 〖解答〗 （1）x=7/3 x=7/3;a=x,b=sym(x),c=sym(x,'d'),d=sym('7/3'), a = P16 Q2: 计算表达式()2 tan arccos x x -在0.25x =和0.78x π =时的函数值。 function y=jie(x) y=tan(-x.^2)*acos(x); >> jie(0.25) ans = -0.0825 >> jie(0.78*pi) ans = 0 + 0.4418i Q3：编写M 命令文件，求5010 2 1 1 1k k k k ==+ ∑∑ 的值。 a=0;b=0; for i=1:50 a=a+i*i; end for j=1:10 b=b+1/j; end c=a+b; >> c c = 4.2928e+004 P27 Q2：矩阵1234567 8 9A ????=??????,4 685563 2 2B ?? ?? =? ????? ,计算A B *,.A B *,并比较两者的区别。 >> A=[1 2 3;4 5 6;7 8 9]; >> B=[4 6 8;5 5 6;3 2 2]; >> A*B ans = 23 22 26 59 61 74 95 100 122 >> A.*B ans = 4 12 24 20 25 36 21 16 18 A*B表示A与B两矩阵相乘。 A.*B表示A与B对应元素相乘。 P34 Q2：编写一个转换成绩等级的程序，其中成绩等级转换标准为：考试分数在[] 90,100显示为优秀；分数在[) 0,60的 60,80的显示为及格；分数在[) 80,90的显示为良好；分数在[) 显示为不及格。 if x>=90 disp('优秀'); elseif x>=80 disp('良好'); elseif x>=60 disp('及格'); else disp('不及格'); end >> x=85 x = 85 良好 第一章 M A T L A B 概况与基本操作 1.选择题（每题2分，共20分）： (1)最初的MATLAB 核心程序是采用D 语言编写的。 A.PASCAL B.C C.BASIC D.FORTRAN (2)即将于2011年9月发布的MATLAB 新版本的编号为C 。 A.MATLAB 2011Ra B.MATLAB 2011Rb C.MATLAB R2011b D.MATLAB R2011a (3)在默认设置中，MATLAB 中的注释语句显示的颜色是B 。 A.黑色 B.绿色 C.红色 D.蓝色 (4)如果要以科学计数法显示15位有效数字，使用的命令是A 。 A.format long e B.format long C.format long g D.format long d (5)在命令窗口新建变量a 、b ，如果只查看变量a 的详细信息，使用的命令为A 。 A.whos a B.who a C.who D.whos (6)如果要清除工作空间的所有变量，使用的命令为C 。 A.clear B.clear all C.两者都可 D.两者都不可 (7)在创建变量时，如果不想立即在命令窗口中输出结果，可以在命令后加上B 。 A.冒号 B.分号 C.空格 D.逗号 (8)如果要重新执行以前输入的命令，可以使用D 键。 A.下箭头↓ B.右箭头→ C.左箭头← D.上箭头↑ (9)如果要查询函数det 的功能和用法，并显示在命令窗口，应使用命令C 。 A.doc B.lookfor C.help D.三者均可 (10)如果要启动Notebook 文档，下列D 操作是可行的。 A.在命令窗口输入notebook 命令 B.在命令窗口输入notebook filename 命令 C.在Word 中启动M-book 文档 D.三者均可 2.填空题（每空1分，共20分）： (1)MATLAB 是matrix 和laboratory 两个单词前三个字母的组合，意为“矩阵实验室”，它的创始人是Cleve Moler 和Jack Little 。 (2)在MATLAB 的默认设置中，关键字显示的字体为蓝色，命令、表达式、计算结果显示的字体为黑色，字符串显示的字体为褐红色，注释显示的字体为绿色，错误信息显示的字体为红色。 (3)在命令窗口中，输出结果显示为各行之间添加空行的命令为format loose ，各行之间不添加空行的命令为format compact 。 (4)在MATLAB 中，各种标点符号的作用是不同的。例如，空格的作用是分隔数组每行各个元素，逗号的作用是分隔数组每行各个元素或函数的各个输入参数，分号的作用是作为不显示命令结果的命令行的结尾或分隔数组各列，冒号的作用是生成一维数组或表示数组全部元素，百分号的作用是引导一行注释，…的作用是连接相邻两行，感叹号的作用是调用操作系统命令。 3.程序设计题（每题10分，共40分） (1)以25m/s 的初速度向正上方投球（g=9.8m/s 2），计算到达最高点的时间tp 以及球从出发点到 最高点的距离hp 。 解：根据物理学知识，物体上抛运动的速度与经过的时间之间的关系为0p p v v gt =-，因此所需要的时间为0p p v v t g -=。而到达最高点时的速度0p v =，因此可根据此公式求出tp ： v0=25;g=9.8;vp=0; tp=(v0-vp)/g tp = 2.5510 M2-3 (1) function yt=x(t) yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5); (2)function yt=x (t) yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5); t=0:0.001:6; subplot(3,1,1) plot(t,x2_3(t)) title('x(t)') axis([0,6,-2,3]) subplot(3,1,2) plot(t,x2_3(0.5*t)) title('x(0.5t)') axis([0,11,-2,3]) subplot(3,1,3) plot(t,x2_3(2-0.5*t)) title('x(2-0.5t)') axis([-6,5,-2,3]) 图像为： M2-5 (3) function y=un(k) y=(k>=0) untiled3.m k=[-2:10] xk=10*(0.5).^k.*un(k); stem(k,xk) title('x[k]') axis([-3,12,0,11]) M2-5 (6) k=[-10:10] xk=5*(0.8).^k.*cos((0.9)*pi*k) stem(k,xk) title('x[k]') grid on M2-7 A=1; t=-5:0.001:5; w0=6*pi; xt=A*cos(w0*t); plot(t,xt) hold on A=1; k=-5:5; w0=6*pi; xk=A*cos(w0*0.1*k); stem(k,xk) axis([-5.5,5.5,-1.2,1.2]) title('x1=cos(6*pi*t)&x1[k]') 【例1.3-5】图示复数i z i z 21,3421+=+=的和。 z1=4+3*i;z2=1+2i; z12=z1+z2 clf,hold on plot([0,z1,z12],'-b','LineWidth',3) plot([0,z12],'-r','LineWidth',3) plot([z1,z12],'ob','MarkerSize',8) hold off,grid on axis equal axis ([0,6,0,6]) text(3.5,2.3,'z1') text(5,4.5,'z2') text(2.5,3.5,'z12') xlabel('real') ylabel('image') shg z12 = a=-8; r_a=a^(1/3) p=[1,0,0,-a]; R=roots(p) MR=abs(R(1)); t=0:pi/20:2*pi; x=MR*sin(t); y=MR*cos(t); plot(x,y,'b:'),grid on hold on plot(R(2),'.','MarkerSize',30,'Color','r') plot(R([1,3]),'o','MarkerSize',15,'Color','b') axis([-3,3,-3,3]),axis square hold off r_a = 1.0000 + 1.7321i R = -2.0000 1.0000 + 1.7321i 【例1.3-10】画出衰减振荡曲线t e y t 3sin 3-=，t 的取值范围是]4,0[π。 t=0:pi/50:4*pi; y=exp(-t/3).*sin(3*t); plot(t,y,'r','LineWidth',2) axis([0,4*pi,-1,1]) xlabel('t'),ylabel('y') 第4章数值运算 习题 4 及解答 1 根据题给的模拟实际测量数据的一组t和)(t y试用数值差分 diff或数值梯度gradient指令计算)(t y'曲线 y'，然后把)(t y和)(t 绘制在同一图上，观察数值求导的后果。（模拟数据从prob_data401.mat获得） 〖目的〗 ●强调：要非常慎用数值导数计算。 ●练习mat数据文件中数据的获取。 ●实验数据求导的后果 ●把两条曲线绘制在同一图上的一种方法。 〖解答〗 （1）从数据文件获得数据的指令 假如prob_data401.mat文件在当前目录或搜索路径上 clear load prob_data401.mat （2）用diff求导的指令 dt=t(2)-t(1); yc=diff(y)/dt; %注意yc的长度将比y短1 plot(t,y,'b',t(2:end),yc,'r') （3）用gradent求导的指令（图形与上相似） dt=t(2)-t(1); yc=gradient(y)/dt; plot(t,y,'b',t,yc,'r') grid on 〖说明〗 ● 不到万不得已，不要进行数值求导。 ● 假若一定要计算数值导数，自变量增量dt 要取得比原有数据相对误差高1、2个量级 以上。 ● 求导会使数据中原有的噪声放大。 2 采用数值计算方法，画出dt t t x y x ? =0sin )(在]10 ,0[区间曲线，并计算)5.4(y 。 〖提示〗 ● 指定区间的积分函数可用cumtrapz 指令给出。 ● )5.4(y 在计算要求不太高的地方可用find 指令算得。 〖目的〗 ● 指定区间的积分函数的数值计算法和cumtrapz 指令。 ● find 指令的应用。 〖解答〗 dt=1e-4; t=0:dt:10; t=t+(t==0)*eps; f=sin(t)./t; s=cumtrapz(f)*dt; plot(t,s,'LineWidth',3) ii=find(t==4.5); s45=s(ii) s45 = 《及应用》实验指导书 《及应用》实验指导书 班级： T1243-7 姓名：柏元强 学号： 20120430724 总评成绩： 汽车工程学院 电测与汽车数字应用中心 目录 实验04051001 语言基础..................... 错误!未指定书签。实验04051002 科学计算及绘图............. 1错误!未指定书签。实验04051003 综合实例编程.. (31) 实验04051001 语言基础 1实验目的 1) 熟悉的运行环境 2) 掌握的矩阵和数组的运算 3) 掌握符号表达式的创建 4) 熟悉符号方程的求解 2实验内容 第二章 1. 创建的变量，并进行计算。 (1) 87，190，计算 、、a*b 。 (87); (190); *b (2) 创建 8 类型的变量，数值与(1)中相同，进行相同的计算。 8(87); 8(190); *b 2．计算： (1) 操作成绩 报告成绩 (2) e3 (3) (60) (3) (3*4) 3．设，，计算： (1) (2) (3) 23; (4*u*v)(v) (((u))^2)/(v^2) ((3*v))/(u*v) 4．计算如下表达式： (1) (2) (3-5*i)*(4+2*i) (2-8*i) 5．判断下面语句的运算结果。 (1) 4 < 20 (2) 4 <= 20 (3) 4 20 (4) 4 20 (5) 'b'<'B' 4 < 20 , 4 <= 20,4 20,4 20,'b'<'B' 6．设，，，，判断下面表达式的值。 (1) (2) (3) (4) (5) (6) 395837; a><>>> 7．编写脚本，计算上面第2题中的表达式。 ('(60)='); ((60)) ('(3)='); ((3)) ('(3*4)='); ((3*4)) 8．编写脚本，输出上面第6题中的表达式的值。395837; 第1章MATLAB概论 1.1与其他计算机语言相比较，MA TLAB语言突出的特点是什么？ MATLAB具有功能强大、使用方便、输入简捷、库函数丰富、开放性强等特点。 1.2 MA TLAB系统由那些部分组成？ MATLAB系统主要由开发环境、MA TLAB数学函数库、MATLAB语言、图形功能和应用程序接口五个部分组成。 1.3 安装MATLAB时，在选择组件窗口中哪些部分必须勾选，没有勾选的部分以后如何补安装？ 在安装MA TLAB时，安装内容由选择组件窗口中个复选框是否被勾选来决定，可以根据自己的需要选择安装内容，但基本平台（即MATLAB选项）必须安装。第一次安装没有选择的内容在补安装时只需按照安装的过程进行，只是在选择组件时只勾选要补装的组件或工具箱即可。 1.4 MA TLAB操作桌面有几个窗口？如何使某个窗口脱离桌面成为独立窗口？又如何将脱离出去的窗口重新放置到桌面上？ 在MA TLAB操作桌面上有五个窗口，在每个窗口的右上角有两个小按钮，一个是关闭窗口的Close按钮，一个是可以使窗口成为独立窗口的Undock按钮，点击Undock按钮就可以使该窗口脱离桌面成为独立窗口，在独立窗口的view菜单中选择Dock ……菜单项就可以将独立的窗口重新防止的桌面上。 1.5 如何启动M文件编辑/调试器？ 在操作桌面上选择“建立新文件”或“打开文件”操作时，M文件编辑/调试器将被启动。在命令窗口中键入edit命令时也可以启动M文件编辑/调试器。 1.6 存储在工作空间中的数组能编辑吗？如何操作？ 存储在工作空间的数组可以通过数组编辑器进行编辑：在工作空间浏览器中双击要编辑的数组名打开数组编辑器，再选中要修改的数据单元，输入修改内容即可。 1.7 命令历史窗口除了可以观察前面键入的命令外，还有什么用途？ 命令历史窗口除了用于查询以前键入的命令外，还可以直接执行命令历史窗口中选定的内容、将选定的内容拷贝到剪贴板中、将选定内容直接拷贝到M文件中。 实验一常用连续时间信号的实现 一、实验目的 (1)了解连续时间信号的特点； （2）掌握连续时间信号表示的向量法和符号法； （3）熟悉MA TLABPlot函数等的应用。 二、涉及的MATLAB函数 1.plot函数 功能：在X轴和Y轴方向都按线性比例绘制成二维图形。 2.ezplot函数 功能：绘制符号函数在一定范围内的二维图形，简易绘制函数曲线。 3．Sym函数 功能：定义信号为符号变量。 4.subplot函数 功能：产生多个绘图区间。 三、实验内容与方法 1.正弦交流信号f(t)=sin(ωt+φ) （1）符号推理法生成正弦交流信号。 MATLAB程序：. t=-0:0.001:1; f=sym('sin(2*pi*t)'); ezplot(f,[0,1]); xlabel('时间(t)'); ylabei('幅值（f)'); title(‘正弦交流信号'); 用符号法生成的正弦交流信号如图所示： （2）数值法生成正弦交流信号。 MATLAB程序：. t=-0:0.001:1; y=sin(2*pi*t); plot(t,y,'k'); xlabel('时间(t)'); ylabei('幅值（f)'); title('正弦交流信号')； 用数值法生成的正弦交流信号如图所示： 2.单边衰减指数信号. MATLAB程序： t1=-1;t2=10;dt=0.1; t=t1:dt:t2; A1=1; %斜率 a1=0.5; %斜率 n=A1*exp(-a1*t); plot(t,n); axis([t1,t2,0,1]); xlabel('时间（t)'); ylabel('幅值(f)'); title('单边衰减指数信号'); 用数值法生成的单边衰减指数信号如图所示： 目录 第一章 (1) 第二章 (5) 第三章 (12) 第四章 (32) 第五章 (47) 第六章 (54) 补充题欧拉法，龙格库塔法解方程，黑板上的题 (57) 第一章 1.创建表达式 %可以用syms先符号运算再带入值 x=1; y=2; z=(sqrt(4*x^2+1)+0.5457*exp(-0.75*x^2-3.75*y^2-1.5*x))/(2*sin(3*y)-1) z = -1.4345 2.计算复数 x=(-1+sqrt(-5))/4; y=x+8+10j y = 7.7500 +10.5590i 3.help命令学三维曲线 x=-5:0.1:5; y=x; [X,Y]=meshgrid(x,y); Z=(sin(sqrt(X.^2+Y.^2)))./(sqrt(X.^2+Y.^2)); subplot(221); surf(X,Y,Z); colormap(cool); subplot(222); plot3(X,Y,Z,'linewidth',4); %绘制三维曲线，也可以随意给定一个三维曲线的函数。如果画这个曲面，那么将绘出一族三维曲线 grid on; subplot(223); meshz(X,Y,Z); %地毯绘图 subplot(224); meshc(X,Y,Z); %等高线绘图 4.peaks等高线（更改原函数） subplot(221); contour(peaks1,20); subplot(222); contour3(peaks1,10); %可以定义等高线条数 subplot(223); contourf(peaks1,10); subplot(224); peaks1; z = 3*(1-x).^2.*exp(-(x.^2) - (y+1).^2) ... - 10*(x/5 - x.^3 - y.^5).*exp(-x.^2-y.^2) ... - 1/3*exp(-(x+1).^2 - y.^2) (1) t=-2:0.001:4; T=2; xt=rectpuls(t-1,T); plot(t,xt) axis([-2,4,-0.5,1.5]) 图象为： (2) t=sym('t'); y=Heaviside(t); ezplot(y,[-1,1]); grid on axis([-1 1 -0.1 1.1]) 图象为： 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]) 图象为： A=2;w0=10*pi;phi=pi/6; t=0:0.001:0.5; xt=abs(A*sin(w0*t+phi)); plot(t,xt) 图象为： (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) 图象为： 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=cos(w0*t).*sinc(t/pi); plot(t,xt) axis([-15,15,-1.1,1.1]) 图象为： （1）function yt=x2_3(t) yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5); （2）function yt=x2_3(t) yt=(t).*(t>=0&t<=2)+2*(t>=2&t<=3)-1*(t>=3&t<=5); t=0:0.001:6; subplot(3,1,1) plot(t,x2_3(t)) title('x(t)') axis([0,6,-2,3]) subplot(3,1,2) plot(t,x2_3(0.5*t)) title('x(0.5t)') axis([0,11,-2,3]) subplot(3,1,3) plot(t,x2_3(2-0.5*t)) title('x(2-0.5t)') axis([-6,5,-2,3]) 图像为： 第二章 1·求下列表达式的值。 （1）w=sqrt(2)*(1+*10^-6) w = （2）a=;b=5;c=; x=(2*pi*a+(c+b)/(pi+a*b*c)-exp(2))/(tan(b+c)+a); x x = （3）a=;b=; y=2*pi*a^2*[(1-pi/4)*b-4)*a]; y y = （4）t=[2,1-3i;5,]; z=1/2*exp(2*t)*log(t+sqrt(1+t^2)); z z = +004 * - - - - 2，已知a，b，求下列表达式的值。 a=[-1,5,-4;0,7,8;3,61,7];b=[8,3,-1;2,5,3;-3,2,0];（1）a+6*b ans = 47 23 -10 12 37 26 -15 73 7 a^2-b+eye(3) ans = -18 -217 17 22 533 109 21 867 526（2）a*b ans = 14 14 16 -10 51 21 125 328 180 a.*b ans = -8 15 4 0 35 24 -9 122 0 b*a ans = -11 0 -15 7 228 53 3 -1 28（3）a/b ans = b\a ans = （4）[a,b] ans = -1 5 -4 8 3 -1 0 7 8 2 5 3 3 61 7 -3 2 0 [a([1,3],:);b^2] ans = -1 5 -4 3 61 7 73 37 1 17 37 13 -20 1 9 3.已知a，完成下列操作。 a=[23,10,,0;41,-45,65,5;32,5,0,32;6,,54,]; （1）输出a在[10,25]范围内的全部元素。 k=find(a>10&a<25) a(k) k = 1 ans = 23 (2)取出a前3行构成矩阵b，前两列构成矩阵c，右下角3*2子矩阵构成矩阵d，b与c的乘积构成矩阵e。 b=a(1:3,:) 汽车理论作业MA TLAB过程 汽车驱动力与阻力平衡图 加速度倒数-速度曲线图 0 10 20 30 40 50 60 70 80 90 100 u 汽车功率平衡图 u/(km/h)最高档等速百公里油耗曲线 Ua/(km/h) 燃油积极性-加速时间曲线 源程序： 《第一章》 m=3880; g=9.8; r=0.367; x=0.85; f=0.013; io=5.83; CdA=2.77; lf=0.218; Iw1=1.798; Iw2=3.598; Iw=lw1+lw2; ig=[6.09 3.09 1.71 1.00]; %变速器传动比 L=3.2; a=1.947; hg=0.9; n=600:1:4000; T=-19.313+295.27* n/1000-165.44*(门/1000)人2+40.874*(门/1000)人3-3.8445*( n/IOOO).%; Ft1=T*ig(1)*io*x/r; %计算各档对应转速下的驱动力 Ft2=T*ig(2)*io*x/r; Ft3=T*ig(3)*io*x/r; Ft4=T*ig(4)*io*x/r; u1=0.377*r*n/(io*ig(1)); u2=0.377*r*n/(io*ig(2)); u3=0.377*r*n/(io*ig(3)); u4=0.377*r*n/(io*ig(4)); u=0:130/3400:130; F仁m*g*f+CdA*u”2/21.15;%计算各档对应转速下的驱动阻力 F2=m*g*f+CdA*u2.A2/21.15; F3=m*g*f+CdA*u3.A2/21.15; F4=m*g*f+CdA*u4.A2/21.15; figure(1); plot(u1,Ft1, '-r' ,u2,Ft2, '-m' ,u3,Ft3, '-k' ,u4,Ft4, '-b' ,u1,F1, '-r' ,u2,F2, '-m' ,u3,F3, ' k' ,u4,F4, '-b' , 'LineWidth' ,2) title( ' 汽车驱动力与阻力平衡图' ); xlabel( 'u_{a}/km.hA{-1}' ) ylabel( 'F/N' ) gtext( 'F_{t1}' ) gtext( 'F_{t2}' ) gtext( 'F_{t3}' ) gtext( 'F_{t4}' ) gtext( 'F_{f}+F_{w}' ) %由汽车驱动力与阻力平衡图知，他们无交点，u4在最大转速时达到最大 umax=u4(3401) Ft1max=max(Ft1); imax=(Ft1max-m*g*f)/(m*g) disp( ' 假设是后轮驱动' ); C=imax/(a/L+hg*imax/L) % 附着率 delta1=1+(Iw1+Iw2)/(m*rA2)+If*ig(1)*rA2*ioA2*x/(m*rA2); delta2=1+(Iw1+Iw2)/(m*rA2)+If*ig(2)*rA2*ioA2*x/(m*rA2); delta3=1+(Iw1+Iw2)/(m*rA2)+If*ig(3)*rA2*ioA2*x/(m*rA2); delta4=1+(Iw1+Iw2)/(m*rA2)+If*ig(4)*rA2*ioA2*x/(m*rA2); a1=(Ft1-F1)/(delta1*m); %加速度 a2=(Ft2-F2)/(delta2*m); a3=(Ft3-F3)/(delta3*m); a4=(Ft4-F4)/(delta4*m); h1=1./a1; %加速度倒数 h2=1./a2; h3=1./a3; h4=1./a4; figure(2); 考核人数______ 考核班次_______________ 任课教员_________ 出题教员签名________ 任课教研室主任签名_______日期_______ 队别__________ 教学班次___________ 学号___________ 姓名____________ …………………………密………………………………封………………………………线……………………………………… 通信系统仿真题目 1．学习电路时已知LC 谐振电路具有选择频率的作用，当输入正弦信号频率与LC 电路的谐 振频率一致时，将产生较强的输出响应，而当输入信号频率适当偏离时，输出响应相对值很弱，几乎为零（相当于窄带通滤波器）。利用这一原理可以从非正弦周期信号中选择所需的正弦频率成分。题图所示RLC 并联电路和电流1()i t 都是理想模型。已知电路的谐振频率为 0100f kHz = =,100R k =Ω谐振电路品质因素Q 足够高（可滤除邻近频率成分） 1()i t 为周期矩形波，幅度为1 mA 当1()i t 的参数(,)T τ为下列情况时，粗略地画出输出电压 2()t υ的波形，并注明幅度值。 （1）510s T s τμμ== （2）1020s T s τμμ== （3）1530s T s τμμ== 2．设()x n 为一限长序列，当0n <和n N ≥时，()0x n =，且N 等于偶数。已知[()]DFT x n = ()X k ，试用()X k 表示以下各序列的DFT 。 (1)1()(1)x n x N n =-- (2)2()(1)()n x n x n =- (3) 3() (01)()()(21)0()x n n N x n x n N N n N n ≤≤-?? =-≤≤-??? 为其他值 (4) 4()()(01) ()2 2 () N N x n x n n x n n ?≠+≤≤ -? =???为其他值 (5) 5()(01)()0 (21)0() x n n N x n N n N n ≤≤-?? =≤≤-??? 为其他值 (6) 6() ()20()n x n x n n ??? ? ?=????? 为偶数为奇数 (DFT 有限长度取2N ，k 取偶数。) (7) 7()(2)x n x n =(DFT 有限长度取 2 N )。 3．已知三角脉冲1() f t 的傅里叶变换为21()24E F Sa τωτω??= ??? 试利用有关定理求210()cos()2f t f t t τω?? =- ??? 的傅里叶变换2()F ω。1()f t 、2()f t 的波形如下图所示。 4．求下图所示半波余弦信号的傅里叶级数。若E=10V ，f=10kHz ，大致画出幅度谱。 5．求下图所示()F ω的傅里叶逆变换()f t 。matlab课后习题解答第二章doc
Matlab课后习题解答
MATLAB第一章作业答案
信号与系统 matlab答案
matlab教程课后作业
matlab课后习题答案第四章
MATLAB基础教程薛山第二版课后习题答案讲解
(完整版)matlab习题及答案(2)
信号与系统作业
MATLAB教程 R2014a 答案 全 张志涌
信号与系统课后matlab作业.
中南大学matlab课后习题
汽车理论课后作业答案MATLAB
信号与系统MATLAB仿真题目
相关主题
文本预览