btcq.net
当前位置:首页 >> mAtlAB 参数方程作图:x=(B+1)*(Exp(A*y)%1)/Exp(A)%1)%B*y,A... >>

mAtlAB 参数方程作图:x=(B+1)*(Exp(A*y)%1)/Exp(A)%1)%B*y,A...

对于这类作图问题,只能通过循环语句来完成.for a=-5:1:5 for b=-0.1:0.1:0.1 endend程序运行结果如有问题,可以进一步私聊.

解答如下:分两个步骤:第一步建立函数模型的M文件打开M编辑器,输入如下代码;function yd=model(beta0,x)%% L/1+a*exp(-b*x)a=beta0(1);b=beta0(2);L=beta0(3);x1=x;yd=L./(1+a.*exp(-b.*x1));保存这个M文件到work文件夹下,名称为:

>> [a,b]=solve('a*x^3-b*y^2=0','exp(a*x^2)-y-1=0','a,b')a =log(y+1)/x^2b =log(y+1)*x/y^2

两边取对数.lny=lna+bx y'=lny,a'=lna,b=b,y'=a'+b*x 求解最小二乘拟合,正则方程组可以得到a',b即:a,b 就可以了

clc;clearx=;%自己给数据y=;%自己给数据p=polyfit(x,log(y),1);b=p(1);a=p(2);yfit=a*exp(b*x);plot(x,y,'r-.')plot(x,yfit,'b-.')legend('拟合前','拟合后')

从图上可以看出,选用y=a*exp(b*x)+c进行拟合是不对的,应改用polynomial中的cubic polynomial,我试着拟合得:Linear model Poly3: f(x) = p1*x^3 + p2*x^2 + p3*x + p4 where x is normalized by mean 2007 and std 3.028Coefficients (with 95%

x会取到0,1/x无意义.主要是有一个无穷大的点,你可以在plot()函数后加一句:axis([-10,10,-10,10]);得到的图形就正常了以下的问题可以参考下:试一下max(y),为Inf,是x = 0时:x(find(y==max(y)));但是在matlab中exp(-1./0)的结果是0,exp(1./(-0.01))结果是Inf,而是矢量时exp(-1./x)得到的是x=-0.01那一项很小,x=0那一项为Inf

syms a b c d t x y ex1 = a*t+b-y; ex2 = c*t+d-x; t = solve(ex2,'t'); y = solve(subs(ex1),'y') y = (-a*d+a*x+b*c)/c

用符号运算f1='…'依次四个式子,然后再用solve(f1,f2…)试试,若不行就要我用@

可以,只要你的数据没有问题.

相关文档
ncry.net | zxpr.net | ldyk.net | msww.net | 596dsw.cn | 网站首页 | 网站地图
All rights reserved Powered by www.btcq.net
copyright ©right 2010-2021。
内容来自网络,如有侵犯请联系客服。zhit325@qq.com