信号与系统程序文档格式.docx
- 文档编号:21045475
- 上传时间:2023-01-27
- 格式:DOCX
- 页数:26
- 大小:221.85KB
信号与系统程序文档格式.docx
《信号与系统程序文档格式.docx》由会员分享,可在线阅读,更多相关《信号与系统程序文档格式.docx(26页珍藏版)》请在冰豆网上搜索。
Q3
b=[1,-0.8];
a=[2,0,-1];
[H1omegal1]=freqz(a,b,4)
[H2omegal2]=freqz(a,b,4,'
whole'
)
H1=
5.0000
2.8200-1.3705i
1.8293-1.4634i
0.9258-0.9732i
omegal1=
0
0.7854
1.5708
2.3562
H2=
0.5556
1.8293+1.4634i
omegal2=
3.1416
4.7124
Q4
n=[0:
63]
a=[ones(1,8)];
x1=repmat(a,1,8);
stem(n,x1);
gridon;
holdon;
b=[ones(1,8)zeros(1,8)];
x2=repmat(b,1,4);
stem(n,x2);
c=[ones(1,8)zeros(1,24)];
x3=repmat(c,1,2);
stem(n,x3);
a1=(1/64)*fft(x1)
figure(4);
stem(n,real(a1));
xlabel('
Timeindexn'
);
ylabel('
Re(a1)'
a2=(1/64)*fft(x2)
figure(5);
stem(n,real(a2));
Re(a2)'
a3=(1/64)*fft(x3)
figure(6);
stem(n,real(a3));
Re(a3)'
figure(7);
stem(n,imag(a1));
Im(a1)'
figure(8);
stem(n,imag(a2));
Im(a2)'
figure(9);
stem(n,imag(a3));
Im(a3)'
Fugure(3)
Figure(4)
Figure(5)
Figure(6)
Figure(7)
Figure(8)
Figure(9)
forn=1:
64
x34(n)=0;
fori=1:
32
x34(n)=x34(n)+a3(i)*exp(j*2*i*pi*n/32);
end
end
y1=abs(x34);
stem(x34,y1)
Afternoon
Q5
functiony=fuliye(t,n)
y=0;
fork=-n:
-1
y=y+exp(j*pi*k*t/2)*(sin(pi/2*k)/(k*pi));
y=y+.5
fork=1:
n
y=2*y-1;
x=-2.5:
0.001:
2,5;
n=input('
pleaseentern:
'
y1=(fuliye(x,n)+1)/2;
T=-2.5:
.001:
y=(square((pi/2)*T+(pi/2))+1)/2;
plot(T,y,x,y1)
axis([-2.52.5-13]);
一阶
三阶
五阶
七阶
九阶
Q6
b=[1,-2];
a=[1,1.5,0.5];
freqs(b,a)
[r,p]=residue(b,a)
r=
6
-5
p=
-1.0000
-0.5000
Q7
b1=[1,5];
a1=[1,2,3];
b2=[2,5,12];
a2=[1,2,10];
b3=[2,5,12];
a3=[1,4,14,20];
b4=[1,2,5];
a4=[1,-3];
zs1=roots(b1);
ps1=roots(a1);
zs2=roots(b2);
ps2=roots(a2);
zs3=roots(b3);
ps3=roots(a3);
zs4=roots(b4);
ps4=roots(a4);
figure
(1)
plot(real(zs1),imag(zs1),'
o'
grid
holdon
figure
(2)
plot(real(ps1),imag(ps1),'
x'
figure(3)
plot(real(zs2),imag(zs2),'
figure(4)
plot(real(ps2),imag(ps2),'
figure(5)
plot(real(zs3),imag(zs3),'
figure(6)
plot(real(ps3),imag(ps3),'
figure(7)
plot(real(zs4),imag(zs4),'
figure(8)
plot(real(ps4),imag(ps4),'
Q8
a=[1-0.90.81];
b=[100];
zplane(b,a);
title('
零、极点分布图'
omega=[0:
511]*pi/256;
[H,omega]=freqz(b,a,omega);
unitcirc=exp(j*omega);
ps1=roots(a);
polevectors=[repmat(unitcirc,2,1)-repmat(ps1,1,512)];
polelength=abs(polevectors);
poleangle=angle(polevectors);
figure
(2),plot(omega'
polelength);
polelength-omega'
figure(3),plot(omega'
poleangle);
poleangle-omega'
ps2=roots(b)
zerovectors=[repmat(unitcirc,2,1)-repmat(ps2,1,512)];
zerolength=abs(zerovectors);
zeroangle=angle(zerovectors);
figure(4),plot(omega'
zerolength);
zerolength-omega'
figure(5),plot(omega'
zeroangle);
zeroangle-omega'
figure(6)
subplot(211),plot(w'
polelength1(:
1).*polelength1(:
2))
频率响应的模'
)
subplot(212),plot(w'
abs(H'
))
频率响应的相角'
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 信号 系统 程序