用雙線性變換法設計巴特沃斯低通IIR數字濾波器。設計指標引數為 在通帶內頻率低於0 2時,最大衰減

時間 2021-08-31 18:47:51

1樓:匿名使用者

fp=100hz,fs=300hz,αp=3db,αs=20db,抽樣頻率為fs=1000hz

fp=100;fs=300;fs=1000;

rp=3;rs=20;

wp=2*pi*fp/fs;

ws=2*pi*fs/fs;

fs=fs/fs %使fs為一wap=tan(wp/2);

was=tan(ws/2);

[n,wn]=buttord(wap,was,rp,rs,'s')[z,p,k]=buttap(n);

[bp,ap]=zp2tf(z,p,k)

[bs,as]=lp2lp(bp,ap,wap)[bz,az]=bilinear(bs,as,fs/2)[h,w]=freqz(bz,az,256,fs*1000);

2樓:手機使用者

1;as=15;

%buttord butter 預設的是採用雙線性變換法

[n,wc]=buttord(wp,ws,rp,as);

[b,a]=butter(n,wc);

wk=0:0.02:pi/2;

subplot(311)

plot(wk,20*log10(abs(freqz(b,a,wk))))

ylabel('db');

xlabel('w(0-pi/2)');

xn=[-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,...

-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,...

0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0] ;

subplot(312)

plot(linspace(0,pi,1024),abs(fft(xn,1024)));

title('原訊號頻譜');

xlabel('w(0-pi)');

subplot(313)

yn=filter(b,a,xn);

plot(linspace(0,pi,1024),abs(fft(yn,1024)))

title('通過低通濾波器後的頻譜');

xlabel('w(0-pi)');

輸入序列 為:

=[-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16

-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0

0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0]

**:%x(n)的心電脈衝函式

x=[-4,-2,0,-4,-6,-4,-2,-4,-6,-6,-4,-4,-6,-6,-2,6,12,8,0,-16,-38,-60,-84,-90,-66,-32,-4,-2,-4,8,12,12,10,6,6,6,4,0,0,0,0,0,-2,-4,0,0,0,-2,-2,0,0,-2,-2,-2,-2,0];

subplot(2,2,1);

n=0:55;

stem(n,x,'.');

xlabel('n');

ylabel('x(n)');

title('x(n)的心電脈衝函式');

%通過濾波器h1(z)後的y1(n)函式

a=0.09036;

b1=[a,2*a,a];

a1=[1,-1.2686,0.7051];

h1=filter(b1,a1,x);

[h1,w]=freqz(b1,a1,100);

%通過濾波器h1(z)、h2(z)後的y2(n)函式

b2=[a,2*a,a];

a2=[1,-1.0106,0.3583];

h2=filter(b2,a2,h1);

[h2,w]=freqz(b2,a2,100);

%通過濾波器h1(z)、h2(z)、h3(z)後的y3(n)函式

b3=[a,2*a,a];

a3=[1,-0.9044,0.2155];

h3=filter(b3,a3,h2);

[h3,w]=freqz(b3,a3,100);

subplot(2,2,2);

stem(n,h3,'.');

xlabel('n');

ylabel('y(n)');

title('通過濾波器h1(z)、h2(z)、h3(z)後的y3(n)函式');

輸出結果:

利用matlab結合雙線性變換法設計一個iir數字butterworth低通濾波器

3樓:匿名使用者

這麼專業的問題去專業論壇去問吧

4樓:

在命令視窗輸入

wp=1k/10k

ws=1.5/10

rp=1.5db

rs=25

[n,wn]=buttord(wp,ws,rp,rs)[b,a]=butter(n,wn)

freqz(b,a)搞定

高等代數線性變換,高等代數關於線性變換的問題!

高等代數是代數學發展到高階階段的總稱,它包括許多分支。如今大學裡開設的高等代數,一般包括兩部分 線性代數初步 多項式代數。很多人把高等代數和線性代數混為一談,但其實高等代數是大學數學專業開設的專業課,線性代數是大學中除了數學專業以外的理科,工科和部分醫科專業開設的課程。代數學 幾何學 分析數學是數學...

線性變換的意義,“線性變換”的物理含義有哪些?

豆村長de草 線性變換的意義 把線性對映寫成具體而簡明的2維數陣形式後,就成了一種矩陣。進而由線性對映的加法規則和複合規則來分別定義矩陣的加法規則和乘法規則是很自然的想法。當空間的基變化 座標系變換 時,線性對映的矩陣也會有規律地變化。在特定的基上研究線性對映,就轉化為對矩陣的研究。利用矩陣的乘法,...

什麼是非退化線性變換?怎麼做

假面 非退化矩陣就是行列式不等於零。若n階矩陣a的行列式 a 0,n階方陣a是非退化的充要條件為a是可逆矩陣。一個n n矩陣是非退化的充要條件是它的秩等於n。設a,b都是數域f上的n n矩陣,矩陣ab為退化的充要條件是a,b中至少有一個是退化的。 延濃綺 比如有一題,4x1x2 2x1x3 2x2x...