tpsk.m

Agilent N5182A信号源自己产生psk fsk ofdm等信号的matlab程序 产生的文件在Agilent N5182A上测试成功

M = 8; % 　进制数
fd=8000;%码元速率    
fs=16000;%采样频率
rolloff=0.5;%滚降系数
% X = randsrc(1,1000,[0 1]); 
% X=randsrc(1,1000,[0 1 2 3 ]);  
X=randsrc(1,1000,[0 1 2 3 4 5 6 7]);  
Y = PSKMOD(X,M);
rrcfilter=rcosine(fd,fs,'fir/sqrt',rolloff);
I_filtered=rcosflt(real(Y),fd,fs,'filter',rrcfilter);
Q_filtered=rcosflt(imag(Y),fd,fs,'filter',rrcfilter);
[num,z]=size(I_filtered);
filtered=I_filtered+i*Q_filtered;
r=fft(filtered)/(num);
r=fftshift(r); 
f=linspace(-num/2,num/2,num);
plot(f,20*log10(abs(r))); 
axis([-num,num,-100,20]);
waveform=change(I_filtered,Q_filtered,num);
filename2 = 'C:\Temp\psksqrt';
[FID, message] = fopen(filename2,'w');   % Open a file to write data
if FID == -1 error('Cannot Open File'); end
fwrite(FID,waveform,'unsigned short');  % write to the file
fclose(FID);