问题见图，希望能给出完整的程序和改善后的matlab程序,即fft变换后的程序和ctz细化后程序，

希望能给出完整的程序和改善后的matlab程序,即fft变换后的程序和ctz细化后程序，希望能尽快解决问题，朋友你好，祝你好运，希望可以加上注释，详细一点最好

clc;
fs=500;%采样频率
N=128;%采样点数
nfft=128;
n=0:1:N-1;%时间序列号
%n/fs:采样频率下对应的时间序列值
n1=fs*(0:nfft/2-1)/nfft;%FFT对应的频率序列
x=sin(2pi78n/fs)+sin(2pi82n/fs)+sin(2pi100n/fs);
XK=fft(x,nfft);%单边幅值谱
subplot(211);plot(n1,abs(XK(1:(nfft/2))));
title('直接利用FFT变换后的频谱');
f1=70;%细化频率段起点
f2=104;%细化频率段终点
M=85;%细化频段的频点数，(这里其实就是细化精度)
w=exp(-1i2pi(f2-f1)/(fsM));%细化频段的跨度(步长)
a=exp(1i2pif1/fs);%细化频段的起始点，这里需要运算一下才能代入czt函数
xk=czt(x,M,w,a);
h=0:1:M-1;%细化频点序列
f0=(f2-f1)/M*h+100;%细化的频率值
subplot(212);plot(f0,abs(xk));
xlabel('f');
ylabel('value');
title('利用CZT变换后的细化频谱');

一点想法，有错的话还望指正

你好，我是有问必答小助手，非常抱歉，本次您提出的有问必答问题，技术专家团超时未为您做出解答


本次提问扣除的有问必答次数，将会以问答VIP体验卡（1次有问必答机会、商城购买实体图书享受95折优惠）的形式为您补发到账户。


因为有问必答VIP体验卡有效期仅有1天，您在需要使用的时候【私信】联系我，我会为您补发。