Matlab 正弦周期信號(hào) 低通濾波器

淓郁q 2016-06-11 23:26:30 700 瀏覽

題目：設(shè)計(jì)不同的低通、高通濾波器，獲取周期信號(hào)sin(n)的高頻和低頻諧波，分析其頻譜希望各位大神幫我用Matlab做做這題，把代碼給我。或者詳細(xì)說說低通高通濾波器怎么搞。很著急，網(wǎng)上很多介紹濾波器的亂七八糟，看不懂。謝謝各位了

參與評(píng)論

評(píng)論

登錄后參與評(píng)論

全部評(píng)論(1條)

就上此網(wǎng)14 2016-06-12 00:00:00

對(duì)1Hz和400Hz的兩個(gè)正弦信號(hào)進(jìn)行低通濾波的matlab程序 ? 對(duì)1Hz和400Hz的兩個(gè)正弦信號(hào)進(jìn)行低通濾波的matlab程序 fs=1000; t=0:1/fs:5; x=cos(2*pi*1*t)+cos(2*pi*400*t); plot(t，x); grid for n=2:length(x); x(n)=0.9999*x(n-1)+x(n); end figure(2); plot(t，x) grid matlab 去除噪聲懸賞分：0 - 解決時(shí)間：2007-4-19 15:49 請(qǐng)問如何用Matlab將圖像中的水平線干擾去除，正弦噪聲應(yīng)該用什么濾波器處理？具體怎樣處理？提問者： sky1_fish - 一級(jí)Z佳答案 1.讀圖象 2.做fft 3.顯示fft 4.肉眼觀察頻譜中的噪聲.(正弦干擾是在頻譜中是幾個(gè)點(diǎn)) 5.頻譜中有噪聲的區(qū)域設(shè)置為0.(或者低通) 6.反fft變換 7.顯示 1.讀圖象 im=imread('文件名'); 2.做fft im_fft=(fftshift(fft2(im)));%fftshift是把低頻成分放到中間，習(xí)慣上要這樣做。 3.顯示fft im_magfft=abs(im_fft); temp1=log(im_magfft+1);%否則幅度太大，不能顯示整個(gè)動(dòng)態(tài)范圍 Xmin = min(min(temp1)); Xmax = max(max(temp1)); scaf = 255/(Xmax-Xmin); im_ftd = floor(scaf*(temp1 - Xmin)); figure; imshow(im_ftd); 4.肉眼觀察頻譜中的噪聲.(正弦干擾是在頻譜中是幾個(gè)點(diǎn)) 5.頻譜中有噪聲的區(qū)域設(shè)置為0.(或者低通) %下面舉例低通 r = 70; %低通的ZX半徑，自己觀察 for i=1:m for j=1:n if (((i-m/2)^2+(j-n/2)^2) < r^2 ) mask(i，j) = 1; else mask(i，j) = 0; end end end im2_fft=im_fft.*mask; 6.反fft變換 im2=ifft(ishiftfft(im2_fft)); 7.顯示 imshow(im2); 參考資料：D.Renshaw， Edinburgh University 關(guān)于matlab在低通濾波器的設(shè)計(jì)... 設(shè)計(jì)一個(gè)讓低于400pi以下頻率的信號(hào)通過的數(shù)字低通濾波器。 ws，wp那些參數(shù)隨便...我實(shí)在做不出來了....謝謝問題補(bǔ)充：或者請(qǐng)高手告訴我，比如我的低通濾波器要讓400pihz一下的頻率通過，那么我的ws，wp那一些參數(shù)究竟要如何選擇或者計(jì)算呢？謝謝大家，分?jǐn)?shù)不多，但是也請(qǐng)看看這個(gè)還得有個(gè)采樣頻率吧，比如根據(jù)奈氏定律，采樣頻率必須高于Z高頻率的2倍，那么你可以設(shè)計(jì)采樣頻率800，或者1000等。此時(shí)數(shù)字頻率就為400pi/1000=0.4pi，那么你可以設(shè)計(jì)wp=0.4pi，ws=0.45pi，還要確定一個(gè)Z小衰減，比如As=40db.這樣就可以直接設(shè)計(jì)了。在MATLAB命令行輸入fdatool，可以調(diào)出一個(gè)濾波器設(shè)計(jì)窗口。在那里面你會(huì)有很大的收獲。例如你這個(gè)400PI的那么你的200HZ，采樣頻率為Fs為1K，F(xiàn)pass設(shè)為200，F(xiàn)stop設(shè)為400，衰減為80db后設(shè)計(jì)濾波器，然活倒出export，數(shù)據(jù)為lowpass1，命名利用，保存文件為*.mat，利用 load *; filter(lowpass1，1，data); 調(diào)用剛才生長(zhǎng)的參數(shù)。data就為你要濾波的信號(hào) MATLAB一階低通濾波器的設(shè)計(jì) 懸賞分：5 - 解決時(shí)間：2009-1-8 11:05 對(duì)f=5sin（2*pi*20t）+5sin（2*pi*60t），截止頻率50Hz，設(shè)計(jì)MATLAB低通濾波器的程序及注釋，謝謝提問者： feiyun0125 - 一級(jí)Z佳答案%濾波前數(shù)據(jù) t = 0:0.01:10; f=5*sin(2*pi*20*t)+5*sin(2*pi*60*t); %濾波器設(shè)計(jì) Fpass = 1; % Passband Frequency Fstop = 50; % Stopband Frequency Apass = 1; % Passband Ripple (dB) Astop = 60; % Stopband Attenuation (dB) Fs = 2000; % Sampling Frequency h = fdesign.lowpass('fp，fst，ap，ast'， Fpass， Fstop， Apass， Astop， Fs); Hd = design(h， 'equiripple'， ... 'MinOrder'， 'any'， ... 'StopbandShape'， 'flat'); %f2濾波后的數(shù)據(jù) f2 = filter(Hd，f); %對(duì)比圖 plot(t，f，'b'，t，f2，'r'); 根據(jù)你的題目設(shè)計(jì)一個(gè)butterworth的3階濾波器 N=3; [z，p，k]=buttap(N);求出該濾波器的零點(diǎn)和極點(diǎn) OmegaC=2*pi*50; 給定的截止頻率 p=p*OmegaC; 由于前面是歸一化的，現(xiàn)在還原 k=k*OmegaC^N; >> k=k*OmegaC^N; B=real(poly(z)); b0=k; b=k*B; 求出的模擬濾波器的分子系數(shù) a=real(poly(p)); 求出的模擬濾波器的分母系數(shù) a，b都是按s的降冪排列的那么該題設(shè)計(jì)的濾波器是分子9.6139e+014 分母s^3+628s^2+197292s+30990544

贊(8)

回復(fù)(0)

評(píng)論

評(píng)論
登錄后參與評(píng)論

登錄或新用戶注冊(cè)

微信登錄
密碼登錄
短信登錄

請(qǐng)用手機(jī)微信掃描下方二維碼
快速登錄或注冊(cè)新賬號(hào)

微信掃碼，手機(jī)電腦聯(lián)動(dòng)

注冊(cè)登錄即表示同意《儀器網(wǎng)服務(wù)條款》和《隱私協(xié)議》

熱門問答

Matlab 正弦周期信號(hào) 低通濾波器: 題目：設(shè)計(jì)不同的低通、高通濾波器，獲取周期信號(hào)sin(n)的高頻和低頻諧波，分析其頻譜希望各位大神幫我用Matlab做做這題，把代碼給我。或者詳細(xì)說說低通高通濾波器怎么搞。很著急，網(wǎng)上很多介紹濾波器的亂七八糟，看不懂。謝謝各位了

如何用MATLAB做一個(gè)低通濾波器？: 幅度|H|=1，0=
2009-12-02 11:02:16 441 3

用MATLAB設(shè)計(jì)一個(gè)低通濾波器: Z好解釋一下各種參數(shù)怎么設(shè)置都是怎么計(jì)算出來的

關(guān)于matlab用IIR低通濾波器對(duì)信號(hào)進(jìn)行濾波去噪問題: 我從網(wǎng)上搜的程序老師說濾波前后的頻譜前后段對(duì)稱、還有什么橫軸要是什么坐標(biāo)（完全聽我也不懂）求大神看我這個(gè)程序怎么把頻譜的后半對(duì)稱部分去掉程序和圖片在下面Ft=8000;Fp=1000;Fs=1200;wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;fp=2*Ft*tan(wp/2);fs=2*Fs*tan(wp/... 我從網(wǎng)上搜的程序老師說濾波前后的頻譜前后段對(duì)稱、還有什么橫軸要是什么坐標(biāo)（完全聽我也不懂）求大神看我這個(gè)程序怎么把頻譜的后半對(duì)稱部分去掉程序和圖片在下面Ft=8000;Fp=1000;Fs=1200;wp=2*pi*Fp/Ft;ws=2*pi*Fs/Ft;fp=2*Ft*tan(wp/2);fs=2*Fs*tan(wp/2);[n11，wn11]=buttord(wp，ws，1，50，'s');%求低通濾波器的階數(shù)和截止頻率[b11，a11]=butter(n11，wn11，'s');%求S域的頻率響應(yīng)的參數(shù)[num11，den11]=bilinear(b11，a11，0.5);%雙線性變換實(shí)現(xiàn)S域到Z域的變換[h，w]=freqz(num11，den11);%根據(jù)參數(shù)求出頻率響應(yīng)plot(w*8000*0.5/pi，abs(h));legend('用butter設(shè)計(jì)');grid;[y，fs，nbits]=wavread ('qzp1.wav');n=length (y);%求出語音信號(hào)的長(zhǎng)度noise=0.01*randn(n，2);%隨機(jī)函數(shù)產(chǎn)生噪聲s=y+noise;%語音信號(hào)加入噪聲S=fft(s);%傅里葉變換z11=filter(num11，den11，s);sound(z11);m11=fft(z11);%求濾波后的信號(hào) subplot(2，2，1);plot(abs(S)，'g');title('濾波前信號(hào)的頻譜');grid;subplot(2，2，2);plot(abs(m11)，'r');title('濾波后信號(hào)的頻譜');grid;subplot(2，2，3);plot(s);title('濾波前信號(hào)的波形');grid;subplot(2，2，4);plot(z11);title('濾波后的信號(hào)波形'); 展開

matlab使用低通濾波器來處理數(shù)據(jù): 想請(qǐng)教各位高人：對(duì)于一組10000個(gè)點(diǎn)的數(shù)據(jù)（接近正弦），如何用MATLAB來進(jìn)行低通濾波，謝謝啦。 PS：10000個(gè)點(diǎn)的值為變壓器原邊一相的電流值

用remez設(shè)計(jì)FIR低通濾波器問題。matlab: 那位高手會(huì)設(shè)計(jì)這個(gè)濾波器，濾波器階數(shù)P=50.通帶[0，0.4pi].阻帶[0.5pi，pi].通帶Z大偏差delta1=0.01.阻帶Z大偏差0.001.在[0，0.4pi]U[0.5pi，pi]間隔中等間隔采樣N=200個(gè)點(diǎn)。謝謝！

MATLAB如何建立RC一階低通濾波器模塊:

關(guān)于matlab離散非周期信號(hào)的頻譜分析一直有問題，麻煩大家?guī)蛶兔Α?: 假設(shè)時(shí)間數(shù)列是t=[0，0.5，1.65，2.02，2.51，3.01，3.49，4.11，4.53]，然后對(duì)應(yīng)的流強(qiáng)值為 I=[15，14.998，14.978，15.012，15.023，14.896，15.022，14.987，15.243]。即對(duì)應(yīng)的信號(hào)為橫軸是時(shí)間，縱軸是流強(qiáng)。對(duì)這個(gè)信號(hào)怎么進(jìn)行頻譜分析？我一直弄不懂的是時(shí)間怎么與流強(qiáng)... 假設(shè)時(shí)間數(shù)列是t=[0，0.5，1.65，2.02，2.51，3.01，3.49，4.11，4.53]，然后對(duì)應(yīng)的流強(qiáng)值為 I=[15，14.998，14.978，15.012，15.023，14.896，15.022，14.987，15.243]。即對(duì)應(yīng)的信號(hào)為橫軸是時(shí)間，縱軸是流強(qiáng)。對(duì)這個(gè)信號(hào)怎么進(jìn)行頻譜分析？我一直弄不懂的是時(shí)間怎么與流強(qiáng)對(duì)應(yīng)起來？展開

用MATLAB畫正弦函數(shù) ，如何分析頻譜圖？: >> clear >> N=1000; >> fs=100; >> n=0:N-1; >> t=n/fs; >> y=3*sin(2*t); >> plot(y); >> y1=fft(y，N); >> y2=abs(y1); >> plot(y2);

幫我設(shè)計(jì)一個(gè)圖像的低通濾波器，用MATLAB: 沒有什么指標(biāo)要求，在截止頻率為50Hz，60Hz，70Hz，80Hz等頻帶上，觀察不同圖像的濾波效果，Z后確定Z優(yōu)。設(shè)計(jì)出特定圖像的低通濾波器，實(shí)現(xiàn)圖像的低通濾波。謝謝大家了.

怎樣用MATLAB畫正弦函數(shù)以及怎樣看頻譜圖？: clear all; x=0:0.01:4*pi; y1=sin(x); figure() subplot(2，1，1) plot(y1) grid on subplot(2，1，2) stem(abs(fft(y1))); 這個(gè)程序運(yùn)行結(jié)果如圖。我的問題1：為什么函數(shù)的周期是六百多？怎樣畫周期是2pi的正弦函數(shù)？2：第二張圖怎樣反應(yīng)它的頻譜呢？它的橫軸和... clear all; x=0:0.01:4*pi; y1=sin(x); figure() subplot(2，1，1) plot(y1) grid on subplot(2，1，2) stem(abs(fft(y1))); 這個(gè)程序運(yùn)行結(jié)果如圖。我的問題1：為什么函數(shù)的周期是六百多？怎樣畫周期是2pi的正弦函數(shù)？2：第二張圖怎樣反應(yīng)它的頻譜呢？它的橫軸和縱軸表示了什么呢？展開

如果周期信號(hào)的周期預(yù)先不知道，如何FFT進(jìn)行頻譜分析:

matlab中如何對(duì)正弦函數(shù)的采樣序列進(jìn)行頻譜分析: clear; Fs=100; N=64; t=(0:1/Fs:1); y=sin(2*pi*10*t); stem(t(1:N)，y(1:N)); 用上面的語句采樣，之后用Y=fft（）進(jìn)行FFT計(jì)算頻譜，那括號(hào)中應(yīng)該是什么？還有FFT之后如何將橫坐標(biāo)轉(zhuǎn)換成頻率。謝謝！

橢圓低通濾波器基于matlab設(shè)計(jì)程序如下，求每句注釋: %連續(xù)信號(hào)的產(chǎn)生及采樣 clear Fs=100; t=(1:100)/Fs; s1=sin(2*pi*t*5); s2=sin(2*pi*t*15); s3=sin(2*pi*t*30); s=s1+s2+s3; plot(t，s); xlabel('時(shí)間(s)'); ylabel('幅值'); grid; %橢圓低通濾波器的設(shè)計(jì) [b，a]=ellip(5，0.1，40，5*2/Fs); [H，w]=freqz(b，a，51... %連續(xù)信號(hào)的產(chǎn)生及采樣 clear Fs=100; t=(1:100)/Fs; s1=sin(2*pi*t*5); s2=sin(2*pi*t*15); s3=sin(2*pi*t*30); s=s1+s2+s3; plot(t，s); xlabel('時(shí)間(s)'); ylabel('幅值'); grid; %橢圓低通濾波器的設(shè)計(jì) [b，a]=ellip(5，0.1，40，5*2/Fs); [H，w]=freqz(b，a，512); plot(w*Fs/(2*pi)，abs(H)); xlabel('頻率(Hz)');ylabel('頻率響應(yīng)圖'); grid; %對(duì)濾波后的信號(hào)進(jìn)行分析和變換 sf=filter(b，a，s); plot(t，sf); xlabel('時(shí)間 (s)'); ylabel('幅值'); axis([0 1 -1 1]); S=fft(s，512); SF=fft(sf，512); w=(0:255)/256*(Fs/2); plot(w，abs([S(1:256)' SF(1:256)'])); xlabel('頻率(Hz)'); ylabel('傅立葉變換圖'); grid; legend({'before'，'after'}); 展開

非周期信號(hào)的頻譜分析如何實(shí)現(xiàn):

matlab中怎么讓信號(hào)通過濾波器: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 巴特沃夫?yàn)V波器 wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30; [N，wc]=buttord(wp，ws，Rp，As，'s')%計(jì)算率波器的階數(shù)和3dB截止頻率 [B，A]=butter(N，wc，'s');%計(jì)算濾波器系統(tǒng)函數(shù)分子分母多項(xiàng)式 fk=0:800/512:8000; wk=2*pi*f... %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 巴特沃夫?yàn)V波器 wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30; [N，wc]=buttord(wp，ws，Rp，As，'s')%計(jì)算率波器的階數(shù)和3dB截止頻率 [B，A]=butter(N，wc，'s');%計(jì)算濾波器系統(tǒng)函數(shù)分子分母多項(xiàng)式 fk=0:800/512:8000; wk=2*pi*fk; Hk=freqs(B，A，wk); figure plot(fk/1000，20*log10(abs(Hk))); grid on，xlabel('頻率（kHz）')，ylabel('幅度（dB）') title('低通濾波器') axis([0，4，-35，5]) %%%%%%%%%%%%%%%%%%%%% 信號(hào)為m=[ones(1，t0/(3*ts))，-2*ones(1，t0/(3*ts))，zeros(1，t0/(3*ts)+1)].*cos(2*pi*fc.*t).*cos(2*pi*fc.*t); cos(2*pi*fc.*t)相干載波展開

matlab 如何消除信號(hào)的高頻分量:

matlab怎樣分析語音信號(hào)的頻譜:

怎么用matlab畫出信號(hào)頻譜分析圖: 請(qǐng)問一下，怎么使用matlab將數(shù)據(jù)采集卡采集到的信號(hào)進(jìn)行頻譜分析，學(xué)校畢業(yè)設(shè)計(jì)要用，但matlab只用來做圖像處理，沒做過信號(hào)分析，進(jìn)行信號(hào)頻域分析，畫出頻譜圖!!!!!

用MATLAB對(duì)信號(hào)并進(jìn)行頻譜分析: 1）查閱相關(guān)資料；（2）采用MATLAB產(chǎn)生的信號(hào)（三種以上頻率，兩種不同的采樣間隔）；（3）對(duì)上述信號(hào)并進(jìn)行頻譜分析，并分析采樣間隔、頻率分辨率對(duì)頻譜的影響；（4）在同一采樣間隔下，將不同頻率的信號(hào)求和，通過頻譜分析，分別設(shè)計(jì)高通、低通、帶通... 1）查閱相關(guān)資料；（2）采用MATLAB產(chǎn)生的信號(hào)（三種以上頻率，兩種不同的采樣間隔）；（3）對(duì)上述信號(hào)并進(jìn)行頻譜分析，并分析采樣間隔、頻率分辨率對(duì)頻譜的影響；（4）在同一采樣間隔下，將不同頻率的信號(hào)求和，通過頻譜分析，分別設(shè)計(jì)高通、低通、帶通、帯阻濾波器對(duì)信號(hào)進(jìn)行濾波處理，實(shí)現(xiàn)某一頻率信號(hào)的提取或?yàn)V除；展開

4月突出貢獻(xiàn)榜

推薦主頁

Matlab 正弦周期信號(hào) 低通濾波器

聯(lián)系我們

關(guān)注我們