仪器网(yiqi.com)欢迎您!

| 注册 登录
网站首页-资讯-专题- 微头条-话题-产品- 品牌库-搜索-供应商- 展会-招标-采购- 社区-知识-技术-资料库-方案-产品库- 视频

问答社区

matlab中升余弦滤波器的程序不知道错在哪了?

忾謉 2013-06-29 11:32:43 399  浏览
  • % zsqEg7_7.m Fd=1e3; %码元时隙为1ms Fs=Fd*20; %在一个码元时隙内采样20次 delay=10; %滤波器延时为10个码元时隙宽度 for r=[0, 0.25, 0.5, 0.75, 1] %滚降系数为0,0.2... % zsqEg7_7.m Fd=1e3; %码元时隙为1ms Fs=Fd*20; %在一个码元时隙内采样20次 delay=10; %滤波器延时为10个码元时隙宽度 for r=[0, 0.25, 0.5, 0.75, 1] %滚降系数为0,0.25,0.5,0.75和1 num = rcosine(Fd,Fs, 'fir/normal',r,delay); t=0:1/Fs:1/Fs*(length(num)-1); figure(1); plot(t,num); axis([0 0.02 -0.3 1.1]);hold on; Hw=abs(fft(num,1000)); f=(1:Fs/1000:Fs)-1; figure(2); plot(f,Hw); axis([0 1500 0 22]);hold on; end 展开

参与评论

全部评论(1条)

  • FyfyY5H1h1 2013-06-30 00:00:00
    Fd=1e3; %码元时隙为1ms Fs=Fd*20; %在一个码元时隙内采样20次 delay=10; %滤波器延时为10个码元时隙宽度 r=[0,0.25,0.5,0.75,1] for i=1:length(r) %滚降系数为0,0.25,0.5,0.75和1 num = rcosine(Fd,Fs, 'fir/normal',r(i),delay); t=0:1/Fs:1/Fs*(length(num)-1); figure(1); plot(t,num); axis([0 0.02 -0.3 1.1]); hold on; Hw=abs(fft(num,1000)); f=(1:Fs/1000:Fs)-1; figure(2); plot(f,Hw); axis([0 1500 0 22]); hold on; end 应该这么写,for循环不能你这么用

    赞(12)

    回复(0)

    评论

获取验证码
我已经阅读并接受《仪器网服务协议》

热门问答

matlab中升余弦滤波器的程序不知道错在哪了?
% zsqEg7_7.m Fd=1e3; %码元时隙为1ms Fs=Fd*20; %在一个码元时隙内采样20次 delay=10; %滤波器延时为10个码元时隙宽度 for r=[0, 0.25, 0.5, 0.75, 1] %滚降系数为0,0.2... % zsqEg7_7.m Fd=1e3; %码元时隙为1ms Fs=Fd*20; %在一个码元时隙内采样20次 delay=10; %滤波器延时为10个码元时隙宽度 for r=[0, 0.25, 0.5, 0.75, 1] %滚降系数为0,0.25,0.5,0.75和1 num = rcosine(Fd,Fs, 'fir/normal',r,delay); t=0:1/Fs:1/Fs*(length(num)-1); figure(1); plot(t,num); axis([0 0.02 -0.3 1.1]);hold on; Hw=abs(fft(num,1000)); f=(1:Fs/1000:Fs)-1; figure(2); plot(f,Hw); axis([0 1500 0 22]);hold on; end 展开
2013-06-29 11:32:43 399 1
急求,matlab自适应滤波器的程序
 
2007-06-11 18:54:50 524 1
MATLAB中设计的滤波器程序是干什么用的?怎么用?
MATLAB中设计的滤波器程序是干什么用的?怎么用? 所有的类似的书中都有这些滤波器的设计~!但设计完这后这些程序怎么滤波啊? 那位高手指点一下啊~1
2007-04-10 15:15:54 467 3
matlab中示波器在哪?
初学matlab全是英文帮帮忙!!... 初学matlab 全是英文 帮帮忙!! 展开
2016-05-24 20:31:20 551 1
matlab中示波器在哪?
 
2018-11-21 18:58:05 471 0
matlab中设计fir滤波器
 
2011-11-30 15:44:48 290 3
matlab中电机模块在哪
 
2018-06-17 05:25:18 302 1
MATLAB中FIR滤波器如何设计
 
2017-09-03 07:45:45 347 1
基于RLS算法和LMS的自适应滤波器的MATLAB程序
 
2011-02-10 03:36:07 429 2
matlab fir滤波器
 
2011-09-08 05:04:49 303 3
matlab做一个滤波器
如题 随便怎样的 Z好有一定的文字说明 不要ctrl+c ctrl+v的 谢谢了
2016-11-11 10:14:37 350 1
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)相干载波 展开
2016-05-20 01:23:03 397 2
基于matlab设计线性相位fir滤波器的程序 急!万分感谢
 
2018-12-04 02:48:33 342 0
使用matlab中fdatool导出的高通butterworth IIR滤波器系数如何在程序中进行滤波
导出的系数为 const int NL = 10; const real64_T NUM[10] = { 0.8558457777315, -7.702611999584, 30.81044799834, -71.89104532945, 107.8365679942, -107.8365679942, 71.89104532945, -30.81044799834, 7.702611999584, ... 导出的系数为 const int NL = 10; const real64_T NUM[10] = { 0.8558457777315, -7.702611999584, 30.81044799834, -71.89104532945, 107.8365679942, -107.8365679942, 71.89104532945, -30.81044799834, 7.702611999584, -0.8558457777315 }; const int DL = 10; const real64_T DEN[10] = { 1, -8.688722481439, 33.55806880257, -75.61737954546, 109.553642522, -105.8293351142, 68.16455775223, -28.22860996292, 6.820250022448, -0.7324719952609 }; 在程序中进行差分计算得到滤波结果,如下: pWave为输入波形,pFilter的内容也复制为输入波形 for(i=16; i32768) { float_Value = 32768; } if (float_Value<-32768) { float_Value = -32768; } *(pFilter+i) = (float)(float_Value);//滤波 } 得到的结果波形不对。 请教问题在哪里? 展开
2013-07-18 13:57:09 466 2
matlab中fdesign滤波器filter design有什么区别
 
2017-03-19 11:04:38 412 1
matlab中想用filter函数实现IIR滤波器
 
2010-12-09 08:31:47 369 2
matlab怎么调用fdatool的滤波器
就是我用fdatool设计了一个滤波器 但是在matlab的主界面上怎么使用它 就是比如说有一个函数y=sin(x) 加入噪声之后的函数叫A,怎么让A通过这个滤波器得到我想要的波形 答得好加分
2017-09-26 04:37:28 362 2
用matlab设计一个滤波器
数字滤波器设计及在语音信号分析中的应用。 步骤: 1、语音信号采集 录制一段课程设计学生的语音信号并保存为文件,要求长度不小于10秒,并对录制的信号进行采样;录制时可以使用Windows自带的录音机,或者使用其它专业的录音软件,录制时需要配备录音硬件... 数字滤波器设计及在语音信号分析中的应用。 步骤: 1、语音信号采集 录制一段课程设计学生的语音信号并保存为文件,要求长度不小于10秒,并对录制的信号进行采样;录制时可以使用Windows自带的录音机,或者使用其它专业的录音软件,录制时需要配备录音硬件(如麦克风),为便于比较,需要在安静、干扰小的环境下录音。 2、语音信号分析 使用MATLAB绘出采样后的语音信号的时域波形和频谱图。根据频谱图求出其带宽,并说明语音信号的采样频率不能低于多少赫兹。 3、含噪语音信号合成 在MATLAB软件平台下,给原始的语音信号叠加上噪声,噪声类型分为如下几种:(1)白噪声;(2)单频噪色(正弦干扰);(3)多频噪声(多正弦干扰);(4)其它干扰,可设置为低频、高频、带限噪声,或Chirp干扰、冲激干扰。绘出叠加噪声后的语音信号时域和频谱图,在视觉上与原始语音信号图形对比,也可通过Windows播放软件从听觉上进行对比,分析并体会含噪语音信号频谱和时域波形的改变。 4、数字滤波器设计及滤波,完成以下题目中的一个 给定滤波器的规一化性能指标(参考指标,实际中依据每个同学所叠加噪声情况而定)例如:通带截止频率wp=0.25*pi, 阻通带截止频率ws=0.3*pi; 通带Z大衰减Rp=1 dB; 阻带Z小衰减Rs=15 dB,每个题目至少设计出5个用不同方法的不同类型滤波器。 题目(1):采用窗函数法与等波纹法分别设计各型FIR滤波器(低通、高通、带通、带阻中的至少3种类型)来对叠加噪声前后的语音信号进行滤波处理,绘出滤波器的频域响应,绘出滤波后信号的时域波形和频谱,并对滤波前后的信号进行对比,分析信号的变化;在相同的性能指标下比较各方法的滤波效果,并从理论上进行分析(或解释)。 题目(2):采用双线性变换法与脉冲响应不变法,分别利用不同的原型低通滤波器(Butterworth型与切比雪夫I型)来设计各型IIR滤波器(低通、高通、带通、带阻中的至少3种类型),绘出滤波器的频域响应;并用这些数字滤波器对含噪语音信号分别进行滤波处理,比较不同方法下设计出来的数字滤波器的滤波效果,并从理论上进行分析(或解释)。 5、回放语音信号 对滤波后的语音信号进行回放,感觉滤波前后语音信号的变化。 6、设计一个语音信号分析与处理系统界面(选作) 利用MATLAB的界面设计功能,设计一个MATLAB环境下的人机交互界面,完成上述滤波器的设计及滤波过程,要求能够接收输入的设计参数,查看设计结果及滤波结果。我的邮箱是471147021@QQ.com 哪位高手知道这方面的问题,给我点提示嘛,或者将答案发到我的邮箱里,谢谢了!!! 展开
2009-06-24 07:58:04 630 2
在Matlab中已知频率响应曲线如何设计滤波器
如题,我已经知道了一个频率响应,存在数组h中,怎样把这个设计成滤波器。
2009-04-09 19:53:24 348 1
滤波器matlab中wp=2*pi*fpass/fs意思
% IIR Lowpass Use Butterworth clear; fs=20; fpass=4; fstop=5; Ap=0.5; As=10; wp=2*pi*fpass/fs; ws=2*pi*fstop/fs; omegap=tan(wp/2); omegas=tan(ws/2); ep=sqrt(10^(Ap/10)-1); es=sqrt(10^(As/10)-1); N=ceil(log(es/e... % IIR Lowpass Use Butterworth clear; fs=20; fpass=4; fstop=5; Ap=0.5; As=10; wp=2*pi*fpass/fs; ws=2*pi*fstop/fs; omegap=tan(wp/2); omegas=tan(ws/2); ep=sqrt(10^(Ap/10)-1); es=sqrt(10^(As/10)-1); N=ceil(log(es/ep)/log(omegas/omegap)); omega0=omegap/ep^(1/N); K=floor(N/2); for i=1:K theta(i)=pi*(N-1+2*i)/(2*N); end for i=1:K G(i)=omega0^2/(1-2*omega0*cos(theta(i))+omega0^2); end for i=1:K a1(i)=2*(omega0^2-1)/(1-2*omega0*cos(theta(i))+omega0^2); end for i=1:K a2(i)=(1+2*omega0*cos(theta(i))+omega0^2)/(1-2*omega0*cos(theta(i))+omega0^2); end if K<(N/2) G0=omega0/(omega0+1); a0=(omega0-1)/(omega0+1); end w=0:pi/300:pi; Hw2=1./(1+(tan(w/2)/omega0).^(2*N)); plot(w/pi,Hw2); grid; 展开
2014-03-24 05:54:13 612 2

5月突出贡献榜

推荐主页

最新话题