wigner.m

内燃机转子仿真

clear all;
clc;
format long;
close all;
load 'bneiranji1';
load 'cneiranji1';

x0=cxyz1(2,2000:2500);
t=(2000*h:h:2500*h);
nn=length(x0);
nfft=2^10;
flag=0;
[wx, waxis] = wig2c(x0,nfft,0.06,flag);
surf(waxis(1,nfft/2+1:nfft)*(1/h),t,wx(:,nfft/2+1:nfft));