文档库

最新最全的文档下载
当前位置:文档库 > 西安邮电学院通信原理实验结果5

西安邮电学院通信原理实验结果5

function [m]=F2T(M,fs)

%-------------------------输入参数

%M:信号的频谱

%fs:系统采样频率

%--------------------输出(返回)参数

%m:傅里叶逆变换后的信号,注意其长度为2的整数次幂,利用其画波形时,要注意选取m的一部分,选取长度和所给时间序列t的长度要一致,plot(t,m(1:length(t))),否则会出错。

m = real(ifft(M))*fs;

function p=power_x(x)

%x:输入信号

%p:返回信号的x功率

p=(norm(x).^2)./length(x);

function [M,m,df]=fftseq(m,ts,df)

%各参数含义与子函数T2F中的完全相同,完成

fs = 1/ts;

if nargin ==2

n1 =0;

else

n1 = fs/df;

end

n2 = length(m);

n = 2^(max(nextpow2(n1),nextpow2(n2)));

M = fft(m,n);

m = [m,zeros(1,n-n2)];

df = fs/n;

function [M,m,df]=fftseq(m,ts,df)

%各参数含义与子函数T2F中的完全相同,完成

fs = 1/ts;

if nargin ==2

n1 =0;

else

n1 = fs/df;

end

n2 = length(m);

n = 2^(max(nextpow2(n1),nextpow2(n2)));

M = fft(m,n);

m = [m,zeros(1,n-n2)];

df = fs/n;

function [M,m,df1,f]=T2F(m,ts,df,fs)

%------------------------输入参数

%m:信号

%ts:系统时域采样间隔

%df:所需的频率分辨率

%fs:系统采样频率

%---------------------输出(返回)参数

%M:傅里叶变换后的频谱序列