刚接触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('预处理：去除趋势项');

xlabel('时间/s');

ylabel('振幅/m')

grid on;

%预处理：小波降噪

[c,l]=wavedec(y1,3,'sym8');

[thr,sorh,keepapp]=ddencmp('den','wv',y1);

y2=wdencmp('gbl',y1,'sym8',3,thr,sorh,keepapp);

figure(3)

plot(t,y2,'r');

hold on

plot(t,y1);

axis([0 1.2 -8 8]*10^-3);

title('预处理：小波降噪');

xlabel('时间/s');ylabel('振幅/m')

grid on;

%降噪后频域图

f =fs*(1:N/2)/N;

S2=abs(fft(y1))*2/N;

figure(4)

plot(f,S2(1:N/2),'r');

hold on

plot(f,S1(1:N/2));

title('频域波形');

ylabel('振幅/m/s^2');

xlabel('频率/Hz');

grid on;

%预处理：带通滤波

Wn=[0.3,0.4];

[b,a]= butter(4,Wn,'s');

y3= filter(b,a,y2);

figure(5);

plot(t,y3);

axis([0 1.6*10^-3 -8*10^-7 8*10^-7]);

title('预处理：带通滤波');

grid on;

%带通后频域图

y4=fft(y3,N);

mag2=abs(y4);

f=(0:N-1)*fs/N;

figure(6)

plot(f,mag2,'r');

hold on

plot(f,S1(1:N));

axis([0 5*10^5 0 3*10^-4]);

title('频域波形');

ylabel('振幅/m/s^2');

xlabel('频率/Hz');

grid on;

%带通滤波后的包络时域图

X = hilbert(y3);%计算信号的hilbert变换

Amp= abs(X); %计算信号的瞬时功率

figure(7);

plot(t,Amp);%画包络时域图

title('带通滤波后的包络时域波形');

axis([0 1.6*10^-3 0 9*10^-7]);

ylabel('振幅/m/s^2');

xlabel('时间/s');

grid on;

%带通滤波后的包络频域图

envelop_fft = abs(fft(Amp))*2/N;

figure(8);

plot(f,10*envelop_fft(1:N));%画包络频谱图

axis([0 5*10^5 0 0.3*10^-6]);

title('包络频域波形');

ylabel('振幅/m/s^2');

xlabel('频率/Hz');

grid on