m文件如下:
function y=modely(beta,X)
beta(1)=1;
beta(2)=1;
beta(3)=1;
beta(4)=1;
beta(5)=0.45;
beta(6)=0.35;
X=[0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
];
x1=X(:,1);
x2=X(:,2);
x3=X(:,3);
y=@(beta,X)beta(1)/(beta(2)*(pi/2)^(1/2))*exp((-abs(beta(3))*(x3-abs(sin(beta(5)))).^2-abs(beta(4))*(x1-abs(sin(beta(6)))).^2)./beta(2)^2);
主程序如下:
X=[0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.643685573 0.138181349 0.218133078
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.475664857 0.216601181 0.307733962
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.313687287 0.314252942 0.372059771
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.193221192 0.348425311 0.458353497
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
0.076187423 0.381624039 0.542188538
];
y=[0.049208998
0.067802007
0.079395017
0.097988027
0.125581036
0.142174046
0.160767055
0.191360065
0.271802133
0.282409663
0.280017192
0.296624722
0.314232251
0.322839781
0.33244731
0.35005484
0.336232701
0.342923873
0.350615045
0.372306217
0.386997389
0.397688561
0.415379733
0.436070905
0.352379986
0.359913653
0.37144732
0.386980987
0.404514653
0.42204832
0.437581987
0.467115653
0.617609347
0.626993082
0.630376818
0.647760553
0.661144288
0.671528024
0.686911759
0.700295494
0.21185162
0.213260648
0.220669676
0.223078704
0.226487732
0.234896761
0.248305789
0.260714817
];
beta0=[0.5,0.5,-2,-2,0.5,0.5];
[beta,r,j] = nlinfit(X,y,model,beta0);希望加上去得约束条件是:
0
0