AM调制解调系统仿真（含软硬结合部分）

【实例截图】

【文件目录】

01 AM调制信号分析

├── 01 Matlab仿真
│   ├── AM_Encode.asv
│   ├── AM_Encode.m
│   ├── main.asv
│   └── main.m
└── 02 软硬结合
    ├── AM_Encode.m
    ├── AM_ROM_.DATA_ch1.mif
    ├── AM_ROM_.DATA_ch2.mif
    ├── DA_OUT.m
    ├── SaveData.m
    └── main.m

2 directories, 10 files

【核心代码】

%%调制信号产生
[y1,y2,y3,y] = AM_Encode(A,F,Fc,t);
%%画图
%调制信号
figure(1)
plot(t,y1);title('调制信号');
xlabel('时间/(ms)');ylabel('幅度/V');
axis([49.96*10^-3 50.04*10^-3,-inf,inf]);
set(gca,'xticklabel',get(gca,'xtick')*1e3)

freq=fft(y1)*2/N;
w=(-N/2:1:N/2-1)*freqPixel; %双边
figure(9)
freq_d=abs(fftshift(freq));
plot(w,freq_d);
xlim([-2*Fc,2*Fc]);
title('调制信号频谱');
xlabel('频率/Hz');ylabel('幅度/V');
axis([-10^5 10^5,-inf,inf]);
%载波
figure(2)
plot(t,y2);title('载波');
xlabel('时间/0.1ms');ylabel('幅度/V');
axis([49.96*10^-4 50.04*10^-4,-inf,inf]);
set(gca,'xticklabel',get(gca,'xtick')*1e4)

freq=fft(y2)*2/N;
w=(-N/2:1:N/2-1)*freqPixel; %双边
figure(10)
freq_d=abs(fftshift(freq));
plot(w,freq_d);
xlim([-2*Fc,2*Fc]);
title('载波信号频谱');
xlabel('频率/Hz');ylabel('幅度/V');

%AM
figure(3)
plot(t,y);title('已调信号');
xlabel('时间/s');ylabel('幅度/V');
axis([49.96*10^-3 50.04*10^-3,-inf,inf]);
set(gca,'xticklabel',get(gca,'xtick')*1e3)

%加直ＡＭ
figure(8)
plot(t,y3);title('加直流分量的已调信号');
xlabel('时间/0.1ms');ylabel('幅度/V');
axis([49.96*10^-3 50.04*10^-3,-inf,inf]);
set(gca,'xticklabel',get(gca,'xtick')*1e3)

freq=fft(y3)*2/N;
w=(-N/2:1:N/2-1)*freqPixel; %双边
figure(11)
freq_d=abs(fftshift(freq));
plot(w,freq_d);
xlim([-2*Fc,2*Fc]);
title('加直流分量信号频谱');
xlabel('频率/Hz');ylabel('幅度/V');
axis([-10^5 10^5,-inf,inf]);


%加噪
y4=awgn(y,0.001);
figure(4)
plot(t,y4);title('已调信号加噪声');
xlabel('时间/s');ylabel('幅度/V');
axis([49.96*10^-3 50.04*10^-3,-inf,inf]);
set(gca,'xticklabel',get(gca,'xtick')*1e3)

freq=fft(y4)*2/N;
w=(-N/2:1:N/2-1)*freqPixel; %双边
figure(12)
freq_d=abs(fftshift(freq));
plot(w,freq_d);
xlim([-2*Fc,2*Fc]);
title('加噪信号频谱');
xlabel('频率/Hz');ylabel('幅度/V');

%求信号频谱
freq=fft(y)*2/N;
w=(-N/2:1:N/2-1)*freqPixel; %双边
figure(5)
freq_d=abs(fftshift(freq));
plot(w,freq_d);
xlim([-2*Fc,2*Fc]);
title('已调信号频谱');
xlabel('频率/Hz');ylabel('幅度/V');
%功率谱
freq_s(1:length(freq_d)/2)=freq_d(length(freq_d)/2 1:length(freq_d));
ypsd=freq_s.*conj(freq_s);
w2=(0:length(freq_s)-1)*fs/length(freq_d); %单边
figure(6)
plot(w2,ypsd);
xlim([0,2*Fc]);
title('已调信号功率谱');
xlabel('频率/Hz');
ylabel('W/Hz');
axis([-10^6 10^6,-inf,inf]);
%自相关函数
figure(7)
[Rx,maxlags]=xcorr(y);  %信号的自相关
plot(maxlags/fs,Rx);
title('信号自相关');
xlabel('时间/ms');
ylabel('幅值/V');
axis([-4*10^-5 4*10^-5,-inf,inf]);
set(gca,'xticklabel',get(gca,'xtick')*1e3)
%均值、方差
aver=mean(y2);
v=var(y2);