信号与系统-信号与系统课程实习报告.doc

文档介绍：HUNAN UNIVERSITY
课程实****报告
题目: 信号与系统
学生姓名梁晴
学生学号 20110801128
专业班级 11级计科一班
指导老师舒维星
完成日期 2013年12月15日
实验一信号的描述:连续离散
1、连续信号
代码:
f=sym('heaviside(t)');
ezplot(f,[-5,5]);
2、离散序列x(n)=RN(n),N=5
代码:
k=0:4;
x=[1 3 2 8 5];
stem(k,x)
实验二卷积
计算以下信号的卷积
代码:
function[f,k] = sconv(f1,f2,k1,k2,p)
%计算连续信号卷积积分f(t)=f1(t)*f2(t)
%f:卷积积分f(t)的对应的非零样值向量;
%k:f(t)对应时间向量
%f1:f1(t) 非零样值向量;
%f2:f2(t) 非零样值向量;
%k1:f(1) 对应时间向量
%k2:f(2) 对应时间向量
%p:取样时间间隔
f = conv(f1,f2);
f = f*p;
k0 = k1(1)+k2(1);
k3 = length(f1)+length(f2)-2;
k = k0:p:k3*p;
subplot 311;
plot(k1-2,f1);
title('f1(t)')
xlabel('t')
axis([-,,-,]);
f2=stepfun(t,0)-stepfun(t,4);
%subplot(2,1,2),plot(t,f2);
ylabel('f1(t)')
subplot 312
plot(k2-2,f2)
title('f2(t)')
xlabel('t')
ylabel('f2(t)')
subplot 313
plot(k-4,f)
h = get(gca,'position');
h(3) = *h(3);
set(gca,'position');
title('f(t)=f1(t)*f2(t)');
xlabel('t')
ylabel('f(t)')
输入参数:p=;
t=0::4;
f1=2*(stepfun(t,1)-stepfun(t,3))
% subplot(2,1,1),plot(t,f1);
axis([-,,-,]);
[f,k]=sconv(f1,f2,t,t,p);
截图:
实验三吉布斯现象
例1Sc(n)=1/2 + 2/pi*cos(pi*n) - 2/3/pi*cos(3*pi*n) + 2/5/pi*cos(5*pi*n)-2/7/pi*cos(7*pi*n)…
代码:
Dfunction WavConstruction
clear
clc
close all
DeT = 1/128*4;
n = -1:DeT:1;
nTao = -1/2:DeT:1/2;
sTao(1:length(nTao)) = ones(1, length(nTao));
s(1:length(n)) = [zeros(1, 1/2*length(nTao)) sTao zeros(1, 1/2*length(nTao))];
plot(n, s)
xlim([- ])
ylim([- ])
figure
sc1 = 1/2 + 2/pi*cos(pi*n);
subplot(2,2,1)
plot(n, sc1)
title('N = 1')
subplot(2,2,2)
sc2 = 1/2 + 2/pi*cos(pi*n) - 2/3/pi*cos(3*pi*n);
plot(n, sc2)
title('N = 3')
subplot(2,2,3)
sc3 = 1/2 + 2/pi*cos(pi*n) - 2/3/pi*cos(3*pi*n) + 2/5/pi*cos(5*pi*n);
plot(n, sc3)
title('N = 5')
subplot(2,2,4)
sc4 =1/2 + 2/pi*cos(pi*n) - 2/3/pi*cos(3*pi*n) + 2/5/pi*cos(5*pi*n) - 2/pi/7*cos(7*pi*n);
plot(n, sc4)
title('N = 7')
figure
subplot(2,1,1)
sc5 = 1/2 + 2/pi*cos(pi*n) - 2/3/pi*cos(3*pi*n) + 2/5/pi*cos(