
请教信号处理高手!!!小波变换除噪相关问题
最近做用到用小波变换分析频谱极值点对应的fmax,初始信号是静电传感器采集到的电压信号,离线分析,导入数据以矩阵形式,在matlab里面,先进行FFT,再进行小波变换,请...
最近做用到用小波变换分析频谱极值点对应的fmax,初始信号是静电传感器采集到的电压信号,离线分析,导入数据以矩阵形式,在matlab里面,先进行FFT,再进行小波变换,请教各位大侠应该用哪个小波处理命令,是一维还是二维?最好提供matlab代码,以供借鉴!小弟感谢不尽!
展开
1个回答
展开全部
load jzxh.txt
s1=jzxh;
fs=500000; %采样频率
t=1:max(size(s1));
t=1000*t./fs;
t=t';
plot(s1),title('plot(s1)'),grid on;
%axis([0,51.2,-10000,10000]);
xlabel('时间(t) 单位:ms');
ylabel('电压(V)');
title('原始信号');
%小波分解
[c,l]=wavedec(s1,4,'db4');
[cd1,cd2,cd3,cd4]=detcoef(c,l,[1 2 3 4]);
ca4=appcoef(c,l,'db4',4);
a4=wrcoef('a',c,l,'db4',4);
d4=wrcoef('d',c,l,'db4',4);
d3=wrcoef('d',c,l,'db4',3);
d2=wrcoef('d',c,l,'db4',2);
d1=wrcoef('d',c,l,'db4',1);
%不需要的系数置零
caa4=zeros(size(ca4));
%消噪后信号重构,并画图
cl=[caa4',cd4',cd3',cd2',cd1']
s2=waverec(cl,l,'db4');
figure(2);
plot(t,s2);title('消噪后重构信号');grid on;
给你个例子,自己对这改改
s1=jzxh;
fs=500000; %采样频率
t=1:max(size(s1));
t=1000*t./fs;
t=t';
plot(s1),title('plot(s1)'),grid on;
%axis([0,51.2,-10000,10000]);
xlabel('时间(t) 单位:ms');
ylabel('电压(V)');
title('原始信号');
%小波分解
[c,l]=wavedec(s1,4,'db4');
[cd1,cd2,cd3,cd4]=detcoef(c,l,[1 2 3 4]);
ca4=appcoef(c,l,'db4',4);
a4=wrcoef('a',c,l,'db4',4);
d4=wrcoef('d',c,l,'db4',4);
d3=wrcoef('d',c,l,'db4',3);
d2=wrcoef('d',c,l,'db4',2);
d1=wrcoef('d',c,l,'db4',1);
%不需要的系数置零
caa4=zeros(size(ca4));
%消噪后信号重构,并画图
cl=[caa4',cd4',cd3',cd2',cd1']
s2=waverec(cl,l,'db4');
figure(2);
plot(t,s2);title('消噪后重构信号');grid on;
给你个例子,自己对这改改
本回答被提问者采纳
已赞过
已踩过<
评论
收起
你对这个回答的评价是?

2025-08-07 广告
广州赛恩科学仪器有限公司(原中大科仪)始创于2064年,是全球领先的精密测量仪器供应商和微弱信号检测方案提供商。公司以锁相放大器为核心产品,陆续推出光学斩波器、源表、功率放大器、电化学工作站、电流源等一系列产品。赛恩科仪推出的锁相放大器,覆...
点击进入详情页
本回答由赛恩科仪提供
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询