对12bit的AD建模后matlab写程序对其进行FFT，计算SNR并输出频谱图，程序该怎么写？？？急

你是不是用simulink将数据输出到workbench？再用matlab计算其SFDR,SNR....等等参数？

最近也在做DAC的动态校准的东西，就把代码给你吧。我注释了一部分，你根据自己的需要自己改回来。

%load ScopeData

dynamic_data=ScopeData.signals.values

data_size=size(dynamic_data)

N=max(data_size)

sample_number=input('输入采样点数：')

if sample_number>N

display('样本数据数量不够,请检查')

else data=dynamic_data(1:sample_number)

N=sample_number

end

F_samples=input('输入采样频率:')

Mean_code=mean(data)

V=data-Mean_code

V=V.*hanning(N)

dout_spect=fft(V)

dout_dB=20*log10(abs(dout_spect))

max_dB=max(dout_dB(1:N/2))

x=(0:N/2-1).*F_samples/N

y=dout_dB(1:N/2)-max_dB

figure(1)

plot(x,y)

axis([0,F_samples/2,-120,0])

grid off

title('FFT图')

xlabel('模拟输入频率（Hz）')

ylabel('幅值（dB）')

F_in=find(dout_dB(1:N/2)==max_dB)

span=max(round(N/200),5)

spanh=2

spectp=(abs(dout_spect)).*(abs(dout_spect))

Pdc=sum(spectp(1:span))

%Ps=sum(spectp(F_in-span:F_in+span))

F_harmonic=[]

P_harmonic=[]

for har_number=1:10

tone=rem((har_number*(F_in-1)+1)/N,1)

if tone>0.5

tone=1-tone

end

F_harmonic=[F_harmonic tone]

har_peak=max(spectp(round(tone*N)-spanh:round(tone*N)+spanh))

har_bin=find(spectp(round(tone*N)-spanh:round(tone*N)+spanh)==har_peak)

har_bin=har_bin+round(tone*N)-spanh-1

P_harmonic=[P_harmonic sum(spectp(har_bin-1:har_bin+1))]

end

disp(F_harmonic)

disp(P_harmonic)

%

% Pd=sum(P_harmonic(2:5))

%

% Pn=sum(spectp(1:N/2))-Pdc-Ps-Pd

format

% SNR=10*log10(Ps/Pn)

%

% SINAD=10*log10(Ps/(Pn+Pd))

SFDR=10*log10(P_harmonic(1)/max(P_harmonic(2:5)))

% THD=10*log10(Pd/P_harmonic(1))

%

% ENOB=(SINAD-1.76)/6.02

text(x坐标,y坐标,['参数a的值为',num2str(a)]

注意上面的第三个参数是字符串，如果要将数值输出，要用num2str函数，这这个数值先转成字符串才可以输出。

最后用中括号括起来，表示字符串的拼接。

以下是一份MATLAB程序，使用筛法求得小于10^8的雷打数：

```matlab

n = 1e8% 设置筛法范围

isPrime = true(n, 1)% 初始化数组为全部为true

isPrime(1) = false% 定义1不是素数

% 首先从2开始，将所有的素数的倍数标记为非素数

for i = 2:sqrt(n)

if isPrime(i)

isPrime(i * i : i : n) = false

end

end

% 输出小于10^8的雷打数

for i = 1:length(isPrime)

if mod(i, 4) == 1 &&isPrime(i)

fprintf('%d\n', i)

end

end

```

在程序中，首先定义了一个筛法范围 `n = 1e8`，然后利用逻辑数组 `isPrime` 来记录每个数是否为素数。初始化时，所有的数默认为素数，将数组全部初始化为 `true`；同时，定义数字1不是素数。之后，从2开始循环，将所有的素数的倍数标记为非素数。具体实现的方法是，如果检测到一个素数i，则将i * i及其之后的所有i的倍数都标记为非素数。

最后，再次循环判断并输出所有小于10^8的雷打数。因为所有的雷打数的形式为4k+1，所以对于满足此条件的素数，都可以视为雷打数。程序使用 `mod(i, 4) == 1` 检查当前数是否满足4k+1的形式，并使用 `isPrime(i)` 检查该数是否为素数。如果两个条件都满足，则输出该数即可。

需要注意的是，在程序中使用了 `sqrt(n)` 函数，以降低时间复杂度。因为n以内的素数最大可能为 `sqrt(n)`，所以只需要进行到 `sqrt(n)` 的筛法，即可得到全部小于n的素数。

---