实在弄不出来了，求教各位MATLAB大神

fsample = 9000;
fsymbol = 1500;
sy = fsample/fsymbol;
Number = 20;
walsh0 = (1:sy)*0+1;
Data_I = sign(randn(Number*sy,1));
Data_I1 = kron(Data_I,walsh0);
Data_I2 = reshape(Data_I',[],1);
rcoscoeff = rcosine(fsymbol,fsample,'sqrt',0.35,5);
Data_I3 = conv(rcoscoeff,Data_I2);
Data_I3 = Data_I3(31:end-30);
fsin = sin(2*pi*(1:Number*sy)*1800/1900);
fqpsk = Data_I3.*fsin';
noise = randn(length(fqpsk),1);
fqpsknoise = fqpsk+noise;
fqpsknorcos = Data_I2.*fsin';
figure(1)    
subplot(3,1,1);
plot(Data_I2)
axis([0 Number*sy 1.1*min(Data_I2) 1.1*max(Data_I2)])
ylabel('amplitude');xlabel('(a)point');
subplot(3,1,2);
plot(Data_I3)
ylabel('amplitude');xlabel('(b)point');
subplot(3,1,3);
plot(rcoscoeff);
ylabel('amplitude');xlabel('(c)point');
figure(2)  
freqz(rcoscoeff);
figure(3)
subplot(3,1,1);
plot(fqpsk);
ylabel('amplitude');xlabel('(a)point');
subplot(3,1,2);plot(fqpsknoise);
ylabel('amplitude');xlabel('(b)point');
subplot(3,1,3);
plot(fqpsknorcos);
ylabel('amplitude');xlabel('(c)point');
figure(4)
n = (1:Number*sy)/(Number*sy)*2;
subplot(3,1,1);
plot(n,20*log(abs(fft(fqpsk))))
ylabel('magnitude(dB)');xlabel('(a)Normalized Frequency(\time\pi rad\sample)');
subplot(3,1,2);
plot(n,20*log(abs(fft(fqpsknoise))))
ylabel('magnitude(dB)');xlabel('(b)Normalized Frequency(\time\pi rad\sample)');
subplot(3,1,3);
plot(n,20*log(abs(fft(fqpsknorcos))))
ylabel('magnitude(dB)');xlabel('(c)Normalized Frequency(\time\pi rad\sample)');