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 展开
全部评论(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; i
32768) { 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月突出贡献榜
推荐主页
最新话题
-
- #DeepSeek如何看待仪器#
- 干体炉技术发展与应用研究
- 从-70℃到150℃:一台试验箱如何终结智能...从-70℃到150℃:一台试验箱如何终结智能调光膜失效风险?解决方案:SMC-210PF-FPC温湿度折弯试验箱的五大核心价值1. 多维度环境模拟,覆盖全生命周期测试需求超宽温域:支持-70℃至+150℃的极限温度模拟(可选配),复现材料在极寒、高温、冷热冲击下的性能表现;控湿:湿度范围20%~98%RH(精度±3%RH),模拟热带雨林、沙漠干燥等复杂工况,暴露材料吸湿膨胀、分层缺陷;动态折弯:0°~180°连续可调折弯角度,支持R1~R20弯曲半径设定,模拟实际装配中的微小应力,提前预警裂纹、断裂风险。
参与评论
登录后参与评论