1樓:電子
%by dynamic
%see also
%contact me matlabsky@gmail.***%2009.2.
%程式編好了,但是初值真的很難選,我到現在還沒有選這個一個合適的初值,你根據實際選擇一個吧,否則程式是執行不起來的
你的函式exp,初值稍微選的不好就會導致0或者∞function mm
clcx=[2030 13 320 1
1850 12 320 1
1850 12 320 1
1970 13 320 1
2090 14 320 1
2210 15 320 1
3010 20 320 1
2540 17 320 1
2710 18 620 2];
p=[1
1 11 1
0.90990991
0.09009009
0.544303797
0.455696203];
b0=[-.005 2 -0.05 8];
b = nlinfit(x,p,@fun,b0)function yhat=fun(bb,x)alpha=x(:,1);
beta=x(:,2);
gamma=x(:,3);
delta=x(:,4);
a=bb(1);
b=bb(2);
c=bb(3);
d=bb(4);
num=exp(a*alpha+b*beta+c*gamma+d*delta);
den=sum(num);
yhat=num./den;
如何用matlab非線性迴歸分析
2樓:兄弟連教育
^%把y=x^a兩邊取bai對數,就有log(y)=a*log(x),如果將log(y)看成
du是log(x)的函zhi數,那麼它們是dao線性的.這時回是方程的個數是10000個,而未知數答是1個(a),這可以用matlab的矩陣除法求解矛盾方程(用的是最小二乘法)
%下面是一個例子,我以a=1.45產生一套10000個的(x,y)序列,當然我加上了一些隨機擾動.(你可以試一下,不加的話,肯定得到的是1.45)
>> a=1.45;
>> x=linspace(1,10,10000)'; %注意這裡你要轉成列向量求解
>> y=x.^a+0.5*rand(10000,1);
>> aa=log(x)\log(y)
如何利用matlab擬合非線性迴歸函式的初始引數
3樓:匿名使用者
推薦你使用
一個函式nlinfit,我簡單給你介紹一下使用方法,以你的模型2為例: 第一步:你需要建立一個function檔案,名字隨便,這裡我們命名為"hougen",在這個檔案內,你要把模型二描述清楚,首先我們假定beta=[a;b;c;d;e;f;g];x=[x,y,m];
matlab 曲線擬合求函式的最佳引數時怎麼確定引數的初始值x=[2:16];
y=[6.42 8.2 9.
58 9.5 9.7 10 9.
93 9.99 10.49 10.
59 10.6 10.8 10.
6 10.9 10.76];
plot(x,y,'*');
% b0=;%待定引數的初始值
fun=inline('b(1)*(1-b(2)*exp(b(3)*x))','b','x');
[b,r,j]=nlinfit(x,y,fun);
br=sum(r.^2)%誤差的平方和
在spss軟體中非線性迴歸的分析的散點圖中可以縱設定截距不
spss統計軟體可以用來做許多資料分析,迴歸分析就是其中之一。迴歸分析就是探索兩種及其以上變數之間的關係,運用十分廣泛,按照自變數和因變數之間的函式關係型別可以分為線性迴歸分析和非線性迴歸分析。迴歸分析不僅可以分析資料,更可以用來 一些資料的發展情況,從而應用非常廣泛。多因素迴歸分析步驟如下 1 開...
MATLAB非線性約束問題,急求
你這個問題的難點在於如何表達目標函式和約束函式,必須把已知量都給出,否則,這兩個函式寫不出來。請把已知量都寫出來,或用資料檔案發過來。matlab遺傳演算法優化問題中,非線性等式約束怎麼新增啊?建立約束條件函式,把非線性的等式約束條件新增加在 c,ceq 中。如function c,ceq nonc...
我用matlab求解約束非線性規劃問題,但是老是出現下面
看到數字了麼,這裡的數字,代表滾動速度 這裡輸入文字 你可以試著修改,回具體需要什答麼樣子的可以根據實際來弄。這個是個最簡單的能自己滾動的 了。根本不需要js運算就能成功的。求助一個用matlab求解非線性規劃的問題,不勝感激!1.把 x,fval fmincon fun x0,a,b,aeq,be...