东华大学信号报告35.docx
- 文档编号:6008233
- 上传时间:2023-01-02
- 格式:DOCX
- 页数:26
- 大小:716.97KB
东华大学信号报告35.docx
《东华大学信号报告35.docx》由会员分享,可在线阅读,更多相关《东华大学信号报告35.docx(26页珍藏版)》请在冰豆网上搜索。
东华大学信号报告35
>>[H7,w7]=eg24(0.25);
>>[H8,w8]=eg24
(1);
>>[H9,w9]=eg24(5);
2.1
(2)
M文件“eg21”代码:
functionF=eg21(N,tao,T,A)
t=-5:
0.01:
5;
F=zeros(size(t));
forn=1:
N
F=F+2*A*sin(n*pi*tao/T)*cos(2*pi*n*t/T)/(n*pi);
end
plot(t,F);
axis([-5.55.5-0.70.7]);
xlabel('t');
title('N');
gridon;
命令窗口:
>>F1=eg21(10,2,4,1);
>>F2=eg21(100,2,4,1);
2.1(3)
M文件“eg212”代码:
functioneg212(N,tao,T,A)
n1=-N:
1:
-1;
c1=A./(pi.*n1).*sin(pi.*n1.*tao/T);
c0=A/T;
n2=1:
1:
N;
c2=A./(pi.*n2).*sin(pi.*n2.*tao/T);
cn=[c1c0c2];
n=-N:
1:
N;
subplot(211);
stem(n,abs(cn),'filled');
xlabel('\omega^omega_0');
title('Magnitudeofcn');
subplot(212);
stem(n,angle(cn),'filled');
xlabel('\omega^omega_0');
title('Phaseofcn');
命令窗口:
>>eg212(10,2,4,1);
>>eg212(10,2,8,1);
>>eg212(10,2,16,1);
>>eg212(10,4,16,1);
>>eg212(10,8,16,1);
3.1
(1)
>>b=[-11];
>>a=[11];
>>[H,w]=freqs(b,a);
>>subplot(211);
>>plot(w,abs(H));
>>set(gca,'xtick',[0:
4:
100]);
>>set(gca,'ytick',[00.40.7071]);
>>xlabel('\omega');
>>ylabel('Magnitude');
>>title('|H(j\omega)|');
>>gridon;
>>subplot(212);
>>plot(w,angle(H));
>>set(gca,'xtick',[0:
4:
100]);
>>set(gca,'ytick',[00.40.7071]);
>>xlabel('\omega');
>>ylabel('Phase');
>>title('\phi(\omega)');
>>gridon;
3.1
(2)
>>symstw;
>>H=(1-j*w)/(1+j*w);
>>e=sin(t)+sin(3*t);
>>E=fourier(e);
>>R=E*H;
>>r=ifourier(R);
>>t1=0:
0.0001:
10;
>>e1=sin(t1)+sin(3*t1);
>>subplot(221);
>>plot(t1,e1);
>>xlabel('t');
>>title('e=sin(t)+sin(3*t)');
>>subplot(222);
>>ezplot(r,[010]);
>>xlabel('t');
>>title('r(t)');
>>gridon;
>>b=[-11];
>>a=[11];
>>[H,w]=freqs(b,a);
>>subplot(223);
>>plot(w,abs(H));
>>set(gca,'xtick',[0:
10]);
>>set(gca,'ytick',[00.40.7071]);
>>xlabel('\omega');
>>ylabel('Magnitude');
>>title('|H(j\omega)|');
>>gridon;
>>subplot(224);
>>plot(w,angle(H));
>>set(gca,'xtick',[0:
10]);
>>set(gca,'ytick',[00.40.7071]);
>>xlabel('\omega');
>>ylabel('Phase');
>>title('\phi(\omega)');
>>gridon;
命令窗口:
>>eg321(500,0)
>>eg321(50,0)
>>eg321(500,100)
3.2
(1)
M文件“eg321”代码:
functioneg321(w,faii)
t=0:
0.001:
20;
x=sinc(t-10);
subplot(221);
plot(t,x);
xlabel('t');
title('原信号');
y=x.*sin(500*t);
subplot(222);
plot(t,y);
xlabel('t');
title('调制后的信号');
z=y.*sin(w.*t+faii);
subplot(223);
plot(t,z);
xlabel('t');
title('解调后的信号');
Rp=0.1;
Rs=80;
wp=40/250;
ws=45/250;
[n,wn]=ellipord(wp,ws,Rp,Rs);
[b,a]=ellip(n,Rp,Rs,wn);
ss=2*filter(b,a,z);
subplot(224);
plot(t,ss);
ylim([-1,1]);
xlabel('t');
title('滤波后的信号');
3.3
(1)
dt=0.005;
t=0:
dt:
10;
f=sinc(t);
subplot(221);
plot(t,f);
xlabel('t');
title('f(t)');
Ts=1;
ts=0:
Ts:
10;
fs=sinc(ts);
subplot(222);
plot(t,f,':
');
holdon;
stem(ts,fs);
xlabel('t');
title('f_s(t)');
w=linspace(-4*pi,4*pi,2001);
F=f.*exp(-j*t.*w)*dt;
subplot(223);
plot(w/pi,abs(F));
xlabel('\omega^pi');
title('|F(\omega)|');
w1=linspace(-4*pi,4*pi,11);
F1=fs*exp(-j*t.*w1);
subplot(224);
plot(w/pi,abs(F1));
xlabel('\omega^pi');
title('|F_s(\omega)|');
命令窗口:
>>eg332
(1);
>>eg332(0.5);
>>eg332(5);
3.3
(2)
M文件“eg332”代码:
functioneg332(Ts)
dt=0.005;
t=0:
dt:
10;
f=sinc(t);
ws=2*pi/Ts;
ts=0:
Ts:
10;
fs=sinc(ts);
subplot(311);
plot(t,f,':
');
holdon;
stem(ts,fs);
xlabel('t');
title('Sampling');
subplot(312);
holdon;
wc=ws/2;
fr=zeros(size(f));
fori=0:
Ts:
10;
sa=fs(i/Ts+1)*sinc(wc/pi*(t-i));
fr=fr+sa;
plot(t,sa,':
');
end;
plot(ts,fs,'o');
plot(t,fr);
xlabel('t');
title('Reconstruction');
err=abs(f-fr);
subplot(313);
plot(t,err);
xlabel('t');
title('Error');
M文件“eg333”代码:
functioneg333(a)
wm=1;
wc=a*wm;
Ts=a*pi/wm;
ws=2*pi/Ts;
n=-100:
100;
nTs=n*Ts;
f=sinc(nTs/pi);
Dt=0.005;
t=-20:
Dt:
20;
fa=f*Ts*wc/pi*sinc((wc/pi)*(ones(length(nTs),1)*t-nTs'*ones(1,length(t))));
error=abs(fa-sinc(t/pi));
t1=-20:
0.5:
20;
f1=sinc(t1/pi);
subplot(311);
stem(t1,f1);
xlabel('kTs');
ylabel('f(kTs)');
title('由sa(t)=sinc(t/pi)的采样信号sa(t)');
subplot(312);
plot(t,fa);
xlabel('t');
ylabel('fa(t)');
title('sa(t)=sinc(t/pi)的采样信号重构sa(t)');
gridon;
subplot(313);
plot(t,error);
xlabel('t');
ylabel('error(t)');
title('采样信号与原信号的误差error(t)');
命令窗口:
>>eg333
(1);
>>eg333(1.1);
>>eg333(2.5);
4.1
>>h=sym('heaviside(t)-heaviside(t-2)');
>>H=laplace(h);
>>e=sym('heaviside(t)');
>>E=laplace(e);
>>R=H*E;
>>r=ilaplace(R)
r=
t-heaviside(t-2)*(t-2)
4.2
(1)
>>b=[1];
>>a=[1221];
>>sys=tf(b,a);
>>[p,z]=pzmap(sys);
>>pzmap(sys);
4.2
(2)
>>clear
>>symss;
>>H=1/(s^3+2*s^2+2*s+1);
>>h=ilaplace(H)
h=
exp(-t)-exp(-t/2)*(cos((3^(1/2)*t)/2)-(3^(1/2)*sin((3^(1/2)*t)/2))/3)
>>ezplot(h)
4.3
>>symss;
>>b=[1];
>>a=[1217];
>>sys=tf(b,a);
>>[p,z]=pzmap(sys);
>>subplot(311);
>>pzmap(sys);%绘制极零图
>>holdon;
>>b=[18];
>>a=[1217];
>>sys=tf(b,a);
>>[p,z]=pzmap(b,a);
>>subplot(312);
>>pzmap(sys);
>>holdon;
>>b=[1-8];
>>a=[1217];
>>sys=tf(b,a);
>>[p,z]=pzmap(sys);
>>subplot(313);
>>pzmap(sys);
>>holdon;
>>symsst;
>>H1=1/(s^2+2*s+17);
>>H2=(s+8)/(s^2+2*s+17);
>>H3=(s-8)/(s^2+2*s+17);
>>h1=ilaplace(H1)
h1=
(sin(4*t)*exp(-t))/4
>>h2=ilaplace(H2)
h2=
exp(-t)*(cos(4*t)+(7*sin(4*t))/4)
>>h3=ilaplace(H3)
h3=
exp(-t)*(cos(4*t)-(9*sin(4*t))/4)
>>subplot(311);
>>ezplot(t,h1)
>>subplot(312);
>>ezplot(t,h2)
>>subplot(313);
>>ezplot(t,h3)
5.1
>>n=-3:
3;
>>x=[0123333];
>>n1=fliplr(n);
>>x1=fliplr(x);
>>n2=n-2;
>>n3=n+2;
>>subplot(221);
>>stem(n,x,'filled');
>>xlabel('n');
>>title('x(n)');
>>subplot(222);
>>stem(n1,x,'filled');
>>xlabel('n');
>>title('x(-n)');
>>subplot(223);
>>stem(n2,x,'filled');
>>xlabel('n');
>>title('x(n+2)');
>>subplot(224);
>>stem(n3,x,'filled');
>>xlabel('t');
>>title('x(n-2)');
5.2
>>n=-18:
18;
>>x=cos(n*pi/6)+sin(n*pi/3)+cos(n*pi/2);
>>stem(n,x,'filled');
>>title('x(n)');
5.3
>>b=[1];
>>a=[132];
>>impz(b,a);
5.4
>>x1=[01211];
>>x2=[11111];
>>x=conv(x1,x2);
>>n1=-2:
2;
>>n2=-2:
2;
>>subplot(311);
>>stem(n1,x1,'filled');
>>xlabel('n');
>>ylabel('x1(n)');
>>gridon;
>>subplot(312);
>>stem(n2,x2,'filled');
>>xlabel('n');
>>ylabel('x2(n)');
>>gridon;
>>n=-4:
4;
>>subplot(313);
>>stem(n,x,'filled');
>>xlabel('n');
>>ylabel('conv(x1,x2)');
>>gridon;
disp(x);
013455421
5.5
>>n=0:
40;
>>h=sin(n/2).*heaviside(n);
>>subplot(311);
>>stem(n,h,'filled');
>>xlabel('n');
>>ylabel('h(n)');
>>x=sin(n/5).*heaviside(n);
>>subplot(312);
>>stem(n,x,'filled');
>>xlabel('n');
>>ylabel('x(n)');
>>y=conv(h,x);
>>n=0:
80;
>>subplot(313);
>>stem(n,y,'filled');
>>xlabel('n');
>>ylabel('y(n)');
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 东华 大学 信号 报告 35