东华大学信号课程设计第二次Word格式.docx

东华大学信号课程设计第二次Word格式.docx

《东华大学信号课程设计第二次Word格式.docx》由会员分享，可在线阅读，更多相关《东华大学信号课程设计第二次Word格式.docx（15页珍藏版）》请在冰豆网上搜索。

）;

T=4;

tao=2;

A=1;

w=2*pi/T;

x=A*tao/T;

forn=1:

1:

N

x=x+（2*A*tao/T）*sin（n*w*tao/2）/（n*w*tao/2）*cos（n*w*t）;

end

plot（t,x）;

xlabel（'

Time（sec）'

title（['

num2str（N）]）;

（3）先令A=1，N=10，改变T和tao

N=10;

T=input（'

T='

tao=input（'

tao='

A=1;

w=2*pi/T;

n1=-N:

-1;

c1=（2*A*tao/T）*sin（n1*w*tao/2）/w/tao./n1*2;

c0=A*tao/T;

n2=1:

N;

c2=（2*A*tao/T）*sin（n1*w*tao/2）/w/tao./n1*2;

cn=[c1c0c2];

n=-N:

subplot（211）;

stem（n,abs（cn）,'

filled'

\omega/\omega_0'

title（'

Magnitudeofck'

subplot（212）;

stem（n,angle（cn）,'

Phaseofck'

T=4,tao=1的频谱图：

T=4,tao=2的频谱图：

T=4,tao=3的频谱图：

T=6,tao=1的频谱图：

T=20,tao=1的频谱图：

2.

（1）求傅立叶变换：

symsAttao;

f=A*（heaviside（t+tao/2）-heaviside（t-tao/2））;

F=fourier（f）

结果：

F=A*（（cos（（tao*w）/2）*i+sin（（tao*w）/2））/w-（cos（（tao*w）/2）*i-sin（（tao*w）/2））/w）

绘制频谱波形：

代码：

令A=1，观察tao变化对波形的影响：

symstw

taoinput（'

A=1；

X=int（A*exp（-j*w*t）,t,-tao/2,tao/2）;

ezplot（X,[-6*pi,6*pi]）;

gridon;

\omega'

ylabel（'

\Magnitude'

X（\omega）tao='

num2str（tao）]）;

3.

（2）

b=[0.0400];

a=[0.040.42];

[H,w]=freqs（b,a）;

plot（w,abs（H））;

set（gca,'

xtick'

[012345678910]）;

ytick'

[00.40.7071]）;

Magnitude'

|H（j\omega）|'

plot（w,angle（H））;

Phase'

\phi（\omega）'

图形：

4.

b=[01];

RC=input（'

RC='

a=[RC1];

freqs（b,a）;

RC=1

RC=10

RC=10^-1

RC=10^-2

RC=10-3

（2）e（t）的波形：

t=0:

0.001:

0.2;

et=cos（100*t）+cos（3000*t）;

plot（t,et）

取RC=1，则H（jw）=1/（1+j*w）;

求滤波后的函数表达式：

symst;

et=cos（100*t）+sin（3000*t）;

E=fourier（et）;

symsw;

H=1/（1+i*w）;

R=E*H;

rt=ifourier（R,t）

rt=（pi*exp（-t*100*i）*（1/10001+（100*i）/10001）+pi*exp（t*100*i）*（1/10001-（100*i）/10001）+pi*exp（-t*3000*i）*（-3000/9000001+i/9000001）+pi*exp（t*3000*i）*（-3000/9000001-i/9000001））/（2*pi）

绘制滤波后r（t）波形：

频率响应曲线：

即

（1）中RC=1的频率响应曲线