Matlab 正弦周期信号低通滤波器

淓郁q 2016-06-11 23:26:30 633 浏览

题目：设计不同的低通、高通滤波器，获取周期信号sin(n)的高频和低频谐波，分析其频谱希望各位大神帮我用Matlab做做这题，把代码给我。或者详细说说低通高通滤波器怎么搞。很着急，网上很多介绍滤波器的乱七八糟，看不懂。谢谢各位了

参与评论

登录后参与评论

全部评论(1条)

就上此网14 2016-06-12 00:00:00

对1Hz和400Hz的两个正弦信号进行低通滤波的matlab程序 ? 对1Hz和400Hz的两个正弦信号进行低通滤波的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 - 解决时间：2007-4-19 15:49 请问如何用Matlab将图像中的水平线干扰去除，正弦噪声应该用什么滤波器处理？具体怎样处理？提问者： sky1_fish - 一级Z佳答案 1.读图象 2.做fft 3.显示fft 4.肉眼观察频谱中的噪声.(正弦干扰是在频谱中是几个点) 5.频谱中有噪声的区域设置为0.(或者低通) 6.反fft变换 7.显示 1.读图象 im=imread('文件名'); 2.做fft im_fft=(fftshift(fft2(im)));%fftshift是把低频成分放到中间，习惯上要这样做。 3.显示fft im_magfft=abs(im_fft); temp1=log(im_magfft+1);%否则幅度太大，不能显示整个动态范围 Xmin = min(min(temp1)); Xmax = max(max(temp1)); scaf = 255/(Xmax-Xmin); im_ftd = floor(scaf*(temp1 - Xmin)); figure; imshow(im_ftd); 4.肉眼观察频谱中的噪声.(正弦干扰是在频谱中是几个点) 5.频谱中有噪声的区域设置为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 关于matlab在低通滤波器的设计... 设计一个让低于400pi以下频率的信号通过的数字低通滤波器。 ws，wp那些参数随便...我实在做不出来了....谢谢问题补充：或者请高手告诉我，比如我的低通滤波器要让400pihz一下的频率通过，那么我的ws，wp那一些参数究竟要如何选择或者计算呢？谢谢大家，分数不多，但是也请看看这个还得有个采样频率吧，比如根据奈氏定律，采样频率必须高于Z高频率的2倍，那么你可以设计采样频率800，或者1000等。此时数字频率就为400pi/1000=0.4pi，那么你可以设计wp=0.4pi，ws=0.45pi，还要确定一个Z小衰减，比如As=40db.这样就可以直接设计了。在MATLAB命令行输入fdatool，可以调出一个滤波器设计窗口。在那里面你会有很大的收获。例如你这个400PI的那么你的200HZ，采样频率为Fs为1K，Fpass设为200，Fstop设为400，衰减为80db后设计滤波器，然活倒出export，数据为lowpass1，命名利用，保存文件为*.mat，利用 load *; filter(lowpass1，1，data); 调用刚才生长的参数。data就为你要滤波的信号 MATLAB一阶低通滤波器的设计悬赏分：5 - 解决时间：2009-1-8 11:05 对f=5sin（2*pi*20t）+5sin（2*pi*60t），截止频率50Hz，设计MATLAB低通滤波器的程序及注释，谢谢提问者： feiyun0125 - 一级Z佳答案%滤波前数据 t = 0:0.01:10; f=5*sin(2*pi*20*t)+5*sin(2*pi*60*t); %滤波器设计 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滤波后的数据 f2 = filter(Hd，f); %对比图 plot(t，f，'b'，t，f2，'r'); 根据你的题目设计一个butterworth的3阶滤波器 N=3; [z，p，k]=buttap(N);求出该滤波器的零点和极点 OmegaC=2*pi*50; 给定的截止频率 p=p*OmegaC; 由于前面是归一化的，现在还原 k=k*OmegaC^N; >> k=k*OmegaC^N; B=real(poly(z)); b0=k; b=k*B; 求出的模拟滤波器的分子系数 a=real(poly(p)); 求出的模拟滤波器的分母系数 a，b都是按s的降幂排列的那么该题设计的滤波器是分子9.6139e+014 分母s^3+628s^2+197292s+30990544

赞(8)

回复(0)

评论

评论
登录后参与评论

热门问答

Matlab 正弦周期信号低通滤波器: 题目：设计不同的低通、高通滤波器，获取周期信号sin(n)的高频和低频谐波，分析其频谱希望各位大神帮我用Matlab做做这题，把代码给我。或者详细说说低通高通滤波器怎么搞。很着急，网上很多介绍滤波器的乱七八糟，看不懂。谢谢各位了

如何用MATLAB做一个低通滤波器？: 幅度|H|=1，0=
2009-12-02 11:02:16 351 3

用MATLAB设计一个低通滤波器: Z好解释一下各种参数怎么设置都是怎么计算出来的

关于matlab用IIR低通滤波器对信号进行滤波去噪问题: 我从网上搜的程序老师说滤波前后的频谱前后段对称、还有什么横轴要是什么坐标（完全听我也不懂）求大神看我这个程序怎么把频谱的后半对称部分去掉程序和图片在下面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/... 我从网上搜的程序老师说滤波前后的频谱前后段对称、还有什么横轴要是什么坐标（完全听我也不懂）求大神看我这个程序怎么把频谱的后半对称部分去掉程序和图片在下面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');%求低通滤波器的阶数和截止频率[b11，a11]=butter(n11，wn11，'s');%求S域的频率响应的参数[num11，den11]=bilinear(b11，a11，0.5);%双线性变换实现S域到Z域的变换[h，w]=freqz(num11，den11);%根据参数求出频率响应plot(w*8000*0.5/pi，abs(h));legend('用butter设计');grid;[y，fs，nbits]=wavread ('qzp1.wav');n=length (y);%求出语音信号的长度noise=0.01*randn(n，2);%随机函数产生噪声s=y+noise;%语音信号加入噪声S=fft(s);%傅里叶变换z11=filter(num11，den11，s);sound(z11);m11=fft(z11);%求滤波后的信号 subplot(2，2，1);plot(abs(S)，'g');title('滤波前信号的频谱');grid;subplot(2，2，2);plot(abs(m11)，'r');title('滤波后信号的频谱');grid;subplot(2，2，3);plot(s);title('滤波前信号的波形');grid;subplot(2，2，4);plot(z11);title('滤波后的信号波形'); 展开

matlab使用低通滤波器来处理数据: 想请教各位高人：对于一组10000个点的数据（接近正弦），如何用MATLAB来进行低通滤波，谢谢啦。 PS：10000个点的值为变压器原边一相的电流值

用remez设计FIR低通滤波器问题。matlab: 那位高手会设计这个滤波器，滤波器阶数P=50.通带[0，0.4pi].阻带[0.5pi，pi].通带Z大偏差delta1=0.01.阻带Z大偏差0.001.在[0，0.4pi]U[0.5pi，pi]间隔中等间隔采样N=200个点。谢谢！

MATLAB如何建立RC一阶低通滤波器模块:

关于matlab离散非周期信号的频谱分析一直有问题，麻烦大家帮帮忙。: 假设时间数列是t=[0，0.5，1.65，2.02，2.51，3.01，3.49，4.11，4.53]，然后对应的流强值为 I=[15，14.998，14.978，15.012，15.023，14.896，15.022，14.987，15.243]。即对应的信号为横轴是时间，纵轴是流强。对这个信号怎么进行频谱分析？我一直弄不懂的是时间怎么与流强... 假设时间数列是t=[0，0.5，1.65，2.02，2.51，3.01，3.49，4.11，4.53]，然后对应的流强值为 I=[15，14.998，14.978，15.012，15.023，14.896，15.022，14.987，15.243]。即对应的信号为横轴是时间，纵轴是流强。对这个信号怎么进行频谱分析？我一直弄不懂的是时间怎么与流强对应起来？展开

用MATLAB画正弦函数，如何分析频谱图？: >> 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);

帮我设计一个图像的低通滤波器，用MATLAB: 没有什么指标要求，在截止频率为50Hz，60Hz，70Hz，80Hz等频带上，观察不同图像的滤波效果，Z后确定Z优。设计出特定图像的低通滤波器，实现图像的低通滤波。谢谢大家了.

怎样用MATLAB画正弦函数以及怎样看频谱图？: 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))); 这个程序运行结果如图。我的问题1：为什么函数的周期是六百多？怎样画周期是2pi的正弦函数？2：第二张图怎样反应它的频谱呢？它的横轴和... 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))); 这个程序运行结果如图。我的问题1：为什么函数的周期是六百多？怎样画周期是2pi的正弦函数？2：第二张图怎样反应它的频谱呢？它的横轴和纵轴表示了什么呢？展开

如果周期信号的周期预先不知道，如何FFT进行频谱分析:

matlab中如何对正弦函数的采样序列进行频谱分析: 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（）进行FFT计算频谱，那括号中应该是什么？还有FFT之后如何将横坐标转换成频率。谢谢！

椭圆低通滤波器基于matlab设计程序如下，求每句注释: %连续信号的产生及采样 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('时间(s)'); ylabel('幅值'); grid; %椭圆低通滤波器的设计 [b，a]=ellip(5，0.1，40，5*2/Fs); [H，w]=freqz(b，a，51... %连续信号的产生及采样 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('时间(s)'); ylabel('幅值'); grid; %椭圆低通滤波器的设计 [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('频率响应图'); grid; %对滤波后的信号进行分析和变换 sf=filter(b，a，s); plot(t，sf); xlabel('时间 (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'}); 展开

非周期信号的频谱分析如何实现:

matlab中怎么让信号通过滤波器: %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 巴特沃夫滤波器 wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30; [N，wc]=buttord(wp，ws，Rp，As，'s')%计算率波器的阶数和3dB截止频率 [B，A]=butter(N，wc，'s');%计算滤波器系统函数分子分母多项式 fk=0:800/512:8000; wk=2*pi*f... %%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%% 巴特沃夫滤波器 wp=2*pi*2500;ws=2*pi*3500;Rp=1;As=30; [N，wc]=buttord(wp，ws，Rp，As，'s')%计算率波器的阶数和3dB截止频率 [B，A]=butter(N，wc，'s');%计算滤波器系统函数分子分母多项式 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]) %%%%%%%%%%%%%%%%%%%%% 信号为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 如何消除信号的高频分量:

matlab怎样分析语音信号的频谱:

怎么用matlab画出信号频谱分析图: 请问一下，怎么使用matlab将数据采集卡采集到的信号进行频谱分析，学校毕业设计要用，但matlab只用来做图像处理，没做过信号分析，进行信号频域分析，画出频谱图!!!!!

用MATLAB对信号并进行频谱分析: 1）查阅相关资料；（2）采用MATLAB产生的信号（三种以上频率，两种不同的采样间隔）；（3）对上述信号并进行频谱分析，并分析采样间隔、频率分辨率对频谱的影响；（4）在同一采样间隔下，将不同频率的信号求和，通过频谱分析，分别设计高通、低通、带通... 1）查阅相关资料；（2）采用MATLAB产生的信号（三种以上频率，两种不同的采样间隔）；（3）对上述信号并进行频谱分析，并分析采样间隔、频率分辨率对频谱的影响；（4）在同一采样间隔下，将不同频率的信号求和，通过频谱分析，分别设计高通、低通、带通、帯阻滤波器对信号进行滤波处理，实现某一频率信号的提取或滤除；展开

5月突出贡献榜