1樓:網友
y=1+a*exp(b*(x1)^2+c*(x2)^2)
可以變成。log(y-1)=log(a)+b*(x1)^2+c*(x2)^2
這樣的話,就很好擬合了。
%先以一定的引數產生一系列的y,再加一些噪音。
以含有噪音的y來反推出a,b,c
a=5;b=;
c=;n=20; %點數。
x1=linspace(-1,1,n);
x2=linspace(,,n);
y=1+a*exp(b*x1.^2+c*x2.^2); 生成資料。
y=y+;加了一些噪音。
下面開始反推abc
a=[ones(size(x1(:)x1(:)2,x2(:)2];
b=log(y(:)1);
para=a\b; %求出來的log(a),b,c組成的列向量。
aa=exp(para(1));
bb=para(2);
cc=para(3);
aa,bb,cc
yy=1+aa*exp(bb*x1.^2+cc*x2.^2);
得到的結果。aa =
bb =cc =
與我們給定的值差不多。
=誰說我上面不禪並是最小二乘法,matlab求解超定方程組的時候,都是用非線性的。
你始果一定要用非線性擬合的話,可以擬合,但是擬合跟你的初值關係太大了。
還是給你乙個非線性擬合的**吧!
a0=5;b0=;
c0=;n=20; %點數。
x1=linspace(-1,1,n);
x2=linspace(,,n);
y=1+a0*exp(b0*x1.^2+c0*x2.^2); 生成資料。
y=y+;加了一廳塌些噪音。
下面開始反推abc
fun=@(p,x)1+p(1)*exp(p(2)*x(:,1).^2+p(3)*x(:,2).^2);
para,res]=lsqcurvefit(fun,[1,,1],[x1(:)x2(:)y(:)
aa=para(1);
bb=para(2);
cc=para(3);
aa,bb,cc
yy=1+aa*exp(bb*x1.^2+cc*x2.^2);
2樓:網友
使用cftool,輸入y,x1,x2,判稿x3資料,在互動介面,自己編公式y=1+a*exp(b*(x1)^2+c*(x2)^2),擬合後逗簡會輸出a,b,c等引數,a,b,c的誤差,殘山衝褲差等。
matlab如何做線性擬合
3樓:信必鑫服務平臺
方法一。1、最常用的是多項式擬合,採用polyfit函式,在命令視窗輸入自變數x和因變數y。
2、以二次多項式擬合為例,輸入p=polyfit(x,y,2),如果想擬合更高次的多項式,更換括號內數字即可。
方法二。1、在matlab自帶的曲線擬合工具包上方工具欄選取apps,點選curve fitting。在上方工具欄選取apps,點選curve fitting
2、輸入自變數x和因變數y。
3、選擇擬合方式,有多項式擬合polynomial,高斯擬合gaussian,冪指數擬合power等等,本次以多項式擬合為例。
4、通過資料計算,可以獲得曲線引數(曲線函式中的各項係數),從而實現曲線擬合。
matlab中多元擬合的函式是什麼?
4樓:long1209龍是我
1、因為這個函式的形式比較特殊對其兩邊取對數後得到 log(y)=log(a)+b*log(x1)+c*log(x2)
2、於是立即就轉換為了線性擬合 [1 log(x1) log(x2)]*log(a) b c]'=log(y)
3、於是[log(a) b c]'=[1 log(x1) log(x2)]\log(y) 這樣就可以得到a,b,c了。不需要進行非線性擬合。
4、我想擬合出乙個公式 y=f(t1,t2),當有任意t1、t2值時,可計算出縱座標。
matlab擬合曲線的方法有幾種
matlab最常bai用擬合曲線的函 du數有 1 regress 和 zhipolyfit 適用於線性函式dao的擬合2 lsqcurvefit 和nlinfit 適用於較復內雜的非容線性函式的擬合 3 fit 適用於最常用函式的擬合,如指數函式,傅立葉函式 高斯函式 多項式函式 冪函式 有理函式...
績效考核的具體方法,績效考核方法
企業進行績效管理如果處理不當會出現比較大的問題,盲目的制定並不一定適合企業實際情況的績效考核制度有時會適得其反。kpi績效管理不錯,但要選好關鍵性績效指標。最好選擇專業機構諮詢一下,先進行企業診斷,再選擇合適的制度方案,只有合適的才是最好的。這方面美華企業管理諮詢機構做的不錯。績效考核方法 績效考核...
煮麵條具體方法,煮麵條的方法和步驟
不管手工面還是掛麵,表面都有一層薄薄的乾粉。在遇到滾燙的開水時會迅速糊化,讓麵條變得黏糊糊的,就不爽滑了,還容易粘連 起坨。開水下鍋麵條受熱過快,表面已經糊化了,但裡面還有白芯,根本沒煮熟,就要延長煮麵的時間,導致麵條發黏 發軟,沒有一點口感。麵館老闆說,想要煮麵條好吃,萬不可開水下鍋。 煮麵條不能...