%?Project1

%?日期:2006.12.28

close?all;clear;clc
%Variable?declare
N=10000;????????????????????????????????????????????????????????????????????%?二進制碼流比特數
T=0.001;????????????????????????????????????????????????????????????????????%?16QAM符號長度（s）
Fd=1/T;?????????????????????????????????????????????????????????????????????%?符號率
Fs=4000;????????????????????????????????????????????????????????????????????%?過采樣率（Hz）
r=Fs*T;?????????????????????????????????????????????????????????????????????%?過采樣因子
alpha=0.25;?????????????????????????????????????????????????????????????????%?滾降因子
Fc=500;?????????????????????????????????????????????????????????????????????%?信道截止頻率（Hz）
SNR=15;?????????????????????????????????????????????????????????????????????%?信道的信噪比

%Bit?Source
source=randint（N1）;????????????????????????????????????????????????????????%?產生二進制碼流

%QAM?Coder
send=bi2de（flipud（reshape（source4N/4））‘）;
QAM=qammod（send16）;
h1=scatterplot（QAM）;
set（h1‘Name‘‘Constellation?Diagram‘）;
title（‘Constellation?of?16QAM?in?source‘）;

%Transmitter
p=rcosine（FdFs‘sqrt‘alpha）;??????????????????????????????????????????????%?均方根升余弦濾波器

len=length（QAM）;
QAMS=zeros（size（1:r*len））;
QAMS（1:r:r*len）=QAM;
QAM_x=real（QAMS）;
QAM_y=imag（QAMS）;
delay=（length（p）-1）/2;
temp1=zeros（size（1:delay））;
T_x_a=[QAM_x?temp1];
T_y_a=[QAM_y?temp1];
T_x_a=filter（p1T_x_a）;
T_y_a=filter（p1T_y_a）;
T_x=T_x_a（1+delay:10000+delay）;
T_y=T_y_a（1+delay:10000+delay）;
figure（2）;
hold?on;
for?i=（length（p）+1）/2:r*2:len-r*2
????plot（T_x（i:i+r*2））;
end;
hold?off;
figure（3）;
freqz（T_x1256）;

%Channal
%[B?A]=butter（22*Fc/Fs）;
B=fir1（82*Fc/Fs）;A=1;
delay=（length（B）-1）/2;
temp1=zeros（size（1:delay））;
C_x_a=[T_x?temp1];
C_y_a=[T_y?temp1];
C_x_a=filter（BAC_x_a）;
C_y_a=filter（BAC_y_a）;
C_x=C_x_a（1+delay:10000+delay）;
C_y=C_y_a（1+delay:10000+delay）;
C_x=awgn（C_xSNR）;
C_y=awgn（C_ySNR）;
figure（4）;
freqz（C_x1256）;

%Equalizer
load?c;
freqz（BA256）;
hold?on;
freqz（c1256）;
hold?off;
delay=（length（c）-1）/2;
temp1=zeros（size（1:delay））;
E_x_a=[C_x?temp1];
E_y_a=[C_y?temp1];
E_x_a=filter（c1E_x_a）;
E_y_a=filter（c1E_y_a）;
E_x=E_x_a（1+delay:10000+delay）;
E_y=E_y_a（1+delay:10000+delay）;

%E_x=C_x;
%E_y=C_y;

%Receiver
delay=（length（p）-1）/2;
temp1=zeros（size（1:delay））;
R_x_a=[E_x?temp1];
R_y_a=[E_y?temp1];
R_x_a=filter（p1R_x_a）;
R_y_a=filter（p1R_y_a）;
R_x=R_x_a（1+delay:10000+delay）;
R_y=R_y_a（1+delay:10000+delay）;
figure（5）;
hold?on;
for?i=（length（p）+1）/2:r*2:len-r*2
????plot（R_x（i:i+r*2））;
end;
hold?off;

%Sample
DQAM_x=R_x（1:4:10000）;
DQAM_y=R_y（1:4:10000）;

Error_x=DQAM_x‘-real（QAM）;
Error_y=DQAM_y‘-imag（QAM）;

DQAM=DQAM_x+j*DQAM_y;
h2=scatterplot（DQAM）;

%Decoder
rcv?=?qamdemod（DQAM16）;

%Error
[Err_numErr_ratio]=symerr（rcvsend‘）

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件?????141824??2007-01-17?19:18??通信原理大作業（1）.doc

?????文件???????2908??2007-12-11?16:47??project1.m

-----------?---------??----------?-----??----

???????????????144950????????????????????3