BPSK QPSK 16QAM 64QAM的matlab仿真程序

编写了BPSK QPSK 16QAM等信号的MATLAB代码，运行得出波形。

function y=ask_2（snr_in）;
snr=10^（snr_in/10）;
x=rand（1001）;%产生随机序列
h=length（x）;%计算序列长度
fc=10000;fs=100000;%设定采样频率以及高频和低频（调制频率） 
ts=10/fc;%设定码元宽度;每码元100个点 
t=0:1/fs:ts*h-1/fs;%产生全部码元总的时间序列起始时刻为0步长为采样频率的倒数即采样时间结束时刻为序列结束时刻
l=length（t）/h;%计算出这个序列的长度
t=reshape（t‘lh）;%将序列处理成矩阵形式，每行元素为一个码元内的所有采样点l，行数为码元数n  ：reshape（lh）
t=t‘;%因矩阵是按列存储，故须作转置处理
y=zeros（hl）;
for i=1:h
    if x（i）>=0.5
         y（i:）=3.5*sin（2*pi*fc.*t（i:））;
    else y（i:）=sin（2*pi*fc.*t（i:））;%赋值
    end
end
y0=reshape（y‘l*h1）;
power=y0‘*y0;
powerzao=power/snr;%以下是产生高斯白噪声序列
a=sqrt（powerzao/（l*h））;
zao=a*randn（l*h1）;
powerzao2=zao‘*zao;
y=y0+zao;%加噪以后的信号
t=reshape（t‘l*h1）;
subplot（311）;plot（ty）;
grid on
title（‘2ASK信号的时域图形‘）

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件        904  2008-12-20 16:02  ask_2.m

     文件        862  2008-11-05 16:13  ask_4.m

     文件        957  2008-12-20 16:03  fsk_2.m

     文件        937  2008-12-20 16:03  fsk_4.m

     文件        607  2008-12-22 09:56  main.m

     文件        683  2008-11-13 14:00  psk_2.m

     文件       2154  2008-11-03 16:57  QAM_16.m

     文件        867  2008-12-03 15:53  qpsk.m

     文件        779  2008-12-20 15:54  TX.m

     文件        865  2008-12-04 10:05  TZSB.m

----------- ---------  ---------- -----  ----

                 9615                    10