MATLAB 2013信号与系统实验考试卷
考试形式：开卷 时间：60分钟（请确保在此之前将程序存到下面指定的地方）
要求：在网络盘Z盘新建一文件夹，名为：学号+姓名（如：060092**张*）。确保考试结束前将考试程序（.m文件）和大作业（代码和报告）存在此目录下。不收移动硬盘等。
请在一个程序内完成以下内容，并要求具有良好的可读性：
1、 编写脚本合成下面信号，并在一帧内显示k=5和k=20 合成的两个波形。
2、 将k=20 时合成的波形通过系统函数为下式的LTI系统，请画出响应波形。
3、 请自己编写脚本（不调用库函数）画出(2)中系统的频率特性曲线。
4、 利用符号运算，求（2）中系统的单位阶跃响应，并画出响应波形。
分析：
函数
function y=f(t,k)
y=0.5;
if k==5
for i=1:5
y=y+cos(i*pi-1)*cos(10*i*t)/(i*i*pi*pi)-cos(i*pi)*sin(10*i*t)/(i*pi);
end
end
if k==20
for i=1:20
y=y+cos(i*pi-1)*cos(10*i*t)/(i*i*pi*pi)-cos(i*pi)*sin(10*i*t)/(i*pi);
end
end
function r=eksm(t)
if t<0
r=0
else
r=1;
end
主程序：
%第1题
t=0:0.01:2;%t的取点
y1=f(t,5);%k=5
y2=f(t,20);%k=20
figure(1);%画第一题的图形
subplot(2,1,1);
plot(t,y1);%k=5时的图形
subplot(2,1,2);
plot(t,y2);%k=20时的图形
%第2题
num=[15,0];
den=[1,10,1];
sys=tf(num,den);%生成系统函数
e=f(t,20);%激励函数是k=20的f(t)
figure(2);%第二题的图形
lsim(sys,e,t);
%第3题
for n=1:400
w(n)=(n-1)*0.05;
H(n)=(15*1j*w(n))/((1j*w(n))^2+10*1j*w(n)+1);
end
figure(3);%第三题的图形
mag=abs(H);%幅频特性
phase=angle(H);%相频特性
subplot(2,1,1);
plot(w,mag);title('幅频特性');
subplot(2,1,2);
plot(w,phase);title('相频特性');
%第4题
figure(4);%第四题的图形
step(sys);%直接取阶跃响应
（注：第4题应用系统函数求）
本文来源：https://www.2haoxitong.net/k/doc/5ad231e95022aaea998f0f92.html
文档为doc格式