求救cot 阶比重采样 问题
%-----数据导入(手动修改)--------------------------------------------%--------------------------------------------------------------------------
N=2200000;
Fs=40000;
order_max=48;
t=(0:N-1)*1/Fs;
array_time_amp=Data1; %时域振动信号导入
pluse=Data2; %脉冲信号导入
%--------------------------------------------------------------------------
%-----计算脉冲发生时刻---------------------------------------------
%--------------------------------------------------------------------------
Num_pluse1=1;
threshold=200; %设定脉冲阈值
for temp2=1:length(pluse)-1;
if (pluse(temp2)<200&&pluse(temp2+1)>=200)
Num_pluse1=;
end
end
if length(Num_pluse1)<2, return,end;
t_pluse=(Num_pluse1-1)/Fs;
%----------------------------------------------------------
%-------数字跟踪滤波-----------------------------------
%-----------------------------------------------------------
i=2;
while i<=length(t_pluse);
ft=2/(t_pluse(i)-t_pluse(i-1));
=ellip(6,3,80,ft*order_max/(Fs/2));
array_time_amp(Num_pluse1(i-1):(Num_pluse1(i)-1))=filter(b,a,array_time_amp(Num_pluse1(i-1):(Num_pluse1(i)-1)));
i=i+1;
end
%----------------------------------------------------------
%---------重采样----------------------------------------
%----------------------------------------------------------
%------等角度时间计算--------------------------------
delt_thet=pi/24; %设定重采样角度
t_angle=[];
for temp3=3:length(t_pluse);
b=inv([1,t_pluse(temp3-2),t_pluse(temp3-2)^2;1,t_pluse(temp3-1),t_pluse(temp3-1)^2;1,t_pl
use(temp3),t_pluse(temp3)^2])*';
if temp3==3;
k=0;
while k<1.5*2*pi/delt_thet;
tt=(sqrt(4*b(3)*(k*delt_thet-b(1))+b(2)^2)-b(2))/(2*(b(3)+eps));
t_angle=;
k=k+1;
end
else
k=pi/delt_thet;
while k>=pi/delt_thet && k<1.5*2*pi/delt_thet;
tt=(sqrt(4*b(3)*(k*delt_thet-b(1))+b(2)^2)-b(2))/(2*(b(3)+eps));
t_angle=;
k=k+1;
end
end
end
array_angle=.*delt_thet;
%------等角度时刻的幅值拟合（三种拟合方式）-----
%------线性插值拟合-----------------------------
% array_angle_amp=[];
% for temp5=1:length(t_angle);
% for temp4=1:length(t);
% if t_angle(temp5)>=t(temp4)&&t_angle(temp5)<t(temp4+1);
% %linear interpolation
% y=;
% x=;
% angle_amp=interp1(x,y,t_angle(temp5),'linear');
% array_angle_amp=;
%
% % -------- 三次多项式插值拟合------------------
% % if temp4==1;
% % y=;
% % x=;
% % angle_amp=interp1(x,y,t_angle(temp5),'linear');
% % array_angle_amp=;
% % else
% %
y=;
% %由前两点与后两点得用三次多项式插值求得此角度对应的幅值
% % x=;
% % angle_amp=interp1(x,y,t_angle(temp5),'cubic');
% % array_angle_amp=;
% % end
% break;
% end
% end
% end
%----------三次样条插值拟合-------------------------------
array_angle_amp=interp1(t,array_time_amp,t_angle,'cubic');
%--------------------------------------------------------
%------------------加窗--------------------------------
%-------------------------------------------------------
w=hanning(length(array_angle_amp))';
array_angle_amp=array_angle_amp.*w;
s%---------------------------------------------------------------
%-------------- FFT ------------------------------------------
%---------------------------------------------------------------
angle_dom_ffty=abs(fft(array_angle_amp))*2/length(array_angle_amp);
delt_order=2*pi/(length(angle_dom_ffty)*delt_thet);
angle_dom_fx=(0:length(angle_dom_ffty)-1)*delt_order;
FFTy=abs(fft(array_time_amp))*2/length(array_time_amp);
fx=(0:length(array_time_amp)-1)/t_sample;
%----------------------------------------------------------------
%---------------2 维图形显示----------------------------------
%----------------------------------------------------------------
figure(1);
subplot(2,1,1),plot(t,array_time_amp),title('timedomianat'),xlabel('time'),ylabel('amplitude');
grid on;
subplot(2,1,2),plot(t,pluse),title('keyphasor'),xlabel('time'),ylabel('amplitude');
grid on;
figure(2);
subplot(2,1,1),plot(fx(1:length(fx)/2),FFTy(1:length(fx)/2)),title('Frequencydomain');
xlabel('f'), ylabel('amplitude');
subplot(2,1,2),plot(array_angle,array_angle_amp), title('angle dominant'),xlabel('angle /rad'),
ylabel('amplitude');
grid on;
figure(3);
subplot(2,1,1),plot(angle_dom_fx(1:length(angle_dom_fx)/2),angle_dom_ffty(1:length(angle_dom_fx)/2));file:///C:/Users/Administrator.L1PAC6T64O1PUEG/Documents/Tencent%20Files/895628164/Image/85Y3B%7B5~EK96Z%609LXLYB@%60M.jpg
但是程序运行出现了这样的结果file:///C:/Users/Administrator.L1PAC6T64O1PUEG/Documents/Tencent%20Files/895628164/Image/Z2A~GC76E%60ZD0@MEKBRW4
我不知道是不是程序有问题
threshold=200; %设定脉冲阈值，这个是怎么设定的
order_max=48;这个是计算出来的吗还是固定数值
整个程序还有哪些需要设定参数的地方
谢谢各位大神了 我是菜鸟
这个程序有点问题，你matlab图上显示的数据是哪来的？ impulse 发表于 2014-3-18 17:24
这个程序有点问题，你matlab图上显示的数据是哪来的？
自己实验室测的能帮我改一下吗 谢谢 最近也在做这个 {:{33}:} 现在成高手了吧！{:{27}:} lbtv 发表于 2014-10-10 16:10
现在成高手了吧！
还行哈哈哈哈哈哈哈哈哈哈哈 impulse 发表于 2014-3-18 17:24
这个程序有点问题，你matlab图上显示的数据是哪来的？
impluse 这个程序哪块出问题了？ 怎么老显示 ??? Error using ==> times Matrix dimensions must agree.调试了好久也没调出 矩阵维数是哪错了？ 谢谢啦 t_pluse 为什么是0？ 369 发表于 2014-12-2 16:14
t_pluse 为什么是0？
这个程序错好多 卡索 发表于 2014-12-2 16:17
这个程序错好多
能加下qq交流下么？ 我的是2899355008 369 发表于 2014-12-3 15:17
能加下qq交流下么？ 我的是2899355008
楼主问题解决了么，可否交流下 369 发表于 2014-12-3 15:17
能加下qq交流下么？ 我的是2899355008
楼主问题解决了么，可否交流下 楼主问题解决了吗
各位，麻烦请教仿真信号怎么产生的 impulse 发表于 2014-3-18 17:24
这个程序有点问题，你matlab图上显示的数据是哪来的？
主任，怎么仿真阶比的脉冲信号啊，谢谢
页:
[1]
2