博主大毕设关于数字下变频（DDC）的CUDA实现，预计工期比较长，所以留下一些文字记录。
主要分为两部分工作，Matlab仿真部分和CUDA实现。
由于很久没有仿真了，所以先用一个简单的AM调制仿真练手。具体代码如下：
代码都是基础的所以就不解释了(环境matlab2016)
1 clc; 2 fm=100; %信号频率 3 fc=500; %载波频率 4 fs=5000; %抽样频率 5 Am=1; 6 A=2; 7 N=512; 8 K=N-1; 9 n=0:N-1; 10 t=(0:1/fs:K/fs); 11 yt=Am*cos(2*pi*fm*t); 12 figure(1) 13 subplot(1,1,1);plot(t,yt);title('频率为100HZ信号时域波形'); 14 15 y0=A+yt; 16 y2=y0.*cos(2*pi*fc*n/fs); 17 y3=fft(y2,N); 18 q1=(0:N/2-1)*fs/N; 19 mx1=abs(y3(1:N/2)); 20 figure(2) 21 subplot(2,1,1); 22 plot(t,y2); 23 title('已调信号时域波形'); 24 subplot(2,1,2); 25 plot(q1,mx1); 26 title('已调信号频谱'); 27 yc=cos(2*pi*fc*t); 28 figure(3); 29 subplot(2,1,1),plot(t,yc),title('载波fc时域') 30 n=0:N-1; 31 yc1 = Am*cos(2*pi*fc*n/fs); 32 y3=fft(yc1,N) 33 q=(0:N/2-1)*fs/N 34 mx=abs(y3(1:N/2)); 35 figure(3) 36 subplot(2,1,2),plot(q,mx),title('载波fc频谱') 37 N=512; 38 n=0:N-1; 39 y4=0.01*randn(1,length(t)); %产生高斯噪声 40 w=y4.^2; %噪声功率 41 figure(4); 42 subplot(2,1,1); 43 plot(t,y4); 44 title('高斯白噪声时域波形'); 45 y5=fft(y4,N); 46 q2=(0:N/2-1)*fs/N; 47 mx2=abs(y5(1:N/2)); 48 figure(4); 49 subplot(2,1,2); 50 plot(q2,mx2); 51 title('高斯白噪声频域波形'); 52 y6=y2+y4; 53 figure(5); 54 subplot(2,1,1); 55 plot(t,y6); 56 title('加噪声后时域信号') 57 q3=q1; 58 mx3=mx1+mx2; 59 subplot(2,1,2); 60 plot(q3,mx3); 61 title('加噪声后频谱') 62 63 64 yv=y6.*yc; %乘以载波想干解调 65 Ws=yv.^2; 66 p1=fc-fm; 67 [k,Wn,beta,ftype]=kaiserord([p1 fc],[1 0],[0.05 0.01],fs);%数字低通过滤波器 68 window=kaiser(k+1,beta);%使用kaiser窗函数 69 b=fir1(k,Wn,ftype,window,'noscale');% 70 yt=filter(b,1,yv); 71 yssdb=yt.*2-2; 72 figure(6) 73 subplot(2,1,1); 74 plot(t,yssdb); 75 title('经过低通已调信号时域波形') ;%解调 76 77 y9=fft(yssdb,N); 78 mx=abs(y9(1:N/2)); 79 subplot(2,1,2); 80 plot(q,mx); 81 title('已调信号频域波形')
实现效果如下：
