 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1981 | 回復(fù): 11 | ||
ww汪紅新蟲 (小有名氣)
|
[求助]
F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4) 已有3人參與
|
|
F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4) 這個公式是把離散點(diǎn)擬合成橢圓的一個公式。有哪位同學(xué)能給我講一下這個公式是怎么來的，我不太理解。 |
回復(fù)此樓» 猜你喜歡
生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338
已經(jīng)有5人回復(fù)
-
332求調(diào)劑
已經(jīng)有3人回復(fù)
-
343求調(diào)劑
已經(jīng)有3人回復(fù)
-
學(xué)碩274求調(diào)劑
已經(jīng)有8人回復(fù)
-
壓國家一區(qū)線，求導(dǎo)師收留，有恩必謝！
已經(jīng)有7人回復(fù)
-
081200-11408-276學(xué)碩求調(diào)劑
已經(jīng)有4人回復(fù)
-
22408 359分調(diào)劑
已經(jīng)有4人回復(fù)
-
本科新能源科學(xué)與工程，一志愿華理能動285求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿華北電力大學(xué)能動專碩，293，求調(diào)劑
已經(jīng)有5人回復(fù)
-
調(diào)劑310
已經(jīng)有7人回復(fù)
» 本主題相關(guān)價值貼推薦，對您同樣有幫助:
- 求大神解釋~~F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4);
已經(jīng)有6人回復(fù)
- ??? Attempted to access C(2); index out of bounds because numel(C)=1.
已經(jīng)有3人回復(fù)
- 1stopt進(jìn)行復(fù)雜動力學(xué)ode多參數(shù)擬合 求助@dingd
已經(jīng)有2人回復(fù)
- matlab
已經(jīng)有5人回復(fù)
- 注冊化工工程師2011年案例
已經(jīng)有43人回復(fù)
- 大家?guī)兔匆粋�matlab的小程序，看看哪里出錯了
已經(jīng)有5人回復(fù)
- matlab問題
已經(jīng)有1人回復(fù)
- matlab問題
已經(jīng)有7人回復(fù)
- matlab程序出錯調(diào)試
已經(jīng)有7人回復(fù)
- 一個MATLAB程序，運(yùn)行不出來
已經(jīng)有1人回復(fù)
- matlab程序問題，總是出錯，麻煩大家看看問題在那里？
已經(jīng)有6人回復(fù)
- 利用Matlab進(jìn)行熱力學(xué)一致性檢驗(yàn)-面積法
已經(jīng)有4人回復(fù)
- 好玩的matlab的3D心形代碼（已全部驗(yàn)證）
已經(jīng)有4人回復(fù)
- 一個HF方法計(jì)算H2鍵長的matlab代碼（一百多行）
已經(jīng)有2人回復(fù)
- 【求助】安裝vasp出錯 make: *** [fftmpi_map.o] 錯誤 1【已解決】
已經(jīng)有6人回復(fù)
- 【資源】Runge-Kutta解ODE的一些程序及應(yīng)用（from MathWorks）
已經(jīng)有3人回復(fù)
- 【求助】求高手幫看一下程序（泡點(diǎn)溫度計(jì)算）
已經(jīng)有2人回復(fù)
- 【分享】關(guān)于Boltzmann 的計(jì)算代碼（MATLAB）
已經(jīng)有9人回復(fù)
guo_zy_123
木蟲 (正式寫手)
- 應(yīng)助: 32 (小學(xué)生)
- 金幣: 4141
- 紅花: 11
- 帖子: 442
- 在線: 197.5小時
- 蟲號: 2601957
- 注冊: 2013-08-17
- 性別: MM
- 專業(yè): 動力學(xué)與控制
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
這應(yīng)該不是擬合橢圓的公式吧， 這是個匿名函數(shù)啊 matlab中的函數(shù)寫法 可以有三種方法創(chuàng)建： �。�1）inline函數(shù)： f=inline('3*x+1','x') 可以得到：f= Inline function : f(x)=3*x+1 輸入t=0:3; f(t) ans = 1 4 7 10 �。�2）匿名函數(shù)： f=@(x)3*x+1 可以得到：f= @(x)3*x+1 輸入t=0:3; f(t) ans = 1 4 7 10 �。�3）創(chuàng)建M-函數(shù) 新建m文件，輸入： function f=equation(x) f=3*x+1; 保存m文件到工作文件夾； 調(diào)用： 輸入t=0:3; f=equation(t) ans = 1 4 7 10 你這個屬于第二種，是個匿名函數(shù) 輸入p和x就能出對應(yīng)函數(shù)值，不過x的格式是個橢圓 p 輸入4個系數(shù) x的第一列是當(dāng)作X，x的第二列當(dāng)作Y 方程就是 P1*X^2 + P2 *X*Y + P3*Y ^2 + P4 這個函數(shù)描述的是 這個東西 |
2樓2015-12-24 08:43:49
guo_zy_123
木蟲 (正式寫手)
- 應(yīng)助: 32 (小學(xué)生)
- 金幣: 4141
- 紅花: 11
- 帖子: 442
- 在線: 197.5小時
- 蟲號: 2601957
- 注冊: 2013-08-17
- 性別: MM
- 專業(yè): 動力學(xué)與控制
3樓2015-12-24 08:47:06
getengqing
木蟲 (正式寫手)
- 應(yīng)助: 104 (高中生)
- 金幣: 3759.8
- 散金: 30
- 紅花: 5
- 帖子: 592
- 在線: 110.5小時
- 蟲號: 2986702
- 注冊: 2014-02-23
- 性別: GG
- 專業(yè): 自然語言理解與機(jī)器翻譯
4樓2015-12-24 09:00:25
dingd
鐵桿木蟲 (職業(yè)作家)
- 應(yīng)助: 1641 (講師)
- 金幣: 15037.3
- 散金: 101
- 紅花: 234
- 帖子: 3410
- 在線: 1223.7小時
- 蟲號: 291104
- 注冊: 2006-10-28
5樓2015-12-24 10:27:40
guo_zy_123
木蟲 (正式寫手)
- 應(yīng)助: 32 (小學(xué)生)
- 金幣: 4141
- 紅花: 11
- 帖子: 442
- 在線: 197.5小時
- 蟲號: 2601957
- 注冊: 2013-08-17
- 性別: MM
- 專業(yè): 動力學(xué)與控制
|
按樓主的最早的描述方程 matlab輸入應(yīng)該是 %匿名函數(shù) F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4)； %輸入橢圓系數(shù) p = [p1,p2,p3,p4];%請將用實(shí)際的數(shù)代換p1等 %輸入采集到的數(shù)據(jù) X = [x11,x12; x21,x22;x31,x32.....];%用采集到的數(shù)據(jù)代換 %求解F ans = F(p,X); %此時求解的 ans是個列向量，每個元素是采集到的點(diǎn)與對應(yīng)P中系數(shù)橢圓方程 %的偏差，這個ans單個元素都趨近與0時，這組系數(shù)p就是你要擬合的橢圓方程系數(shù) |
6樓2015-12-24 10:27:47
guo_zy_123
木蟲 (正式寫手)
- 應(yīng)助: 32 (小學(xué)生)
- 金幣: 4141
- 紅花: 11
- 帖子: 442
- 在線: 197.5小時
- 蟲號: 2601957
- 注冊: 2013-08-17
- 性別: MM
- 專業(yè): 動力學(xué)與控制
7樓2015-12-24 10:29:16
guo_zy_123
木蟲 (正式寫手)
- 應(yīng)助: 32 (小學(xué)生)
- 金幣: 4141
- 紅花: 11
- 帖子: 442
- 在線: 197.5小時
- 蟲號: 2601957
- 注冊: 2013-08-17
- 性別: MM
- 專業(yè): 動力學(xué)與控制
【答案】應(yīng)助回帖
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
xzhdty: 金幣+10, 代樓主發(fā)放金幣 2015-12-24 20:41:09
ww汪紅: 金幣+40, ★★★★★最佳答案 2015-12-26 15:58:23
xzhdty: 金幣+10, 代樓主發(fā)放金幣 2015-12-24 20:41:09
ww汪紅: 金幣+40, ★★★★★最佳答案 2015-12-26 15:58:23
|
%給你個更加詳細(xì)橢圓函數(shù)的擬合 %復(fù)制到matlab運(yùn)行 %記得給金幣 %記得給金幣 %記得給金幣（重要的事情說三遍） clear; clc; close all; % 設(shè)出圓錐曲線方程 F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4)*x(:,1)+p(5)*x(:,2)+p(6); % 離散數(shù)據(jù)點(diǎn) x=[1.7729 1.9228 1.7338 1.9072 2.0539 1.6137 2.0656 1.6412 1.8611 1.48765 1.9005 1.4971 2.0732 1.6546 1.8338 1.9405 1.9375 1.5104 1.6878 1.5177 1.7031 1.5097 1.9577 1.5201 1.9872 1.5437 2.0341 1.5805 2.0723 1.6546 2.0681 1.8284 2.0557 1.8483 2.0491 1.5651 ]; p0=[1 1 1 1 1 1]; warning off % 擬合系數(shù)，最小二乘方法 p=nlinfit(x,zeros(size(x,1),1),F,p0); plot(x(:,1),x(:,2),'ro'); hold on; xmin=min(x(:,1)); xmax=max(x(:,1)); ymin=min(x(:,2)); ymax=max(x(:,2)); % 作圖 ezplot(@(x,y)F(p,[x,y]),[-1+xmin,1+xmax,-1+ymin,1+ymax]); title('曲線擬合'); legend('樣本點(diǎn)','擬合曲線') |
8樓2015-12-24 10:37:07
ww汪紅
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 141
- 散金: 5
- 帖子: 160
- 在線: 43.6小時
- 蟲號: 2455769
- 注冊: 2013-05-09
- 專業(yè): 凝聚態(tài)物性 II ：電子結(jié)構(gòu)
|
F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4); % 離散數(shù)據(jù)點(diǎn) x=load('G:\Users\Daniel\Desktop\Total.txt'); p0=[1 1 1 1]; warning off % 擬合系數(shù),最小二乘方法 p=nlinfit(x,zeros(size(x,1),1),F,p0); plot(x(:,1),x(:,2),'r.'); axis square; hold on; xmin=min(x(:,1)); xmax=max(x(:,1)); ymin=min(x(:,2)); ymax=max(x(:,2)); % 作圖 ezplot(@(x,y)F(p,[x,y]),[-1+xmin,1+xmax,-1+ymin,1+ymax]); title('曲線擬合'); legend('樣本點(diǎn)','擬合曲線') set(gca,'xtick',[-1:0.2:1]); set(gca,'ytick',[-1:0.2:1]); xlim([-1 1]); ylim([-1 1]); line([-1 1],[0 0]); line([0 0],[-1 1]); 麻煩請問能不能設(shè)置擬合出來的橢圓長軸沿y方向，短軸沿x方向？ 還有這個公式F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4) 哪個代表長軸哪個代表短軸呀~？小白一個跪謝大神 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺，不對該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對版權(quán)有異議，請聯(lián)系郵箱：xiaomuchong@tal.com - 附件 1 : Total.txt
2015-12-26 15:52:15, 1.47 K
9樓2015-12-26 15:54:06
guo_zy_123
木蟲 (正式寫手)
- 應(yīng)助: 32 (小學(xué)生)
- 金幣: 4141
- 紅花: 11
- 帖子: 442
- 在線: 197.5小時
- 蟲號: 2601957
- 注冊: 2013-08-17
- 性別: MM
- 專業(yè): 動力學(xué)與控制
|
去掉交叉項(xiàng)就是xy在坐標(biāo)軸上了，不過貌似擬合效果不好就是x*y那個項(xiàng)去掉，對應(yīng)p也少系數(shù)就行了 發(fā)自小木蟲IOS客戶端 |
10樓2015-12-26 16:29:50
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ww汪紅 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 22408 359分調(diào)劑 +3 | Qshers 2026-03-27 | 4/200 |
|
|---|---|---|---|---|
|
|
[考研] 304求調(diào)劑 +6 | 曼殊2266 2026-03-27 | 6/300 |
|
|
[考研] 一志愿太原理工安全工程300分，求調(diào)劑 +4 | 0857求調(diào)劑. 2026-03-24 | 5/250 |
|
|
|
[考研] 277跪求調(diào)劑 +5 | 1915668 2026-03-27 | 9/450 |
|
|
|
[考研]
 材料與化工考研調(diào)劑
+17
|
孅華 2026-03-22 | 17/850 |
|
|
|
[考研] 0703化學(xué)/290求調(diào)劑/本科經(jīng)歷豐富/工科也可 +9 | 丹青奶蓋 2026-03-26 | 10/500 |
|
|
|
[考研] 352分 化工與材料 +5 | 海納百川Ly 2026-03-27 | 5/250 |
|
|
|
[考研] 一志愿211院校 344分 東北農(nóng)業(yè)大學(xué)生物學(xué)學(xué)碩，求調(diào)劑 +5 | 丶風(fēng)雪夜歸人丶 2026-03-26 | 8/400 |
|
|
|
[考研] 274求調(diào)劑 +17 | 顧九笙要謙虛 2026-03-24 | 23/1150 |
|
|
|
[考研] 334求調(diào)劑 +3 | 雨清天晴 2026-03-21 | 3/150 |
|
|
|
[考研] 324求調(diào)劑 +8 | hanamiko 2026-03-26 | 10/500 |
|
|
|
[考研]
材料調(diào)劑
 5+4
|
想要一壺桃花水 2026-03-25 | 10/500 |
|
|
|
[考研] 291 求調(diào)劑 +7 | 化工2026屆畢業(yè)?/a> 2026-03-21 | 8/400 |
|
|
|
[考研] 309求調(diào)劑 +4 | gajsj 2026-03-25 | 5/250 |
|
|
|
[考研] 網(wǎng)絡(luò)空間安全0839招調(diào)劑 +4 | w320357296 2026-03-25 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +3 | 李李不服輸 2026-03-25 | 3/150 |
|
|
|
[考研] 材料專碩331求調(diào)劑 +4 | 鮮當(dāng)牛 2026-03-24 | 4/200 |
|
|
|
[考研] 307求調(diào)劑 +3 | 余意卿 2026-03-21 | 6/300 |
|
|
|
[考研] 336化工調(diào)劑 +4 | 王大坦1 2026-03-23 | 5/250 |
|
|
|
[考研] 石河子大學(xué)（211、雙一流）碩博研究生長期招生公告 +3 | 李子目 2026-03-22 | 3/150 |
|