文章编号:11360时间:2024-09-30人气:
仔细看看以下程序,慢慢体会。>> format rat>> x=1./magic(3)x = 1/811/6 1/31/51/7 1/41/91/2 >> x=1./magic(3)*pix = 355///678 355/339 71//791 355//1017 355/226 >> y=sin(x)y = 579/1513 *1/2 1170///7449 985//5564 1
假设数据集为x, y,则 2^y = 1 + x,转化为多项式拟合。程序示例:
y = log2(1+x);
plot(x, y, o); hold on;
t = polyfit(x, 2.^y, 1);
y_ = log2(t(1)*x+t(2));
plot(x, y_)
相关论文:Vc++下如何利用Matlab工具箱进行数字信号处理摘要: 本文详述了在Vc环境下如何利用Matlab工具箱进行数字信号处理,全文以Matlab工具箱中功率谱密度分析函数为例,介绍了通过Matlab自带的引擎、Matlab自身的编译器以及利用MathTools公司的Matcom进行对工具箱函数的调用。 关键词:Matlab M-文件 引擎 编译器 Matcom Vc++Matlab的信号处理工具箱是信号算法文件的集合,它处理的基本对象是信号与系统,信号处理工具箱位于目录、Toolbox\Signal下,利用工具箱中的文件可以实现信号的变换、滤波、谱估计、滤波器设计等。 在其它的环境如Vc下如果能调用Matlab工具箱中的文件,会大大地加快一些算法的实现,同时其可靠性也很高。 利用Matlab引擎 Matlab引擎采用客户和服务器计算方式,在运用中,Vc的C语言或C++语言的程序作为前端客户机,它向Matlab引擎传递命令和数据信息,并从Matlab引擎接收数据信息,它提供了下列几个函数: engOpen, engGetArray, engPutArray, engEvaString, engOutputBuffer ,engClose与客户机进行交互。 下面例程是在Vc下建一个基于对话框的应用程序,在对话框中设置一个Button控件OnMatlabEngine.,在对话框 文件中加入”engine.h” 和“math.h” 头文件,下面给出部分程序清单。 Void CtestmatlabDlg::OnMatlabEngine(){Engine *ep;mxArray* T=NULL,*result=NULL,*mFs=NULL,*mnfft= NULL;double datax[1024];char buffer[1024];for(int j=0;j<1024;j++)//注:如通过采集卡采集数据可将采集的数据放在datax[]数组中,此循环就不需要{double samt=(double)(1.0/1024); datax[j]=sin(2.0*63.0*samt*3.+1.15*3.);}double *pPxx,*pFxx;if(!(ep=engOpen( \0))){//打开Matlab引擎,建立与本地Matlab的连接fprintf(stderr,\n Cant start MATLAB engine\n);exit(-1);} double Fs[1]={1024};//因为Matlab所有参与运算的参数都是矩阵的形式,因而下列几行将参数转变double nfft[1]={1024};//成Matlab可接受的矩阵形式。 T=mxCreateDoubleMatrix(1,1024,mxREAL);mnfft=mxCreateDoubleMatrix(1,1,mxREAL);mFs=mxCreateDoubleMatrix(1,1,mxREAL);mxSetName(T,T);mxSetName(mnfft,mnfft);mxSetName(mFs,mFs);memcpy((char*)mxGetPr(T),(char*)datax, 1024*sizeof(double));memcpy((char*)mxGetPr(mnfft),(char*)nfft, sizeof(double));memcpy((char*)mxGetPr(mFs),(char*)Fs,1*sizeof(double));engPutArray(ep,T); //将转化的参数放入引擎中,此时可在Matlab command窗口下查看此参数engPutArray(ep,mnfft); engPutArray(ep,mFs);engEvalString(ep,[pxx,fo]=psd(T,mnfft,mFs);); //利用引擎执行工具箱中文件engOutputBuffer(ep,buffer,512); //如只想看显示图形,可将返回参数去掉,psd无返回参数缺省情况下会自动画图形result=engGetArray(ep,pxx);//取出引擎中的数据放在所指的区域中供后续处理pPxx=mxGetPr(result); result=engGetArray(ep,fo);pFxx=mxGetPr(result);engEvalString(ep,plot(fo,10*log10(pxx)););//利用引擎画图engEvalString(ep,title(功率谱分析););engEvalString(ep,xlabel(Hz););engEvalString(ep,ylable(db););mxDestroyArray(T); //释放内存mxDestroyArray(mFs);mxDestroyArray(mnfft);mxDestroyArray(result);engEvalString(ep,close;);engClose(ep);}上述程序在Vc下编译需要将 和两个动态库利用以下的命令:lib/def:<自己的Matlab的安装路径,下同>e:\ Matlab\extern\include\* /Machine:ix86 /out:*来生成程序所需的静态连接库和,将和所在的目录加入Vc++ project/link/object/library modules下即可。 利用Matlab自身的编译器调用工具箱中的函数 Matlab的编译器可将Matlab的M文件转换为为C或C++的源代码以产生完全脱离Matlab运行环境的独立的运用程序,但Matlab本身的资料说明编译器如用来建立独立的运用程序,不能调用Matlab工具箱中的函数,这非常不利于搞一些特殊的算法。 本人研究了一段时间发现,工具箱中的函数既然是M文件就一定可以用编译器来编译,以提供如Vc的调用函数,但是编译器只能编译一个独立的M文件,即这个 M文件不依赖于其他的M文件。 如果M文件中又调用了其他的M文件,可将被调用的M文件拷贝到调用M文件的相应位置,作适当的改动就可以用于编译器编译。 编译器不支持图形函数,所以M文件中如有图形函数需注释掉。 当Matlab的编译器mcc加入适当的参数-e(mcc –e *.*)或-p(mcc –p *.*)就可生成将输入的M文件转换为适用于特定运用的C或C++源代码。 这样如果要在Vc下编译通过,还需连入以下几个库, , , . . 以及Matlab C MATH库,建议采用前述的方法将动态连接改为静态连接。 对于C/C++编译环境的设置,在Matlab command窗口下运行mex –setup 然后依提示操作,而对于C/C++连接环境的设置,运行mbuild –setup依提示操作即可。 下面给出利用编译器将Matlab工具箱中psd.m文件生成可供Vc调用的函数。 将psd.m文件拷贝一份至Matlab\bin目录下,改写相应调用的M文件如nargchk.m, hanning.m等。 为生成的代码简洁,对于采集数据处理输入参数很明了的情况下可作大量的删减,最终使psd.m成为一个不依赖于其他M文件的独立的M文件,注意千万注释掉作图代码,最终改成如下形式,限于篇幅给出关键的几步: function [Pxx,f]=psd(Fs,nfft,noverlap,x)window=o,5*(1-cos(2*pi*(1:nfft)’/(nffft+1)));//hanning 窗dflag=’none’;window=window(;)………………………………….以上只要稍懂Matlab语言和信号处理知识就可完成这项工作。 假设上述代码重新存为testwin.m,在Matlab command 窗口下设置好环境参数运行mcc –e testwin,则可在Matlab\bin下生成testwin.c ,如运行mcc –p testwin 则生成. Vc下建立一个基于对话框的文件,然后在对话框里加一个Button控件OnButtonPsd 将上述生成的.c文件的头文件加入到工程的中,且将#ifdef_cplusplus extern “c”{#end ifc代码声明加入Vc的包含文件和生成的.C的包含文件之间将#ifdef_cplusplus}#end if加入文件未尾为了简洁且便于处理将生成的c函数稍改动,给出部分代码如下: void CTestpsdwinDlg::OnButtonPsd(){mxArray* x_rhs_;//指向采集数据存放区Fs=;//数据采集的频率 nfft=1024;//1024点的fftdouble datax[1024]//采集的数据x_rhs_mxCreateDoubleMatrix(1,1024,mxReal);memcpy(mxGetPr(x_rhs_),datax,1024*sizeof(double));noverlap=512;……………….………………(&Pxx,&Spec);mccCopy(&f,&frevgg_vector);for(int j=0;j<(int)(nfft/2+1);j++){datap[j]=mccGetRealVectorElement(&Pxx, (j+1));//功率谱密度存于datap[]数组dataf[j]=mccGetRealVectorElement(&f, (j+1));//相应频率存于数组dataf[]中}mccFreeMatrix(&Pxx);………………(WM_PAINT,0,0);//利用Vc下的图形函数画图Return;}如上生成的程序可读性不太好,而生成的c++代码则可读性较好,但千万注意只能用 Matlab的MATH库,不可用c++的MATH库,否则编译会出错,限于篇幅在此不述。 3)利用Matcom调用工具箱中的函数Matcom编译M文件,先将M文件按照与Matcom的cpp库的对应关系翻译为cpp源代码,然后用对应版本的c编译器将cpp文件编译成相应的exe或dll文件,所以第一次运行要指定c编译器的路径,否则无法编译,指定好的编译信息就写在Matcom\bin\文件中,不过这一步按装matcom时,它自动寻找编译器并将其写入文件中,matcom4.5版中使用TeeChart3.0 OCX控件,因而它支持图形操作。 我们依然用上述的testwin.m文件,不要将图形函数注释掉,利用Mideva来生成可被Vc调用的信号处理程序。 运行Mideva在主界面上直打开M文件,在菜单中选择compile to dll,输入testwin..在Matcom debug目录下可以找到这样的几个文件,testwin.c ,testwin.h,,,,等。 将上述和testwin.h加入工程中,project/add to project/files并且在相应的文件中加入”stdafx.h” 加连接库:Tools\option\directory\ , 选include选项,加入e:\matcom45\lib (包含matcom.h) library选项,加入e:\matcom45\lib4) project\add to project\files 文件类型选项选()将e:\matcom45\lib\加入工程中编译运行。 相应代码如下:void CtestmatcomDlg::OnpsdButton(){double datap[512],dataf[512];initM(MATCOM_VERSION);//初始化matcom库Mm Fs,nfft,noverlap;//创建矩阵Mm x=zeros(1,1024);Fs=1024;nfft=1024;noverlap=128;dMm(Pxx_o);dMm(f_o);//创建并命名矩阵datax[];//数据采集的数据存于此数组中for(int i=1;i<=1024;i++){x.r(1,i)=datax[i+1];//给x阵赋值}testwin(Fs,nfft,noverlap,x,i_o,Pxx,f_o);//matcom生成的函数for(i=0;i<513;i++){//取出功率谱密度分析结果dataf[i]=f_o.r(i+1,1);datap[i]=Pxx_o.r(i+1,1);}exitM();return;}可见利用Matcom进行M文件转换非常的容易,生成的代码可读性很好,以上的转换同时生成了可供Vc调用的动态连接库,其使用和一般的动态库一样使用。 同时需指明Matcom不仅可转换独立的不依赖于其它M文件的M文件,同时可转换调用其它M文件的M文件嵌套。 条件是这此M文件在同一个目录下面,如前所述的psd.m可直接用上述方法转换,生成了多个重载形式的psd函数结论: 利用Mtlab引擎调用工具箱中的函数可节省大量的系统资源,应用程序整体性能较好,但不可脱离Matlab 的环境运行。 用Matlab编译器进行工具箱函数的调用,须转换相应的M文件使其成为独立的M文件,且不支持图形函数,转换的代码可读性不太好。 用Matcom 进行转换非常方便,生成的代码可读性很好,支持图形函数,且代码执行的速度比不转换平均要快1.5倍以上。 以上程序在Vc++ 6.0,Matlab5.2,Matcom4.5中调试通过,以上方法在工程实践中已得到很好的运用。 仅供参考,请自借鉴希望对您有帮助
两种方法,一种是把X轴设置为plot(log10(x),y);
第二种semilogx(x,y)
利用plot函数画出图像,图像为空白,解决方法如下:
直接按照你的程序运行,y的值除了第一个都是-Inf,所以无法作图。 如果你需要计算精确值得话,可以考虑计算y的对数。 如果只是想看看x和y的关系曲线,那可以把t的值减小一些,比如改成:
这样就可以画图了,只是这样x和y的值都会改变。 而x和y都是t的非线性的表达式,这样减小t的值,y对x的曲线形状肯定也改变了,简单来说,就是只画了曲线的前分之一。 这样画成的曲线是这样的:
可以大致看出这个曲线的走势,就是往右边越来越陡峭的下降,直到-inf。 看到这里肯定有疑问,为什么y的值其实不大(-0.3),却会成为-inf?仔细看下题主给出的这个式子,为了便于描述,我把这些分式表示的数字都计算成小数了(造成的误差不大,不会影响结果):y=0.1997*t-0.0040*(log(6.2687*(387.5939*exp(49.4697*t)-387.1945).^2);再简化一下:y=A*t-B*(log(C*(D*exp(E*t)-F).^2);可以看出我们的计算经过了:先求exp->减去常数->平方->取对数 这个过程,其实计算前后数量级应该变化不大,但是在计算过程中求exp再平方这个步骤的中间值很大,很有可能导致溢出,所以就会产生-inf了。 找到了原因,接下来我们想办法解决它。 首先根据对数运算的性质,平方可以取出来:y=A*t-2*B*(log(C*(D*exp(E*t)-F));然后,系数C乘在大数的前面也可能导致溢位,所以也取出来:y=A*t-2*B*(log(D*exp(E*t)-F)+log(C));系数D的值也不小,一样提出来: y=A*t-2*B*(log(exp(E*t)-F/D)+log(D)+log(C));在这里D的值为387.5939,F的值为387.1945,F/D约等于1(小于1),由于t>0,这里exp(E*t)总是大于1的。 当exp(E*t)的值不是很大时,F/D不能忽略,而当exp(E*t)很大时,F/D就可以忽略了。 现在让我们看看这个分界点取在多少比较好。 先画出 :y=A*t-2*B*(log(exp(E*t)-F/D)+log(D)+log(C));
可以看出虽然t的取值是从0到20,但是画到t=14.34的时候,就没有了,这说明后面的数值是-inf了。 这个时候exp(E*t)等于多少呢?exp(E*14.34)=1.2205e+的308次方!当然这是Matlab在我机器上所能计算的最大的double型数值,这个数字比1大太多太多了,完全可以忽略掉这个1了。 那我们就取这个分界点在14.34吧,当t大于14.34时,忽略F/D,这样y的计算式可以化简为:y=A*t-2*B*(E*t+log(D)+log(C));这个形式就不会有巨大的中间数值,可以放心计算了。 通过这些分析,可以把y的计算公式写成如下的分段形式:if (t<14.34)y=A*t-2*B*(log(exp(E*t)-F/D)+log(D)+log(C));elsey=A*t-2*B*(E*t+log(D)+log(C));end这样就可以画出x-y曲线:
这就是题主想要画的曲线了。 由于t的取值步长是1000,所以从第二个点开始,就是按照后半段公式计算的y,y就是t的线性函数。 而x是t的对数函数,于是就画出了上面这样的图形。 其实这个图形在t的值比较小的时候,是有其它的形态变化的,参见第一个图。
内容声明:
1、本站收录的内容来源于大数据收集,版权归原网站所有!
2、本站收录的内容若侵害到您的利益,请联系我们进行删除处理!
3、本站不接受违法信息,如您发现违法内容,请联系我们进行举报处理!
4、本文地址:http://www.jujiwang.com/article/5b67a880fa40ec1703ad.html,复制请保留版权链接!
牛腩新闻发布系统是一款功能强大的工具,可帮助企业与潜在客户建立联系并推动业务增长,以下是一些方法,您可以利用牛腩新闻给您的电子邮件列表或潜在客户清单,6.跟踪您的结果并进行调整牛腩新闻发布系统提供详细的分析,使您可以跟踪您的结果并根据需要进行调整,这将帮助您优化您的活动并最大化您的投资回报,结论牛腩新闻发布系统是一个强大的工具,可以帮...。
本站公告 2024-09-25 17:07:48
前言作为两种流行的Web技术,IIS,Internet信息服务,和PHP,超文本预处理器,在不同的领域发挥着至关重要的作用,IIS是一个强大的Web服务器,以其稳定性、可靠性和对Windows操作系统的原生支持而闻名,另一方面,PHP是一种通用脚本语言,广泛用于开发动态且交互式Web应用程序,将IIS与PHP集成在一起,可以充分利用这...。
本站公告 2024-09-25 13:18:44
在线投资组合或网站是向潜在投资者展示您的投资策略和业绩的一种绝佳方式,它可以帮助您吸引新客户、建立您的声誉并增加您的业务,创建在线投资组合或网站的好处它可以帮助您吸引新客户,它可以帮助您建立您的声誉,它可以增加您的业务,它可以使您更容易向潜在投资者展示您的投资策略和业绩,它使您能够接触到全球受众,创建在线投资组合或网站的步骤选择一个平...。
本站公告 2024-09-23 08:41:27
引言React是一个用于构建用户界面的流行JavaScript库,近年来,React团队不断发布新特性和改进,以增强开发人员体验并提升应用程序性能,在文章中,我们将探讨三个重要的React技术前沿,ReactHooks、Suspense和Fiber,ReactHooksReactHooks是在Reactv16.8中引入的一项重大特性,...。
本站公告 2024-09-16 13:19:57
BSDC库MicrosoftVisualC,forLinuxEclipseCDTJetBrainsCLion使用C语言库和工具在C语言项目中使用库和工具时,请记住以下几点,选择满足您的项目需求的库和工具,按照库或工具的文档进行适当的安装和配置,使用库或工具时,遵守其许可条款,不断研究和了解C语言库和工具的最新发展,结论海量的C语言库...。
最新资讯 2024-09-13 11:22:08
简介Rank函数是一个强大的数据转换函数,在数据建模中有着广泛的应用,它可以根据指定的列或表达式对行进行排序,并返回每个行的排名,Rank函数不仅可以用于构建预测模型,还可以提升模型的预测精度,构建预测模型1.数据准备使用Rank函数构建预测模型的第一步是准备数据,这包括清洗、转换和标准化数据,2.特征工程下一步是进行特征工程,这涉及...。
最新资讯 2024-09-12 12:59:45
前言监听循环是服务器端不断监听服务器套接字上的连接请求的过程,它是服务器端与客户端通信的基础,也是确保服务器端能够及时响应客户端请求的关键,监听循环的工作原理监听循环通常在以下步骤中进行,服务器端创建一个服务器套接字,并将其绑定到一个特定的IP地址和端口号,服务器端调用listen,函数,指定服务器套接字可以挂起的最大连接请求数,服...。
本站公告 2024-09-11 07:57:03
淘宝导航代码是淘宝店铺流量的重要来源之一,它可以通过多种方式将用户引流到您的店铺,但是,许多卖家对淘宝导航代码的使用存在疑惑,这可能会影响店铺的流量,本文将对淘宝导航代码进行详细讲解,扫除店铺流量的拦路虎,1.什么是淘宝导航代码,淘宝导航代码是一段可以放置在任意网页中的代码,当用户点击代码中的链接时,将跳转到您的淘宝店铺,淘宝导航代码...。
技术教程 2024-09-10 22:03:31
引言作为一名VB开发人员,能够连接到SQL数据库至关重要,本手册旨在提供有关如何使用VB.NET连接到SQLServer数据库的分步指南,先决条件VisualStudio2019或更高版本SQLServer数据库步骤1,添加ADO.NET引用1.在VisualStudio中打开VB.NET项目,2.右键单击,引用,文件夹,然后选择,添...。
本站公告 2024-09-10 08:24:45
在本文中,我们将引导您逐步创建强大的面部检测和识别算法,我们将使用Java编程语言,它是一种适用于广泛应用的通用语言,通过按照这些步骤操作,您可以创建自己的面部检测和识别系统,用于各种任务,如安全、医疗和娱乐,第1步,安装OpenCvOpenCV是一个开源计算机视觉库,它提供了用于面部检测和识别的功能,要安装OpenCV,请按照下面给...。
最新资讯 2024-09-07 20:25:31
欢迎来到学校网站构建指南,本指南将介绍学校网站构建过程中你需要了解的一切,从选择正确的平台到设计和开发你的网站,选择正确的平台在开始构建网站之前,你需要选择一个适合你的学校的平台,有许多不同的平台可供选择,每个平台都有其优点和缺点,以下是一些最受欢迎的学校网站平台,WordPress,WordPress是一个免费的开源平台,被广泛用于...。
最新资讯 2024-09-05 20:15:11
加强抵御措施的意义2013年10月28日发生的北京恐怖袭击事件震惊了中国和世界,此次袭击造成5人死亡,40多人受伤,凸显了恐怖主义的严重威胁,为了应对这一威胁,中国政府采取了一系列措施来加强抵御措施,防止未来发生类似事件,加强抵御措施具有以下重要意义,保护人民生命安全,抵御措施旨在保护人民免受恐怖袭击的伤害,确保他们的生命和安全,维护...。
互联网资讯 2024-09-04 03:19:26