matlab对多组数据拟合,
一内容介绍时频分析作为非平稳信号处理的一种重要方法，它提供了时域与频域的联合分布信息，清楚的刻画了信号频率随时间的变化关系，能给出每个时刻下信号的瞬时频率及其幅值。本文把时频分析作为研究对象，在对传统时频分析方法做出深入分析的基础上，重点研究同步压缩连续小波变换的时频分析方法，并将其应用到解决多分量非平稳信号的时频聚集性、分离以及瞬时频率估计的问题上。
2部分代码清楚的
load(vib_data1.mat )
fs=12000N=1200%采样频率和点数
时间=(1:N)/fs；%时间顺序
fre=(fs/2)/(N/2):(fs/2)/(N/2):(fs/2)；%频率序列
数据=数据(1:N)；
[TFR山]=MSST _ Y _新（数据，95,50);
伽马=10^(-2);
西格玛=0.025;
[~，~，~，~，~，Ts1，Ts2，Ts3，Ts4]=sstn(data，gamma，sigma)；
.........
数字
suptitle(图19 );
支线剧情(211)
情节（时间，数据);
ylabel(安培/克);
xlabel( Time/s )；
轴([0 0.1-4 4]);
子情节(212)
ft=abs(fft(数据))/600;
plot(fre，ft(1:end/2))；
ylabel(安培/克);
xlabel(频率/赫兹);
轴([0 6000 0 0.3]);
..........
x1=0.06x2=0.065
y1=2.6y2=3.7
dd=0.1
数字
suptitle(图20 );
支线剧情(511);
imagesc(时间，fre/1000，ABS(TFR))；
%xlabel(时间/秒);
ylabel( Fre/kHz )；
正常男性染色体组型轴；轴([0 0.1 y1-DD y2 DD])；
矩形(位置，[x1 y1 x2-x1 y2-y1]， EdgeColor ， red ， Linewidth ，1)；
轴(位置,[0.92,0.805,0.07,0.15]);
imagesc(时间，fre/1000，ABS(TFR))；
%xlabel(时间/秒);
ylabel( Fre/kHz )；
正常男性染色体组型轴
轴关闭；
%set(ha， xlim ，[x1 x2]， ylim ，[y1 y2])；
xlim([x1，x2])；ylim([y1，y2])；
支线剧情(512);
imagesc(time，fre/1000，ABS(Ts1))；
%xlabel(时间/秒);
ylabel( Fre/kHz )；
正常男性染色体组型轴；轴([0 0.1 y1-DD y2 DD])；
矩形(位置，[x1 y1 x2-x1 y2-y1]， EdgeColor ， red ， Linewidth ，1)；
轴(位置,[0.92,0.61,0.07,0.15]);
imagesc(time，fre/1000，ABS(Ts1))；
%xlabel(时间/秒);
ylabel( Fre/kHz )；
正常男性染色体组型轴
轴关闭；
%set(ha， xlim ，[x1 x2]， ylim ，[y1 y2])；
xlim([x1，x2])；ylim([y1，y2])；
支线剧情(513);
imagesc(time，fre/1000，ABS(Ts2))；
%xlabel(时间/秒);
ylabel( Fre/kHz )；
正常男性染色体组型轴；轴([0 0.1 y1-DD y2 DD])；
矩形(位置，[x1 y1 x2-x1 y2-y1]， EdgeColor ， red ， Linewidth ，1)；
轴(位置,[0.92,0.44,0.07,0.15]);
imagesc(time，fre/1000，ABS(Ts2))；
%xlabel(时间/秒);
ylabel( Fre/kHz )；
正常男性染色体组型轴
轴关闭；
%set(ha， xlim ，[x1 x2]， ylim ，[y1 y2])；
xlim([x1，x2])；ylim([y1，y2])；
支线剧情(514);
imagesc(time，fre/1000，ABS(Ts4))；
%xlabel(时间/秒);
ylabel( Fre/kHz )；
矩形(位置，[x1 y1 x2-x1 y2-y1]， EdgeColor ， red ， Linewidth ，1)；
正常男性染色体组型轴；轴([0 0.1 y1-DD y2 DD])；
轴(位置,[0.92,0.26,0.07,0.15]);
imagesc(time，fre/1000，ABS(Ts4))；
%xlabel(时间/秒);
ylabel( Fre/kHz )；
轴关闭；
%set(ha， xlim ，[x1 x2]， ylim ，[y1 y2])；
xlim([x1，x2])；ylim([y1，y2])；
支线剧情(515);
imagesc(time，fre/1000，ABS(MTs))；
%xlabel(时间/秒);
ylabel( Fre/kHz )；
正常男性染色体组型轴；轴([0 0.1 y1-DD y2 DD])；
矩形(位置，[x1 y1 x2-x1 y2-y1]， EdgeColor ， red ， Linewidth ，1)；
轴(位置,[0.92,0.09,0.07,0.15]);
imagesc(time，fre/1000，ABS(MTs))；
%xlabel(时间/秒);
ylabel( Fre/kHz )；
正常男性染色体组型轴
轴关闭；
%set(ha， xlim ，[x1 x2]， ylim ，[y1 y2])；
xlim([x1，x2])；ylim([y1，y2])；
..........
MTs1=MTs
MTs1(1:310，)=0；
MTs1(360:end，)=0；
MTs2=MTs
MTs2(1:260，)=0；
MTs2(310:end，)=0；
[Cs2]=brevridge_mult(abs(MTs2)，(1:2)/13.9，2，1，5)；
[Cs1]=brevridge_mult(abs(MTs1)，(1:2)/13.9，2，1，5)；
IF1=[Cs1(1，1:1130) Cs1(2，1131:end)]；
IF2=[Cs2(2，1:176) Cs2(1，177:end)]；
Cs(1，)=IF1Cs(2，)=IF2
n=1200
ds=3；
对于k=1:2
对于j=1:n
MTs_sig(k，j)=sum(real(MTs(max(1，Cs(k，j)-ds):min(round(n/2)，Cs(k，j) ds)，j))；
目标
目标
数字
suptitle(图21 );
子情节(411)
plot(time，fre(IF1)/1000， b-)；坚持住；plot(time，fre(IF2)/1000， b-)；
轴([0 0.1 y1-DD y2 DD])；
ylabel( Fre/kHz )；
子情节(412)
plot(time，MTs_sig(1， b-)；
ylabel(安培/克);
轴([0 0.1-2.5 2.5]);
子情节(413)
plot(time，MTs_sig(2， b-)；
ylabel(安培/克);
轴([0 0.1-1.5 1.5]);
子情节(414)
情节(时间，数据， k-，线宽,1);坚持住；
plot(time，sum(MTs_sig)， r-)；
ylabel(安培/克);
轴([0 0.1-2.5 2.5]);
图例(测量信号，‘S1S2’)；
.......
数字
suptitle(图22 );
支线剧情(211)
plot(time，real(MTs_sig(1，))/max(real(MTs_sig(1，))， b-， Linewidth ，2)；坚持住；
plot(time，detrend(fre(Cs(1，)))/max(detrend(fre(Cs(1，)))/1.25， r-，线宽,1.5);
轴([0 0.1-1 1]);
xlabel( Time/s )；
子情节(212)
plot(time，real(MTs_sig(2，))/max(real(MTs_sig(2，))， b-， Linewidth ，2)；坚持住；
plot(time，detrend(fre(Cs(2，)))/max(detrend(fre(Cs(2，)))/1.25， r-，线宽,1.5);
轴([0 0.1-1 1]);
图例("单组元模式"、"检测条件");
xlabel( Time/s )；
3运行结果
四参考文献[1]王智慧。小波时频同步压缩变换方法研究及应用d .西安电子科技大学， 2020.
郑重声明：本文由网友发布，不代表盛行IT的观点，版权归原作者所有，仅为传播更多信息之目的，如有侵权请联系，我们将第一时间修改或删除，多谢。