%E4_8_NoiseAndCarrier.M
f4=100;
f1=300;???????%信號(hào)1頻率為200Hz
f3=400;
f2=500;???????%信號(hào)2頻率為800Hz
Fs=2000;??????%采樣頻率為2KHz
N=12;?????????%量化位數(shù)
%Pn=Pn0*Fs/2;?%噪聲信號(hào)功率
%width=0.5;???%函數(shù)SAWTOOTH（）的尺度參數(shù)為0.5
%duty=50;?????%函數(shù)SQUQRE（）的尺度參數(shù)為50
%產(chǎn)生信號(hào)
t=0:1/Fs:1;
c1=2*pi*f1*t;
c2=2*pi*f2*t;
c3=2*pi*f3*t;
c4=2*pi*f4*t;
%sq=square（cduty）;%產(chǎn)生方波
%tr=sawtooth（cwidth）;%產(chǎn)生三角波
s1=sin（c1）;%產(chǎn)生正弦波
s2=sin（c2）;%產(chǎn)生正弦波
s3=sin（c3）;
s4=sin（c4）;
s=s1+s2+s3+s4;???%產(chǎn)生兩個(gè)單載波合成后的信號(hào)
%產(chǎn)生隨機(jī)序列信號(hào)
noise=randn（1length（t））;%產(chǎn)生高斯白噪聲序列

%歸一化處理
noise=noise/max（abs（noise））;
s=s/max（abs（s））;

%12比特量化
Q_noise=round（noise*（2^（N-1）-1））;
Q_s=round（s*（2^（N-1）-1））;

%調(diào)用自已設(shè)計(jì)的濾波器函數(shù)對(duì)信號(hào)進(jìn)行濾波
hn=E4_8_Fir8Serial;
Filter_noise=filter（hn1Q_noise）;
Filter_s=filter（hn1Q_s）;

%求信號(hào)的幅頻響應(yīng)
m_noise=20*log（abs（fft（Q_noise1024）））/log（10）;?m_noise=m_noise-max（m_noise）;
m_s=20*log（abs（fft（Q_s1024）））/log（10）;?m_s=m_s-max（m_s）;
%濾波后的幅頻響應(yīng)
Fm_noise=20*log（abs（fft（Filter_noise1024）））/log（10）;?Fm_noise=Fm_noise-max（Fm_noise）;
Fm_s=20*log（abs（fft（Filter_s1024）））/log（10）;?Fm_s=Fm_s-max（Fm_s）;
%濾波器本身的幅頻響應(yīng)
m_hn=20*log（abs（fft（hn1024）））/log（10）;?m_hn=m_hn-max（m_hn）;

%設(shè)置幅頻響應(yīng)的橫坐標(biāo)單位為Hz
x_f=[0:（Fs/length（m_s））:Fs/2];
%只顯示正頻率部分的幅頻響應(yīng)
mf_noise=m_noise（1:length（x_f））;
mf_s=m_s（1:length（x_f））;
Fmf_noise=Fm_noise（1:length（x_f））;
Fmf_s=Fm_s（1:length（x_f））;
Fm_hn=m_hn（1:length（x_f））;

%繪制幅頻響應(yīng)曲線
subplot（211）
plot（x_fmf_noise‘-.‘x_fFmf_noise