DDS 期末作业总结

2019-11-08 00:47:58

字体：大中小

来源：转载

供稿：网友

姓名	绘制RTL图	工具生成RTL图	频率字推导	频率字准确	相位累加器波形截图	波表ROM地址截图	正弦波形截图	Matlab频谱截图
田栋	有	有	有	否	无	无	有	有、频率轴不正确
殷松	有	有	无	是	无	有	有	无
王婷琦	有	有	有	是	有	有	有	有、频率轴不正确
李京会	有	有	无	否	无	有	有	有、频率轴不正确、幅度不正常
鞠汉	有	有	无	是，但误差分析不对	看不清	看不清	有	有、频率轴不正确
刘权	有	有	无	是	无	有	有	无
刘溪	有	有	无	否	无	有	有	有、频率轴不正确
马德壮	有	有	有	是	有	被省略	有	有、频率轴不正确
夏健力	有	有	有	否	无	看不清	有	有、频率轴不正确
郭云海	有	有	有、不正确	否	无	看不清	有	有、频率轴不正确
刘阳	有	有	无	是	有	被省略	有	有、频率轴不正确
郭治卫	有	有	无	否	无	无	有	有、频率轴不正确

重复情况统计

相位累加器截图有出现重复

截图时间不同，但是信号名称及图样非常相似，差异性很小

重复波形图-1

重复波形图-2

Matlab频谱分析代码的重复

总结：频谱分析代码的设计思路非常一致，而且关于频率轴的错误出现的错误图样也完全一致。教训：以后在期末大作业提交之前，不能公开各个同学的博客地址

CODE 1

clear;sin_wave_out;signal = reshape(sin_wave, 1024, 1); %矩阵转置fs=5e6; %采样频率N=1024; %采样点数t=0:1/fs:(N-1)/fs; %采样时刻figure(1);plot(t,signal);xlabel('Time/s'); ylabel('Magnitude');Y = fft(signal,N); %做FFT变换Ayy = abs(Y); %取模Ayy=Ayy/(N/2); %换算成实际的幅度F=(0:N-1)*fs/N; %换算成实际的频率值，Fn=(n-1)*Fs/Nfigure(2);stem(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果axis([0 5E6 0 2500]);title('幅度-频率曲线图'); xlabel('Frequency (Hz)'); ylabel('Magnitude');

CODE 2

clear;signal_tap_data_3M;signal=transpose(signal_out); %矩阵转置fs=50E6; %采样频率N=1024; %采样点数t=[0:1/fs:(N-1)/fs]; %采样时刻figure(1);plot(t,signal);xlabel('Time (s)');ylabel('Magnitude');Y = fft(signal,N); %做FFT变换Ayy = abs(Y); %取模Ayy=Ayy/(N/2); %换算成实际的幅度Ayy(1)=Ayy(1)/2;F=([1:N]-1)*fs/N; %换算成实际的频率值，Fn=(n-1)*Fs/Nfigure(2);stem(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果axis([0 5E6 0 2500]);title('幅度-频率曲线图');xlabel('Frequency (Hz)');ylabel('Magnitude');

CODE 3

将signal tap采集的正弦波数据生成.txt文件保存，导入到matlab中分析频谱。clear all; close all；varnum; %调用signal tap 采集的正弦波数据，列矢量signal=transpose(VarName); %转置fs=50E6; %采样频率N=1024; %采样点数t=[0:1/fs:(N-1)/fs]; %采样时刻figure(1);plot(t,signal);title('正弦波信号');xlabel('Time (s)');ylabel('Magnitude');Y = fft(signal,N); %做FFT变换Ayy = abs(Y); %取模Ayy=Ayy/(N/2); %换算成实际的幅度Ayy(1)=Ayy(1)/2;F=([1:N]-1)*fs/N; %换算成实际的频率值，Fn=(n-1)*Fs/Nfigure(2);stem(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果axis([0 5E6 0 2500]);title('正弦信号频谱图');xlabel('Frequency (Hz)');ylabel('Magnitude');

CODE 4

clear;signal_tap_data_500KHZ;fs=50E6; %采样频率N=1024; %采样点数t=[0:1/fs:(N-1)/fs]; %采样时刻xlabel('Time (s)');ylabel('Magnitude(dB)');signal=transpose(signal_out); %矩阵转置 y=fft(signal,N);%对信号进行傅里叶变换yy=abs(y); %求得傅里叶变换后的振幅yy=yy/(N/2); %换算成实际的幅度figure(1);plot(t,signal);F=([1:N]-1)*fs/N; %换算成实际的频率值，Fn=(n-1)*Fs/Nfigure(2);stem(F(1:N/2),yy(1:N/2)); %显示换算后的FFT模值结果axis([0 5E6 0 2500]);title('幅度-频率曲线图');xlabel('Frequency (Hz)');ylabel('Magnitude');

CODE 5

clear;%%选频3MHz,不补零，原数据f_3M%调用f_3M文件signal=transpose(ROMOUT); %矩阵转置：将列矩阵转置成行矩阵fs=50E6; %采样频率N=1024; %采样点数t=[0:1/fs:(N-1)/fs]; %采样时刻figure(1);subplot(2,2,1);plot(t,signal);title({'3M-未补零-原数据长度 ';'时域波形图'});xlabel('Time (s)');ylabel('Magnitude');Y = fft(signal,N); %做FFT变换Ayy = abs(Y); %取模Ayy=Ayy/(N/2); %换算成实际的幅度Ayy(1)=Ayy(1)/2;F=([1:N]-1)*fs/N; %换算成实际的频率值，Fn=(n-1)*Fs/Nfigure(1);subplot(2,2,3);stem(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果axis([0 5E6 0 2500]);title('幅度-频率曲线图');xlabel('Frequency (Hz)');ylabel('Magnitude');

CODE 6

san; data=san(:,5);%取数据中的第五列 signal=transpose(data); %矩阵转置 fs=50E6; %采样频率 N=1024; %采样点数 t=[0:1/fs:(N-1)/fs]; %采样时刻 figure(1);plot(signal); xlabel('Time (s)'); ylabel('Magnitude'); Y = fft(signal,N); %做FFT变换 Ayy = abs(Y); %取模 Ayy=Ayy/(N/2); %换算成实际的幅度 Ayy(1)=Ayy(1)/2; F=([1:N]-1)*fs/N; %换算成实际的频率值，Fn=(n-1)*Fs/N figure(2); stem(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果 axis([0 5E6 0 2500]); title('幅度-频率曲线图'); xlabel('Frequency (Hz)'); ylabel('Magnitude');

CODE 7

clear;fre_500k_1024; %调用fre_500k_1024文件signal=transpose(RAOUT); %矩阵转置：将列矩阵转置成行矩阵fs=50E6; %采样频率N=1024; %采样点数t=[0:1/fs:(N-1)/fs]; %采样时刻figure(1);plot(t,signal);xlabel('Time (s)');ylabel('Magnitude');Y = fft(signal,N); %做FFT变换Ayy = abs(Y); %取模Ayy=Ayy/(N/2); %换算成实际的幅度Ayy(1)=Ayy(1)/2;%第一个点是直流分量，该点模值是直流分量的N倍，即实际幅度为模值的1/N倍。上式已经将所有点除以（N/2），在此基础上将该点除以2即可。F=([1:N]-1)*fs/N; %换算成实际的频率值，Fn=(n-1)*Fs/Nfigure(2);stem(F(1:N/2),Ayy(1:N/2)); %显示换算后的FFT模值结果axis([0 5E6 0 2500]);title('幅度-频率曲线图');xlabel('Frequency (Hz)');ylabel('Magnitude');

CODE 8

signal = transpose(VarName6) / 2048; %数据归一化处理CLK = 50 * 1000000;N = 512;t=[0:1/CLK:(N-1)/CLK]; %持续时间figure(1);plot(t,signal);xlabel('时间(s)');ylabel('value');title('时域图');%绘制时域图形Y = fft(signal,N); %FFT变换size = abs(Y); %取模size=size/(N/2); %换算成实际的幅度figure(2);stem(F(1:N/2),size(1:N/2)); %fft频域图axis([0 5E6 0 1.5]);title('频谱图');xlabel('频率(Hz)');ylabel('幅度');

上一篇：C# 双击全窗口/全屏显示和恢复

下一篇：PAT 1092