機械测试MATLAB的应用.docVIP

首先，采样点数不变，改变采样频率 信号频率取50Hz 采样频率取50、100、200、1000四种情况。 采样点数N=128. 程序： fs=50（100、200、1000）; %设定采样频率 N=128; n=0:N-1; t=n/fs; f0=50; %设定正弦信号频率 %生成正弦信号 x=sin(2*pi*f0*t); figure(1); subplot(111); plot(t,x);%作正弦信号的时域波形 xlabel(t); ylabel(y); title(正弦信号y=2*pi*50t时域波形); grid; f0=50Hz时，产生混叠 f0=100Hz时，理论上恰好不产生混叠，但实际上已经错误。 f0=200Hz时，比较尖，不平滑 f0=1000Hz时，比较平滑 下面做FFT变换 两种：一种是MATLAB自带的FFT变换程序，另一种是自编的FFT变换程序。 先介绍自编程序： y=myfft（x，n） function y= myfft(xr,n) p=0:n-1;% 开始倒位序 nu=log2(n); p1=p; b=zeros(1,n); for t=1:nu; p2=floor(p1/2); b=b*2+(p1-2*p2); p1=p2; end; yr(p+1)=xr(b+1); xr=yr; % 倒位序结束 t=0:n/2-1; %计算因子w开始 (只计算w0到w n/2-1) for v=0:n/2-1; w=exp(-2*i*pi*t/n); end; %计算因子w结束 for m=1:nu;% 计算x(k)开始 h=2^(m-1); k=1; while(kn+1) for t=1:h; y=bitshift(k-1,nu-m,nu)+1; %求w的幂次数 xch(k)=xr(k)+w(y)*xr(k+h); k=k+1; end; for t=1:h; y=bitshift(k-1-h,nu-m,nu)+1; %求w的幂次数 xch(k)=xr(k-h)-xr(k)*w(y); k=k+1; end; end; xr=xch; end;% 计算x(k)结束 y=xr%输出变换后的结果 给出画图程序： fs=(自取);%设定采样频率 N=(自取); n=0:N-1; t=n/fs; f0=50;%设定正弦信号频率 %生成正弦信号 x=sin(2*pi*f0*t); %进行FFT变换并做频谱图 y=fft(x,N)（或者是myfft）;%进行fft变换 mag=abs(y);%求幅值 f=(0:length(y)-1)*fs/length(y);%进行对应的频率转换 figure(1); subplot(111); plot(f,mag);%做频谱图 axis([0,100,0,80（可改）]); xlabel(频率(Hz)); ylabel(幅值); title(正弦信号y=sin(2*pi*f0*t)幅频谱图); grid; 给定信号频率50Hz，采样频率1000Hz，改变采样点数N N=128 N=512 N=1024 给定信号频率50Hz，采样频率200Hz，采样点数N=128 给定信号频率50Hz，采样频率120Hz，改变采样点数N=128。发生错误了。 给定信号频率50Hz，采样频率100Hz，改变采样点数N=128。（？）