信号与线性系统分析
实验报告
学院：xxxxxxxxxxxxxxx
班级： xxxxxxxxxxxxxx
学号： xxxxxxxxxxxx
姓名： xxxxxxxx
2011-12-13
实验一1. 产生-10<t<10之间的序列和f=sint/t的值
代码：
t=-10:0.5:10;
f=sin(t)./t;
plot(t,f)
波形图：
2. 绘制阶跃函数的波形
定义阶跃函数代码：
function f=heaviside(t);
f=(t>0);
调用阶跃函数代码：
f=heaviside(t);
plot(t,f)
axis([-1,3,-0.2,1.2])
阶跃波形图：
3.画出f=exp(-2*t) .*heaviside(t).
代码：
f=exp(-2*t) .*heaviside(t);
plot(t,f)
axis([-1,5,-0.1,0.4])
波形图：
3. 正弦函数程序函数单数代码：t=-pi:pi/40:pi;
f=sin(2*pi*50*t); plot(t,f)
axis([-3,3,-1.5,1.5]) 波形图：
实验二
连续信号的时域描述与运算
一.信号的平移和反转
1.将函数u（t）=heaviside(t);
代码：
function f=u(t);
f=heaviside(t);
2.画出f(t)=t*[u(t)-u(t-1)]
代码：
f=t.*[u(t)-u(t-1)];
plot(t,f)
axis([-3,3,-0.1,1.2])
波形图：
定义initialsignal（t）= t*[u(t)-u(t-1)]; 代码：
function f=initialsignal(t);
f=t.*[u(t)-u(t-1)];
波形的平移和反转过程：
代码：
t=-2:0.01:2;
f=initialsignal(t);
subplot(231)
plot(t,f)
f1=initialsignal(t+1);
subplot(232)
plot(t,f1)
f2=initialsignal(t-1); subplot(233)
plot(t,f2)
f3=initialsignal(-t); subplot(234)
plot(t,f3)
f4=initialsignal(-t+1); subplot(235)
plot(t,f4)
f5=initialsignal(-t-1); subplot(236)
plot(t,f5)
波形图：
二．信号的尺度变换
代码定义f=initialsignal(t)；
代码：
function f=initialsignal(t);
f（t）=t.*[(u(t)-u(t-1))+(u(t-1)-u(t-2))]; 波形initialsignal(t);
2.绘出f（-2t+2）
代码：
t=-3:0.01:3;
f=initialsignal(t); subplot(231)
plot(t,f)
f1=initialsignal(t+2); subplot(232)
plot(t,f1)
f2=initialsignal(2*t+2); subplot(233)
plot(t,f2)
f3=initialsignal(2*t); subplot(234)
plot(t,f3)
f4=initialsignal(-2*t);
subplot(235)
plot(t,f4)
f5=initialsignal(-2*t+2);
subplot(236)
plot(t,f5)
波形图：
三．卷积积分
F(t)=f1(t)*f2(t)=[u(t)-u(t-1)]*[u(t)-u(t-1)]代码：
符号运算方法代码；
syms tao;
t=sym('t','positive');
ft1=sym('heaviside(t)-heaviside(t-1)'); ft2=sym('heaviside(t)-heaviside(t-1)'); ft_tao=subs(ft1,t,tao)*subs(ft2,t,t-tao); ft=int(ft_tao,tao,0,t);
ft=simplify(ft);
ezplot(ft,[0 2]);
grid on
数值运算方法：
t=-1:0.01:3;
f1=u(t)-u(t-1);
f=0.01*conv(f1,f1);
tmin=-2;
tmax=6;
t1=tmin:0.01:tmax;
plot(t1,f)
grid on
波形图：
实验三
离散信号的时域描述与运算一：1.单边指数序列f（n）=a^n*u(n)
定义u（n）
代码: n=0:10;
a1=1.2;a2=-1.2;a3=0.8;a4=-0.8;
f1=a1.^n;f2=a2.^n;f3=a3.^n;f4=a4.^n;
subplot(221)
stem(n,f1,'fill');xlabel('n');grid on
subplot(222)
stem(n,f2,'fill');xlabel('n');grid on
subplot(223)
stem(n,f3,'fill');xlabel('n');grid on subplot(224)
stem(n,f4,'fill');xlabel('n');grid on 序列图：
二．序列的平移,反转，尺度变换1.序列插零值源程序：
clf;
n=0:50;
x=sin(2*pi*0.12*n);
y=zeros(1,3*length(x));
y([1:3:length(y)])=x;
subplot(211)
subplot(212)
m=0:3*length(x)-1;
stem(m,y,'.')
波形图：
2.序列线性插值源代码
clf;
L=input('intput the interp number L,L='); n=0:49;
x=sin(2*pi*0.12*n);
y=interp(x,L);
subplot(211)
subplot(212)
m=0:50*L-1;
stem(m,y(1:50*L),'.')
运行：intput the interp number L,L=2
序列图：
实验四
连续信号的频域分析1. 周期信号的三角形式的傅里叶级数
源程序：先定义fourierseries函数
function y=fourierseries(m,t);
y=0.5;
for n=1:m
y=y+2*sin(n*pi/2)/(n*pi).*cos(n*pi.*t); end
调用函数画图代码：
t=-6:0.01:6;
d=-6:2:6;
fxx=pulstran(t,d,'rectpuls');
f1=fourierseries(3,t);
f2=fourierseries(9,t);
f3=fourierseries(21,t);
f4=fourierseries(45,t);
subplot(221)
plot(t,fxx,'r',t,f1,'b');
grid on
axis([-6 6 -0.1 1.1]);
subplot(222)
plot(t,fxx,'r',t,f2,'b');
grid on
axis([-6 6 -0.1 1.1]);
subplot(223)
plot(t,fxx,'r',t,f3,'b');
grid on
axis([-6 6 -0.1 1.1]); subplot(224)
plot(t,fxx,'r',t,f4,'b'); grid on
axis([-6 6 -0.1 1.1]);。