马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。
您需要 登录 才可以下载或查看,没有账号?我要加入
x
:@Q 1》》选择合适的窗函数设计一个FIR数字低通滤波器,要求:通带截至频率为w =0.2 ,R =0.25dB;阻带截至频率为w =0.3 ,A =50dB。求脉冲响应并描绘一张滤波器的频率响应图
wp=0.2*pi;ws=0.3*pi; %输入设计指标
tr_width=ws-wp; %计算过渡带的宽度
M=ceil(6.6*pi/ tr_width)+1; %由海明窗的过渡带宽度选择合适的阶数
n=[0:1:M-1];
wc=(ws+wp)/2; %理想低通的截止频率(取平均值)
hd=ideal_lp(wc,M); %理想低通滤波器冲激响应
w_ham=(haming(M))’;
h=hd .* w_ham; %给理想低通加海明窗
[db,mag,pha,grd,w]= freqz_m(h,[1]); %求出幅频响应
delta_w=2*pi/1000;
Rp=-(min(db(1:1:wp/delta_w+1))); %检验通带波动
As=-round(max(db(ws/delta_w+1:1:501)));%检验最小阻带衰减
% plots(绘图)
subplot﹙2,2,1﹚,stem﹙n,hd﹚; title('Ideal Impulse Response');
axis﹙0 M-1 -0.1 0.3﹚;xlabel﹙'n'﹚;ylabel('hd(n)');
subplot﹙2,2,2﹚;stem﹙n,w_ham);title('Hamming Window') ;
axis﹙[0 M-1 0.1 1]﹚; xlabel﹙'n'﹚;ylabel('w(n)');
subplot﹙2,2,3﹚;stem﹙n,h﹚; title('Actual Impulse Response');
axis﹙[0 M-1 -0.1 0.3]﹚; xlabel﹙'n'﹚;ylabel('h(n)');
subplot﹙2,2,4﹚;plot(w/pi,db);title('Magnitude Response in dB');grid
axis﹙[0 1 -100 10]﹚; xlabel﹙'frequency in pi units'﹚;
ylabel('Decibels');
??? Undefined command/function 'ideal_lp'.
2》》用MATLAB的关系运算式来产生单位阶跃序列u﹙n﹚。
n1=-2;n2=8;n0=0;
n=n1:n2; %生成离散信号的时间序列
x=[n>=n0]; %生成离散信号x﹙n﹚
stem﹙n,x, 'filled'﹚;
axis﹙[n1,n2,0,1.1*max﹙x﹚]﹚;
title﹙'单位阶跃序列'﹚;
xlabel﹙'时间﹙n﹚'﹚;ylabel﹙'幅度x﹙n﹚'﹚;
??? stem﹙n,x, 'filled'﹚;
|
Error: Missing variable or function.
3》》编写产生a=1/2和a=2实指数连续信号和离散序列的程序
n1=-10;n2=10;a1=0.5;a2=2;
na1=n1:0;x1=a1.^na1;
na2=0:n2;x2=a2.^na2;
subplot﹙2,2,1﹚,plot﹙na1,x1﹚;
title﹙'实指数原信号﹙a<1﹚'﹚;
subplot﹙2,2,3﹚,stem﹙na1,x1, 'filled'﹚;
title﹙'实指数原信号﹙a<1﹚'﹚;
subplot﹙2,2,2﹚,plot﹙na2,x2﹚;
title﹙'实指数原信号﹙a>1﹚'﹚;
subplot﹙2,2,4﹚, stem﹙na2,x2, 'filled'﹚;
title﹙'实指数原信号﹙a>1﹚'﹚;
??? na1=n1:0;x1=a1.^na1;
|
Error: Missing variable or function.
还有很多这样的例子,好久都没弄明白!不管怎么弄,就是不出图xiang!!!! |