xy506_sess6and7.m

来自「一个利用matlab simulink仿真简化通信系统的模型。里面m文件为sim」· M 代码 · 共 26 行

global Fs audio
global M stddev partition codebook predictor_num;
global constellation eq
global insig outsig;

load TEST_AUD.MAT

chan=[1 (10^(-5/20))*exp(j*1000000*1/Fs) (10^(-8/20))*exp(j*1000000*2/Fs)];

chan_m=[chan(1) chan(2) chan(3) 0 0; 0 chan(1) chan(2) chan(3) 0; 0 0 chan(1) chan(2) chan(3); 0 0 0 chan(1) chan(2); 0 0 0 0 chan(1)];
eq=[1 0 0 0 0]*(inv(chan_m));

M=8; % the number of constellation point

constellation=[1 0.707*(1+j) 0.707*(-1+j) j 0.707*(1-j) -j -1 0.707*(-1-j)];

signal_power=1; % From the constellation diagram i.e.modulation scheme
snr=12;
stddev=sqrt(signal_power/(10^(snr/10)));

bits=log2(M); % 2^bits=M;
[predictor_num, codebook, partition]= dpcmopt(audio, 1, (2^bits));

sim('xy506_Ses6and7');

sound(outsig,Fs);