Signal_No?=?1;
Interference_No?=?2;
S_No?=?Signal_No?+?Interference_No;
sensor_No?=?8;
azimuth?=?[0*pi/180?10*pi/180?70*pi/180];
Fs?=?30e6;
Fc?=?10e6;
F1?=?[8e6?20e6];
theta?=?[10?45?52?93?43]*pi/180;
RB?=?2e6;
M?=?2;
wavelength?=?3e8/Fc;??%波長
d?=?wavelength/2;
K?=?30;
Ps?=??-30;
SNR?=?10;
SIR?=?-20;
M_No?=?40;
Data_No?=?M_No*Fs/RB;
t?=?1/Fs:1/Fs:Data_No/Fs;
Ps_1?=?10^（Ps/10）;
Pi_1?=?10^（（Ps-SIR）/10）;
Pn_1?=?10^（（Ps-SNR）/10）;
bit?=?randint（1M_No）;
bitstream?=?[];
?for?i?=?1:M_No
???????if?bit（i）==1
?????????????bitstream?=?[bitstream?ones（1?Fs/RB）];
???????else
?????????????bitstream?=?[bitstream?-ones（1?Fs/RB）];
???????end
end
?Carrier_R?=?cos（2*pi*Fc*t）;
?S_R?=?Carrier_R?.*?bitstream;
?Carrier_I?=?sin（2*pi*Fc*t）;
?S_I?=?Carrier_I?.*bitstream;
?
?b?=?sqrt（（S_R*S_R‘/length（S_R）））;
?
?Signal_R?=?sqrt（Ps_1）*S_R/sqrt（（S_R*S_R‘/length（S_R）））;
?Signal_I?=?sqrt（Ps_1）*S_I/sqrt（（S_I*S_I‘/length（S_I）））;
?Signal（1:）?=?comp