matlab 拟合实验数据

来源:百度知道 编辑:UC知道 时间:2024/05/21 09:14:46
实验数据:g1=1.65,g2=1.78,g3=2.0023,A1=252,A2=255

g1=(g3+2n)[cos(a)]^2-g3*[sin(a)]^2;
g2=g3*[sin(a)]^2-2n*sin(a)cos(a)
A1=f1(P,K,a),A2=f2(P,K,a)

a=f3(D,Ds,Dt),D=n*654,Ds=f4(y),Dt=f5(y),f12345形式已知,

n,y,K为可调参量,拟合实验数据,应该怎么做
R0=1.95*10^(-8),R1=R0*(1+0.7*y^2),R2=R0*(1+0.5*y^2),
D1=(4/7)*14930*((R0/R2)^3-(R0/R1)^3),D2=(2/7)*1493*((R0/R2)^5-(R0/R1)^5),
T0=654,P0=314*10^(-4),T1=n*T0,p=n*P0,a=0.5*atan(2^(0.5)*T1)/(0.5*T1-(3*D1-5*D2)),
g0=2.0023,g1=(g0+2*n)*cos(a)^2-g0*sin(a)^2,
g2=g0*sin(a)^2-2.828*n*sin(a)*cos(a),
A1=p*(k*cos(2*a)+2*((3/9.898)*sin(2*a)-(2/7)-(6/7)*cos(a)^2)),
A2=p*(k*sin(a)^2-2*((1/7)*sin(a)^2+(11/19.796)*sin(2*a))),
n,k,y为可调参量,abs(g1)=1.65,abs(g2)=1.78,abs(A1)=252*10^-4,abs(A2)=255*10^-4

http://www.blog.edu.cn/user2/joyfull/archives/2006/1417064.shtml

给你个提示,
输入:

x=[1 2 3 4 5 10 15 20 25 30 35 40 45 50 55 60 65 70 75 80 85 90 95 96 97 98 99 100 ]

y=[ 0 0 0 0.5 1.5 5.5 8 10 11 12.5 14.5 15.5 16.5 17 18 18.5 19 18 16 13 10 6 4 2.5 2.5 1 1 0]

plot(x,y,'ro')

这样就画出图来了。

在图形窗口中,选Tools > Basic Fitting 打开简单拟合对话框。

可以选择几个公式。
比如quadratic得下面结果

y = p1*x^2 + p2*x +
p3

Coefficients:
p1 = -0.0075168
p2 = 0.77989
p3 = -2.2647

Norm of residuals =

6.0295

或者cubic
y = p1*x^3 + p2*x^2 +
p3*x + p4

Coefficients:
p1 = -4.0802e-005
p2 = -0.0013296
p3 = 0.54468
p4 = -0.91163

Norm of residuals =
3.9283

加油吧,希望有所帮助

ma