m文件如下:
function y=modely(beta,X);
X=[0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.643685573 0.138181349 0.218133078
0.475664857 0.216601181 0.307733962
0.313687287 0.314252942 0.372059771
0.193221192 0.348425311 0.458353497
0.113377277 0.378605746 0.508016977
0.877547538 0 0.122452462
0.674765579 0.096568907 0.228665514
0.497203845 0.181127428 0.321668727
0.340432788 0.255785004 0.403782208
0.201002799 0.32218441 0.476812791
0.076187423 0.381624039 0.542188538
];
beta(1)=1;
beta(2)=1;
beta(3)=1;
beta(4)=1;
beta(5)=0.45;
beta(6)=0.35;
x1=X(:,1);
x2=X(:,2);
x3=X(:,3);
y=@(beta,X)beta(1)/(beta(2)*(pi/2)^(1/2))*exp((beta(3)*(x3-beta(5)).^2+beta(4)*(x1-beta(6)).^2)./beta(2)^2);
主程序如下:
y=[0.029114418
0.047518884
0.05892335
0.077327817
0.104732283
0.12113675
0.139541216
0.169945682
0.223360053
0.233194052
0.230028051
0.24586205
0.262696049
0.270530048
0.279364047
0.296198046
0.264832119
0.270100439
0.276368758
0.296637077
0.309905396
0.319173716
0.335442035
0.354710354
0.262048045
0.266815661
0.275583276
0.288350892
0.303118507
0.317886122
0.330653738
0.357421353
0.508638154
0.515434762
0.51623137
0.531027978
0.541824586
0.549621194
0.562417802
0.57321441
0.092238749
0.089985015
0.084731281
0.091477547
0.100223812
0.099970078
0.106716344
0.12846261
0.056867221
0.090869063
0.106870905
0.119872747
0.130874588
0.13887643
0.171878272
0.194880114
0.277725025
0.298843819
0.308962614
0.309081408
0.312200203
0.318318998
0.326437792
0.334556587
0.193122638
0.205416712
0.213710787
0.214004862
0.208298936
0.200593011
0.207887086
0.19418116
0.37063685
0.384010093
0.387383336
0.382756579
0.370129822
0.358503066
0.360876309
0.351249552
0.28988524
0.297163091
0.306440942
0.307718792
0.306996643
0.306274494
0.314552344
0.320830195
0.086052776
0.084954916
0.089857057
0.089759197
0.090661338
0.096563478
0.107465619
0.117367759
];
X=[0.83439788 0.049170736 0.116431385
0.83439788 0.049170736 0.116431385
0.643685573 0.138181349 0.218133078
0.475664857 0.216601181 0.307733962
0.313687287 0.314252942 0.372059771
0.193221192 0.348425311 0.458353497
0.113377277 0.378605746 0.508016977
0.877547538 0 0.122452462
0.674765579 0.096568907 0.228665514
0.497203845 0.181127428 0.321668727
0.340432788 0.255785004 0.403782208
0.201002799 0.32218441 0.476812791
0.076187423 0.381624039 0.542188538
];
beta0=[0.5,0.5,-2,-2,0.5,0.5];
[beta,r,j] = nlinfit(X,y,model,beta0);
(希望加上的限制条件是beta(3)<0;beta(4)<0;)但不知道怎么加,希望大家帮个忙,老板那里催了,救命。。。