有点相似DFT运算呢,matlab R2010如何读取NC文件,如果不斟酌运算速度的话仍是好实现的。 首先函数等号左边应当是X(e^jw)吧? 等号右边e的负jwn次方能够依据欧拉公式 e^jx=cosx+jsinx 变换成正弦和余弦函数,mapinfo的图如何能够在一张纸上打印出来?,比方说w=-20则X(e^jw)=∑x[n]*[cos(20n)+jsin(20n)] n ∈(-20,20)根据你供给的前提x[n]是一个实数序列,则程序须要两个变量分辨存储实部跟虚部,采取轮回方法从-20加到20,此为第一层循环,得到X(e^j(-20)),然后在将w从-20循环到20顺次得到各个w值对应的函数值,最后便可以用plot函数将这些值话出来,当然程序履行到这一步函数的成果将虚部和实局部开来寄存的,可以求把实部虚部门别平方乞降后开放得到复数的幅值然后画点,得到幅频谱,也可以求出不同w对应的相位而后画点得到相位谱。 第一,对你的问题中w的取值距离还有待商议。第二,这种盘算方式比拟费时,不适于实时剖析,当然耗时也是绝对而言,如果x[n]是实数序列的话程序执行完也要不了多久。 假如还不清楚的话可以参考一下数字信号处置的DFT(离散傅里叶变更)部分和FFT疾速算法部分。 这个谜底只是探讨问题的解决计划,mapgis里特殊线型如何设置,如下程序供参考: function Y=myfun(X,w) Y=[],MAPGIS属性结构如何修改?;n=1:length(X); for m=1:length(w) Y=[Y,X.*exp(-1*w(m)*sqrt(-1)*(n-1)/2)]; end 保留此函数为myfun,mapgis如何建库呀,用什么软件.m文件,当你给定数组X和w时,调用函数输出的就是所求,mapgis图层不能合并问题及拓扑错误如何修改?。 恳请验证!
(责任编辑:admin)
|