刚接触MATLAB,想请教老师们,如何编写一个噪声和冲击的模拟信号,还有关于共振解调的程序想请专业人士帮看看哪里有问题
fs = 10^6;%采样频率
N=4096;%采样长度
t=(0:N-1)/fs;%产生时间序列
x=load('jiaozai.m');
subplot(211)
plot(t,x)
title('时域波形');
axis([0 1.2 -8 8]*10^-3);
ylabel('振幅/m/s^2');
xlabel('时间/s');
grid on;
%初始信号的fft变换
f =fs*(1:N/2)/N;
S1=abs(fft(x))*2/N;
subplot(212)
plot(f,S1(1:N/2));
title('频域波形');
ylabel('振幅/m/s^2');
xlabel('频率/Hz');
grid on;
%对信号的预处理
m=1;
t1=reshape(t,4096,1);
a=polyfit(t1,x,m);
y1=x-polyval(a,t1);
figure(2)
plot(t1,y1,'r');
hold on
plot(t1,x);
axis([0 1.2 -8 8]*10^-3);
title('预处理:去除趋势项'