matlab条件限制曲线,带限制条件的参数拟合

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

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值