MATLAB作业3.docx
- 文档编号:5039600
- 上传时间:2023-05-07
- 格式:DOCX
- 页数:11
- 大小:157.76KB
MATLAB作业3.docx
《MATLAB作业3.docx》由会员分享,可在线阅读,更多相关《MATLAB作业3.docx(11页珍藏版)》请在冰点文库上搜索。
MATLAB作业3
《信号与系统》课程研究性学习手册
本组成员王舟帆09271107
朱亚楠09211150
高智勇09211131
林怡09211135
指导教师魏杰
时间2010.11.21
专题研讨信号与系统的频域分析
研讨题目:
4-8
(1)k=1:
10;x=cos(0.2*pi*k);
X=fft(x);m=0:
9;stem(m,real(X));title('X(m)的实部');xlabel('m');figure;
stem(m,imag(X));title('X(m)的虚部')xlabel('m');
(2)
k=1:
2;x=cos(pi*k);X=fft(x);m=0:
1;
stem(m,real(X));
title('X(m)的实部');xlabel('m');figure;stem(m,imag(X));title('X(m)的虚部')xlabel('m');
4-4
(1)
functiony=sf1(t,w)
y=(t>=0&t<=2).*
(1).*exp(-1i*w*t);
w=linspace(-6*pi,6*pi,512);
N=length(w);X=zeros(1,N);
fork=1:
N
X(k)=quadl('sf1',0,2,[],[],w(k));
end
figure
(1);
plot(w,real(X));
(2)
functiony=sf2(t,w)
y=(t>=0&t<=2).*(cos(2*t)).*exp(-1i*w*t);
w=linspace(-6*pi,6*pi,512);
N=length(w);X=zeros(1,N);
fork=1:
N
X(k)=quadl('sf2',0,2,[],[],w(k));
end
figure
(1);
plot(w,real(X));
(3)
functiony=sf3(t,w)
y=(t>=0).*exp(-t).*exp(-1i*w*t);
w=linspace(-6*pi,6*pi,512);
N=length(w);X=zeros(1,N);
fork=1:
N
X(k)=quadl('sf3',0,100000,[],[],w(k));
end
figure
(1);
plot(w,real(X));
(4)
functiony=sf4(t,w)
y=exp(-li*w*t);
w=linspace(-6*pi,6*pi,512);
N=length(w);X=zeros(1,N);
fork=1:
N
X(k)=quadl('sf4',-1000,1000,[],[],(k));
end
figure
(1);
plot(w,real(X));
4-2
(1)
计算得P=0.25
t=-4:
0.001:
4;
N=3;
C0=0;
xN=C0*ones(1,length(t));
forn=1:
2:
N
xN=xN+sinc(n/2).^2.*n*pi*0.25.*sin(pi*n*t);
end
plot(t,xN);
title('N=3');
(2)
计算得P=0.33
t=-4:
0.001:
4;
N=1;
C0=0.5;
xN=C0*ones(1,length(t));
forn=1:
2:
N;
xN=xN+(-0.5*sinc(n/2).^2+sinc(n)).*cos(pi*n*t);
end
plot(t,xN);
title('N=1');
(3)
t=-4:
0.001:
4;
N=10;
C0=0;
xN=C0*ones(1,length(t));
forn=1:
2:
N
xN=xN+sinc(n/2).^2.*n*pi*0.25.*sin(pi*n*t);
end
plot(t,xN);
title('N=10');
N取10时的图形:
t=-4:
0.001:
4;
N=50;
C0=0;
xN=C0*ones(1,length(t));
forn=1:
2:
N
xN=xN+sinc(n/2).^2.*n*pi*0.25.*sin(pi*n*t);
end
plot(t,xN);
title('N=50');
N取50时的图形:
t=-4:
0.001:
4;
N=100;
C0=0;
xN=C0*ones(1,length(t));
forn=1:
2:
N
xN=xN+sinc(n/2).^2.*n*pi*0.25.*sin(pi*n*t);
end
plot(t,xN);
title('N=100');
N取100时的图形:
该实验验证了Gibbs现象。
4-1
(1)
实部
n=-20:
20;
Cn=-sinc(n/2).^2.*n*pi*0.25;
stem(n,Cn);
title('real(Cn)');
虚部
n=-20:
20;
Cn=0.5*pi.*ones(1,41);
stem(n,Cn);
title('imag(Cn)');
(2)
n=-10:
10;
Cn=-0.5*sinc(n/2).^2+sinc(n);
stem(n,Cn);
title('Cn');
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- MATLAB 作业
![提示](https://static.bingdoc.com/images/bang_tan.gif)