频域积分的VC++.docVIP

频域积分的VC++程序 解决以下几个问题 频域积分的原理 matlab程序 完成 将频域积分的matlab涉及到的函数，翻译成c++ 加入到类中 完成 加入的函数、变量的调用和传递 这是重点 Updata(true)； Updata(false); 2、matlab 和 c++、本系统变量的程序对应 x=x_accr; % 原始信号 % 系统的通道数据 m_Data[0]和m_Data[1]; n=length(x); % 原始信号的数据长度 % 接受从类Cshowview传递的数据长度和采样频率变量，及变量在两个类中的传递。 t=0:1/sf:(n-1)/sf; % 模拟信号的时间间隔 % 系统中是从现场采集的数据。 nfft=2^nextpow2(n); % 适用于FFT分析的数据长度 % 系统中有这个函数，已经添加好，注意调用方式。 y=fft(x,nfft); % FFT变换 % 系统中有FFT的c++ 函数，但如何处理这个返回值？ df=sf/nfft; % 频率分辨率 %系统中已经添加 ni=round(fmin/df+1); % 计算指定频带对应频率数组的下标 na=round(fmax/df+1); % 四舍五入 取最近的整数 %系统的截止频率设置对话框的操作-设定截止频率、积分次数 dw=2*pi*df; % 计算原频率间隔 w1=0:dw:2*pi*0.5*sf; % 建立正的离散原频率向量 w2=-2*pi*(0.5*sf-df):dw:-dw; % 建立负的离散原频率向量 w=[w1,w2]; %将正负组成一个数组 % 以上向量、数组如何翻译成C++？ it=2; %积分次数 w=w.^it; % 数组的平方，如何翻译成c++? c=1000*1000; a=zeros(1,nfft); % 建立一个长度为nfft的零向量 a(2:nfft-1)=y(2:nfft-1)./w(2:nfft-1); % 如何翻译成c++ if it==2 % 下面的是相位、幅值的变换，翻译成c++ y=-a; else a1=imag(a); a2=real(a); y=a1-a2*i; end a=zeros(1,nfft); a(ni:na)=y(ni:na); a(nfft-na+1:nfft-ni+1)=y(nfft-na+1:nfft-ni+1); y=ifft(a,nfft); % 调用逆傅里叶变换 y=real(y(1:n))*c; figure(2); subplot(2,1,1); plot(t,y); % 画出时域波形 a(2:nfft-1)=y(2:nfft-1)./w(2:nfft-1); % 如何翻译成c++ 上述这条语句困扰我很久了。 一、对于FFT的调用，参考rotview6ok系统，单通道频谱分析 1、注意私有 变量的定义： complexdouble * m_pFreData; 对于视图类的初始化，调用的消息为：OnInitialUpdate(); 对于对话框类的初始化调用的消息为：WM_INintilog 2、对于公有函数，关注： InitShow()和ChangeDataLength()、DrawFftGraph() 比如 void CSingleFrequencyDlgView::ChangeDataLength() // 改变数据长度之后，调用fft { delete [] m_pFreData; delete [] m_pFreData1; m_pFreData=NULL; m_pFreData1=NULL; m_pFreData=new complexdouble[m_FftLength]; FFTTrans(m_pData,m_FftLength,m_pFreData,m_Ampty); // 包含头文件，返回值是 DrawFftGraph(); // 画频谱图 } 3、使 编辑框（控件）显示数值： if (myDoc-m_data.GetDataLength()!=0) { m_AccurateFre=m_frequency-m_AccurateFre; CString tmp; tmp.Format(%.2f,