MATALA杂波代码

clear all;close all;clc;
p=1.5q=0.5;
fr=1000;
num=8000;%采样点数数据点数
n=0:1:num-1;
%-----由功率谱和相关系数确定相关高斯序列；
fs=1000;%采样频率
f=fs/num*（-num/2:num/2-1）;%对称的功率谱
sigmaf=50;
amp=1/（sqrt（2*pi）*sigmaf）;
Pw=amp*exp（-（f）.^2/（2*sigmaf^2））;%高斯功率谱


sn=ifft（Pw）;%由最终的高斯功率分布得到相关系数序列
r=sqrt（sn*（gamma（1+2/p）-gamma（1+1/p）.^2）/gamma（1+1/p）.^2）;%这里这个式子是怎么得出的这个很是不理解 怎么解出的？并且原式中应该还应该与韦布尔分布的参数有关系这里完全是两个相关系数的关系不知道为什么原式中还有高斯超几何函数怎能么就没了呢 从两个2一个3来看应该是默认的p=1

Xw=fft（rnum）;     %相关高斯随机序列功率谱
abs_Hw=sqrt（abs（Xw）/r（1））;%这个有一个流程图可以表示

w=rand（1num）;%生成随机序列但是均匀分布的
ww=2*pi*w;
ww=ww/max（abs（ww））;%应该是归一化用
nois=exp（j*（ww））;  %加入的随机相位 相干信号  这里确实在一篇论文中有提起 但是确实还不是很明白
Hw=abs_Hw.*nois; %成形滤波器频谱
%%%%%%%%%%%%%%%
vi1=randn（1num）;
vq1=randn（1num）;
v1=vi1+j*vq1;
%std（v1）
subplot（221）
plot（nabs（v1））
title（‘独立不相关高斯随机序列‘）
xlabel（‘a‘）
V1=fft（v1）;
W1=Hw.*V1;
w1=ifft（W1）;
w1sigma=std（w1）;
w1mu=mean（w1）;
w1=（w1-w1mu）/w1sigma;
%--------------
vi2=randn（1num）;%为什么原来是rands 这里不明白
vq2=randn（1num）;
v2=vi2+j*vq2;
V2=fft（v2）;
W2=Hw.*V2;
w2=ifft（W2）;
w2sigma=std（w2）;
w2mu=mean（w2）;
w2=（w2-w2mu）/w2sigma;

%-------------
w1=w1*