请教MATLAB编程高手,如何用matlab自定义函数拟合曲线,求详细程序,万分感谢
第一张图中T即是曲线的拟合表达式,第二张图是n和α的一般表达式,且n和α只与λ有关,s,d都是常数,即通过拟合n和α来达到拟合T,数据只有(λ,T),还想请问MATLAB...
第一张图中T即是曲线的拟合表达式,第二张图是n和α的一般表达式,且n和α只与λ有关,s,d都是常数,即通过拟合n和α来达到拟合T,数据只有(λ,T),还想请问MATLAB高手,如何编写程序呢?万分感谢
展开
展开全部
实现代码:
clc,clear all,close all
a1=?;a2=?;b1=?;b2=?;d=?;s=?; %需要已知
lambda=[.....]; %数据
T=[.....]; %数据
n=a1./lambda.^2+a2;
alpha=10^(b1./lambda.^2+b2);
A=16*n^2*s;
B=(n-l)^3*(n+s^2);
C=2*(n^2-l*(n^2-s^2));
D=(n-l)^3*(n-s^2);
phi=4*pi*n*d./lambda;
x=exp(-alpha*d);
a(1)=n;a(2)=alpha;
func=@(a,lambda)(A.*x./(B-C*x*cos(phi)+D.*x.^2));
a0=rand(1,2)
a= nlinfit(lambda,T,func,a0);
或
a=lsqcurvefit(func,a0,lambda,T)
vpa(a,10)
上述代码,需要你提供数据来验证。
clc,clear all,close all
a1=?;a2=?;b1=?;b2=?;d=?;s=?; %需要已知
lambda=[.....]; %数据
T=[.....]; %数据
n=a1./lambda.^2+a2;
alpha=10^(b1./lambda.^2+b2);
A=16*n^2*s;
B=(n-l)^3*(n+s^2);
C=2*(n^2-l*(n^2-s^2));
D=(n-l)^3*(n-s^2);
phi=4*pi*n*d./lambda;
x=exp(-alpha*d);
a(1)=n;a(2)=alpha;
func=@(a,lambda)(A.*x./(B-C*x*cos(phi)+D.*x.^2));
a0=rand(1,2)
a= nlinfit(lambda,T,func,a0);
或
a=lsqcurvefit(func,a0,lambda,T)
vpa(a,10)
上述代码,需要你提供数据来验证。
追问
万分感谢你详细的回答,可能我表达的不够清晰,让你理解有误。这里的a1,a2,b1,b2是需要拟合来确定的,不是已知的,给出的数据仅是(lambda,T),即通过该数据拟合确定a1,a2,b1,b2。也就是说只能通过拟合来确定a1,a2,b1,b2,就可以求出T的具体表达式,ba wu qi wu san yi si san wu 这是我的koukou,期待详细交流
展开全部
X=[65.31 44.88 39.92 34.73 29.51 24.17 18.20 12.33 6.02 0.64];Y=[0.98 2.02 2.89 2.87 2.81 5.02 5.11 4.75 5.584 5.17];fun=inline('a(1)*X+a(2)./X+a(3)','a','X');b=[-1 0 1]a = nlinfit(X,Y,fun,b);vpa(a,10)A=a(1),B=a(2),C=a(3) A= -....
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询