主题:[讨论]有谁帮忙看下程序怎么改?
clear;
%====定义待仿真序列的维数 N
global N
N=320;
%====定义产生‘1’的概率为 p
global p
p=0.5;
%==============================
%首先产生随机二进制序列
source=randsrc(1,N,[1,0;p,1-p]);
%==============================
%对产生的二进制序列进行QAM调制
[source1,source2]=Qam_modulation(source);
%===============================
%画出星座图
figure(1);
plot_astrology(source1,source2);
%===============================
%两路信号进行插值
sig_insert1=insert_value(source1,8);
sig_insert2=insert_value(source2,8);
%===============================
%画出两路信号的波形图
figure(2);
plot_2way(sig_insert1,sig_insert2,length(sig_insert1),0.5);
title('两路信号的波形图');
%===============================
%通过低通滤波器
[sig_rcos1,sig_rcos2]=rise_cos(sig_insert1,sig_insert2,0.25,2);
%===============================
%画出两路信号的波形图
figure(3);
plot_2way(sig_rcos1,sig_rcos2,length(sig_rcos1)/4,0.5);
title('通过低通滤波器后两路信号波形图');
%stem_2way(sig_insert1,sig_insert2,length(sig_insert1)/4,0.5);
%===============================
%====将基带信号调制到高频上
[t,sig_modulate]=modulate_to_high(sig_rcos1,sig_rcos2,0.25,2.5);
figure(4);
plot(t(1:500),sig_modulate(1:500));
%===============================
%====将滤波后的信号加入高斯白噪声
snr=10;
[x1,x2]=generate_noise(sig_rcos1,sig_rcos2,snr);
sig_noise1=x1';
sig_noise2=x2';
end;
figure(5)
plot_2way(sig_noise1,sig_noise2,length(sig_noise1)/4,0.5);
%===============================
%====经过匹配滤波器
% [x1,x2]=match_flt(sig_noise1,sig_noise2,0.25,2);
% sig_match1=x1';
% sig_match2=x2';
[sig_match1,sig_match2]=rise_cos(sig_noise1,sig_noise2,0.25,2);
figure(6);
plot_2way(sig_match1,sig_match2,length(sig_match1)/4,0.5);
%===============================
%采样
[x1,x2]=pick_sig(sig_match1,sig_match2,8);
sig_pick1=x1;
sig_pick2=x2;
%画出星座图
figure(7)
plot_astrology(sig_pick1,sig_pick2);
%解调
signal=demodulate_sig(sig_pick1,sig_pick2);
%画出误码率曲线图
figure(8)
plot_snr;
哪个帮忙看下?
%====定义待仿真序列的维数 N
global N
N=320;
%====定义产生‘1’的概率为 p
global p
p=0.5;
%==============================
%首先产生随机二进制序列
source=randsrc(1,N,[1,0;p,1-p]);
%==============================
%对产生的二进制序列进行QAM调制
[source1,source2]=Qam_modulation(source);
%===============================
%画出星座图
figure(1);
plot_astrology(source1,source2);
%===============================
%两路信号进行插值
sig_insert1=insert_value(source1,8);
sig_insert2=insert_value(source2,8);
%===============================
%画出两路信号的波形图
figure(2);
plot_2way(sig_insert1,sig_insert2,length(sig_insert1),0.5);
title('两路信号的波形图');
%===============================
%通过低通滤波器
[sig_rcos1,sig_rcos2]=rise_cos(sig_insert1,sig_insert2,0.25,2);
%===============================
%画出两路信号的波形图
figure(3);
plot_2way(sig_rcos1,sig_rcos2,length(sig_rcos1)/4,0.5);
title('通过低通滤波器后两路信号波形图');
%stem_2way(sig_insert1,sig_insert2,length(sig_insert1)/4,0.5);
%===============================
%====将基带信号调制到高频上
[t,sig_modulate]=modulate_to_high(sig_rcos1,sig_rcos2,0.25,2.5);
figure(4);
plot(t(1:500),sig_modulate(1:500));
%===============================
%====将滤波后的信号加入高斯白噪声
snr=10;
[x1,x2]=generate_noise(sig_rcos1,sig_rcos2,snr);
sig_noise1=x1';
sig_noise2=x2';
end;
figure(5)
plot_2way(sig_noise1,sig_noise2,length(sig_noise1)/4,0.5);
%===============================
%====经过匹配滤波器
% [x1,x2]=match_flt(sig_noise1,sig_noise2,0.25,2);
% sig_match1=x1';
% sig_match2=x2';
[sig_match1,sig_match2]=rise_cos(sig_noise1,sig_noise2,0.25,2);
figure(6);
plot_2way(sig_match1,sig_match2,length(sig_match1)/4,0.5);
%===============================
%采样
[x1,x2]=pick_sig(sig_match1,sig_match2,8);
sig_pick1=x1;
sig_pick2=x2;
%画出星座图
figure(7)
plot_astrology(sig_pick1,sig_pick2);
%解调
signal=demodulate_sig(sig_pick1,sig_pick2);
%画出误码率曲线图
figure(8)
plot_snr;
哪个帮忙看下?