matlab用隐函数拟合的程序怎么编写
发布网友
发布时间：2022-09-25 02:33
共2个回答
热心网友
时间：2023-09-15 03:51
clear;clc
F=@(p,x)p(1)*exp((x(:,1)-x(:,2)*p(4))/0.026/p(6))+p(2)*exp((x(:,1)-x(:,2)*p(4))/0.026/p(7))+(x(:,1)-x(:,2)*p(4))/p(5)-p(1)-p(2)+p(3)-x(:,2);
x=[0,4.84827E-9
0.05,7.20861E-9
0.1,9.486E-9
0.15,1.33346E-8
0.2,1.85467E-8
0.25,2.54804E-8
0.3,3.57885E-8
0.35,4.94846E-8
0.4,6.86021E-8
0.45,9.34068E-8
0.5,1.36416E-7
0.55,2.18504E-7
0.6,4.51033E-7
0.65,1.2785E-6
0.7,4.57114E-6
0.75,1.70904E-5
0.8,5.62476E-5
0.85,1.45646E-4
0.9,2.93975E-4
0.95,4.91251E-4
1,7.21464E-4 ];
p0=[3 5 4 10 9 10 4];%%因为你这个函数很*，初值选择不好得不到好结果，这个初值不错
warning off
p=nlinfit(x,zeros(size(x,1),1),F,p0);
disp('I1、I2、IL、Rs、Rsh、n1、n2分别为：');
disp(num2str(p));
plot(x(:,1),x(:,2),'ro');hold on;
ezplot(@(x,y)F(p,[x,y]),[0,1,-1e-3,1e-3]);
title('曲线拟合');legend('样本点','拟合曲线')
热心网友
时间：2023-09-15 03:52
等式右侧可以看成双自变量x(1)=x,x(2)=y,这样就成了多变量拟合了。
【补充】初始值不好确定，但方法就是这样的：
clc; clear
A=[...
0 4.84827E-9
0.05 7.20861E-9
0.1 9.486E-9
0.15 1.33346E-8
0.2 1.85467E-8
0.25 2.54804E-8
0.3 3.57885E-8
0.35 4.94846E-8
0.4 6.86021E-8
0.45 9.34068E-8
0.5 1.36416E-7
0.55 2.18504E-7
0.6 4.51033E-7
0.65 1.2785E-6
0.7 4.57114E-6
0.75 1.70904E-5
0.8 5.62476E-5
0.85 1.45646E-4
0.9 2.93975E-4
0.95 4.91251E-4];
format short g
x=A(:,1);y=A(:,2);
X=[x,y];
%y=I1*exp((x-y*Rs)/0.026/n1)+I2*exp((x-y*Rs)/0.026/n2)+(x-y*Rs)/Rsh-I1-I2+IL
%I1=a(1),Rs=a(2),n1=a(3),I2=a(4),n2=a(5),Rsh=a(6),IL=a(7)
fun=inline('a(1)*exp((x(:,1)-x(:,2)*a(2))/0.026/a(3))+a(4)*exp((x(:,1)-x(:,2)*a(2))/0.026/a(5))+(x(:,1)-x(:,2)*a(2))/a(6)-a(1)-a(4)+a(7)','a','x')
a=nlinfit(X,y,fun,[ 1 1 1 1 1 1 1])
%{
figure(1)
plot(x,y)
figure(2)
semilogy(x,y)
figure(3)
loglog(x,y)
%}
