function?Signal=QAM（RsFsFc）
%?this?file?used?to?produce?qam?signal
%?clear?all;

%?simulation?time
Temulate?=?3;
%roll?factor
alfa?=?1;
%carrier?intial?phase
Qi?=?pi/8;
%snr
nos?=?21;

%?%?symbol?rate
%?Rs?=?9600;
%?%?sample?rate
%?Fs=?Rs*8;
%?%?carrier?frequency
%?Fc=?19200;
%?%?Fc=Rs*2;

string?=?‘signal_16qamdata.dat‘;
bit_number?=?floor（Temulate*Rs）;

rawdata?=?rand（1bit_number）;
for?i=1:bit_number
???if（rawdata（i）>=0.5）
??????rawdata（i）=1;
???else
??????rawdata（i）=0;
???end
end

h?=?rcosfir（alfa[-33]811）;%;%1rcosfir（alfa1121‘sqrt‘）

In?=?zeros（1bit_number*8/4+8）;
Qn?=?zeros（1bit_number*8/4+8）;

j?=?1;
for?i?=?1:4:bit_number
????if（（rawdata（i）?==?0）?&?（rawdata（i+1）?==?0））
????????In（1+j*8）?=?3;%*0.41;
????elseif（（rawdata（i）?==?0）?&?（rawdata（i+1）?==?1））
????????In（1+j*8）?=?1;%0.41;
????elseif（（rawdata（i）?==?1）?&?（rawdata（i+1）?==?0））
????????In（1+j*8）?=?-1;%0.41;
????elseif（（rawdata（i）?==?