2016-10-17 25 views
0

私は、生態系CO2交換(NEE)を予測するために、3つの独立変数(光、温度、蒸気圧不足(VPD))を使用して非線形モデルに適合させようとしています。Matlab:非線形回帰、2基準

私はnlinfit機能を使用する方法を知っているが、私の問題は、私は2つの基準を使用したいということです: -

1. if VPD < 1.3 
NEE = (Param(1).*Param(2).*Ind_var(:,1))./(Param(1).*Ind_var(:,1)+Param(2)) + Param(3).*(1.6324.^((Ind_var(:,2)-18)./10)); 
2. if VPD >= 1.3 
NEE = (Param(1).*(Param(2).*exp(-Param(4).*(Ind_var(:,3)-1.3))).*Ind_var(:,1))./(Param(1).*Ind_var(:,1)+(Param(2).*exp(-Param(4).*(Ind_var(:,3)-1.3))))+Param(3).*(1.6324.^((Ind_var(:,2)-18)./10)); 

基本的には、独立変数VPD(蒸気圧赤字は)1.3以下であれば、私は私のParam(4) = 0を強制したいです。

しかし、私はそれを行う方法がわかりません。

お手伝いできますか? ありがとう、 Alexis

答えて

0

Param(4)Param(4)*(VPD<1.3)に置き換えることができます。 条件式は、真なら1に、Matlabでは偽なら0に評価されます。

+0

ありがとうございます!それは完全に動作します。 –