matlab wp,matlab程序求助,不懂怎么定义没有定义的函数第一个：wp=0.2*pi;ws=0.3*pi;tr_width=ws-wp;M=ceil(6.6*pi/tr_width)+1;n...

matlab程序求助,不懂怎么定义没有定义的函数第一个：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=(hamming(M))";h=hd.*w_ham;freqz (h,[1])figure(2);subplot

2016-08-08

matlab程序求助,不懂怎么定义没有定义的函数

第一个：

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=(hamming(M))';

h=hd.*w_ham;

freqz (h,[1])

figure(2);

subplot(2,2,1),stem(n,hd);title('理想脉冲响应')

axis([0 M-1 -0.3 0.3]);xlabel('n');ylabel('hd(n)')

xa=0.*n;

hold on

plot(n,xa,'k');

hold off

subplot(2,2,2),stem(n,w_ham);title('hamming 窗')?

axis([0 M-1 -0.3 1.2]);xlabel('n');ylabel('w(n)')

subplot(2,2,3),stem(n,h);title('实际脉冲响应')

axis([0 M-1 -0.3 0.3]);xlabel('n');ylabel('h(n)')

hold on

plot(n,xa,'k');

hold off

第二个：

ws1=0.2*pi; wp1=0.35*pi;

ws2=0.65*pi;wp2=0.8*pi;

As=60;

tr_width=min((wp1-ws1),(wp2-ws2));

M=ceil(11*pi/tr_width)+1;

n=[0:1:M-1];

wc1=(ws1+wp1)/2;

wc2=(ws2+wp2)/2;

hd=ideal_lp(wc2,M)-ideal_lp(wc1,M); %需要定义

w_bla=(blackman(M))';

h=hd.*w_bla

freqz(h,[1])

figure(2);

subplot(2,2,1),stem(n,hd);title('idael?impulse?response')?

axis([0 M-1 -0.4 0.5]);xlabel('n');ylabel('hd(n)')

xa=0.*n;

hold on

plot(n,xa,'k');

hold off