 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 查看: 1982 | 回復(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) 這個(gè)公式是把離散點(diǎn)擬合成橢圓的一個(gè)公式。有哪位同學(xué)能給我講一下這個(gè)公式是怎么來的，我不太理解。 |
回復(fù)此樓» 猜你喜歡
食品工程專碩一志愿中海洋309求調(diào)劑
已經(jīng)有8人回復(fù)
-
求調(diào)劑
已經(jīng)有7人回復(fù)
-
311（085601）求調(diào)劑
已經(jīng)有4人回復(fù)
-
305求調(diào)劑
已經(jīng)有4人回復(fù)
-
一志愿北化085600材料專碩275|有文章專利｜求調(diào)劑
已經(jīng)有7人回復(fù)
-
一志愿華理，數(shù)一英一285求A區(qū)調(diào)劑
已經(jīng)有12人回復(fù)
-
289求調(diào)劑
已經(jīng)有11人回復(fù)
-
生物學(xué)學(xué)碩，一志愿湖南大學(xué)，初試成績338
已經(jīng)有6人回復(fù)
-
321求調(diào)劑
已經(jīng)有7人回復(fù)
-
343求調(diào)劑
已經(jīng)有4人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 求大神解釋~~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ù)雜動(dòng)力學(xué)ode多參數(shù)擬合 求助@dingd
已經(jīng)有2人回復(fù)
- matlab
已經(jīng)有5人回復(fù)
- 注冊(cè)化工工程師2011年案例
已經(jīng)有43人回復(fù)
- 大家?guī)兔匆粋�(gè)matlab的小程序，看看哪里出錯(cuò)了
已經(jīng)有5人回復(fù)
- matlab問題
已經(jīng)有1人回復(fù)
- matlab問題
已經(jīng)有7人回復(fù)
- matlab程序出錯(cuò)調(diào)試
已經(jīng)有7人回復(fù)
- 一個(gè)MATLAB程序，運(yùn)行不出來
已經(jīng)有1人回復(fù)
- matlab程序問題，總是出錯(cuò)，麻煩大家看看問題在那里？
已經(jīng)有6人回復(fù)
- 利用Matlab進(jìn)行熱力學(xué)一致性檢驗(yàn)-面積法
已經(jīng)有4人回復(fù)
- 好玩的matlab的3D心形代碼（已全部驗(yàn)證）
已經(jīng)有4人回復(fù)
- 一個(gè)HF方法計(jì)算H2鍵長的matlab代碼（一百多行）
已經(jīng)有2人回復(fù)
- 【求助】安裝vasp出錯(cuò) make: *** [fftmpi_map.o] 錯(cuò)誤 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小時(shí)
- 蟲號(hào): 2601957
- 注冊(cè): 2013-08-17
- 性別: MM
- 專業(yè): 動(dòng)力學(xué)與控制
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
這應(yīng)該不是擬合橢圓的公式吧， 這是個(gè)匿名函數(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 你這個(gè)屬于第二種，是個(gè)匿名函數(shù) 輸入p和x就能出對(duì)應(yīng)函數(shù)值，不過x的格式是個(gè)橢圓 p 輸入4個(gè)系數(shù) x的第一列是當(dāng)作X，x的第二列當(dāng)作Y 方程就是 P1*X^2 + P2 *X*Y + P3*Y ^2 + P4 這個(gè)函數(shù)描述的是 這個(gè)東西 |
2樓2015-12-24 08:43:49
guo_zy_123
木蟲 (正式寫手)
- 應(yīng)助: 32 (小學(xué)生)
- 金幣: 4141
- 紅花: 11
- 帖子: 442
- 在線: 197.5小時(shí)
- 蟲號(hào): 2601957
- 注冊(cè): 2013-08-17
- 性別: MM
- 專業(yè): 動(dòng)力學(xué)與控制
3樓2015-12-24 08:47:06
getengqing
木蟲 (正式寫手)
- 應(yīng)助: 104 (高中生)
- 金幣: 3759.8
- 散金: 30
- 紅花: 5
- 帖子: 592
- 在線: 110.5小時(shí)
- 蟲號(hào): 2986702
- 注冊(cè): 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小時(shí)
- 蟲號(hào): 291104
- 注冊(cè): 2006-10-28
5樓2015-12-24 10:27:40
guo_zy_123
木蟲 (正式寫手)
- 應(yīng)助: 32 (小學(xué)生)
- 金幣: 4141
- 紅花: 11
- 帖子: 442
- 在線: 197.5小時(shí)
- 蟲號(hào): 2601957
- 注冊(cè): 2013-08-17
- 性別: MM
- 專業(yè): 動(dòng)力學(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];%請(qǐng)將用實(shí)際的數(shù)代換p1等 %輸入采集到的數(shù)據(jù) X = [x11,x12; x21,x22;x31,x32.....];%用采集到的數(shù)據(jù)代換 %求解F ans = F(p,X); %此時(shí)求解的 ans是個(gè)列向量，每個(gè)元素是采集到的點(diǎn)與對(duì)應(yīng)P中系數(shù)橢圓方程 %的偏差，這個(gè)ans單個(gè)元素都趨近與0時(shí)，這組系數(shù)p就是你要擬合的橢圓方程系數(shù) |
6樓2015-12-24 10:27:47
guo_zy_123
木蟲 (正式寫手)
- 應(yīng)助: 32 (小學(xué)生)
- 金幣: 4141
- 紅花: 11
- 帖子: 442
- 在線: 197.5小時(shí)
- 蟲號(hào): 2601957
- 注冊(cè): 2013-08-17
- 性別: MM
- 專業(yè): 動(dòng)力學(xué)與控制
7樓2015-12-24 10:29:16
guo_zy_123
木蟲 (正式寫手)
- 應(yīng)助: 32 (小學(xué)生)
- 金幣: 4141
- 紅花: 11
- 帖子: 442
- 在線: 197.5小時(shí)
- 蟲號(hào): 2601957
- 注冊(cè): 2013-08-17
- 性別: MM
- 專業(yè): 動(dòng)力學(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
|
%給你個(gè)更加詳細(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小時(shí)
- 蟲號(hào): 2455769
- 注冊(cè): 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]); 麻煩請(qǐng)問能不能設(shè)置擬合出來的橢圓長軸沿y方向，短軸沿x方向？ 還有這個(gè)公式F=@(p,x)p(1)*x(:,1).^2+p(2)*x(:,1).*x(:,2)+p(3)*x(:,2).^2+p(4) 哪個(gè)代表長軸哪個(gè)代表短軸呀~？小白一個(gè)跪謝大神 |
» 本帖附件資源列表
-
歡迎監(jiān)督和反饋：小木蟲僅提供交流平臺(tái)，不對(duì)該內(nèi)容負(fù)責(zé)。
本內(nèi)容由用戶自主發(fā)布，如果其內(nèi)容涉及到知識(shí)產(chǎn)權(quán)問題，其責(zé)任在于用戶本人，如對(duì)版權(quán)有異議，請(qǐng)聯(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小時(shí)
- 蟲號(hào): 2601957
- 注冊(cè): 2013-08-17
- 性別: MM
- 專業(yè): 動(dòng)力學(xué)與控制
|
去掉交叉項(xiàng)就是xy在坐標(biāo)軸上了，不過貌似擬合效果不好就是x*y那個(gè)項(xiàng)去掉，對(duì)應(yīng)p也少系數(shù)就行了 發(fā)自小木蟲IOS客戶端 |
10樓2015-12-26 16:29:50
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 ww汪紅 的主題更新
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 321求調(diào)劑 +6 | 璞玉~~ 2026-03-25 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 332求調(diào)劑 +3 | @MZB382400 2026-03-28 | 3/150 |
|
|
[考研] 070305高分子化學(xué)與物理 304分求調(diào)劑 +4 | c297914 2026-03-28 | 4/200 |
|
|
|
[考研] 070300求調(diào)劑306分 +4 | 26要上岸 2026-03-27 | 4/200 |
|
|
|
[考研] 266分，求材料冶金能源化工等調(diào)劑 +7 | 哇呼哼呼哼 2026-03-27 | 9/450 |
|
|
|
[考研] 085405 考的11408求各位老師帶走 +3 | Qiu學(xué)ing 2026-03-28 | 3/150 |
|
|
|
[考研] 333求調(diào)劑 +3 | question挽風(fēng) 2026-03-23 | 3/150 |
|
|
|
[考研] 286求調(diào)劑 +4 | lim0922 2026-03-26 | 4/200 |
|
|
|
[考研] 調(diào)劑求收留 +7 | 果然有我 2026-03-26 | 7/350 |
|
|
|
[考研] 325求調(diào)劑 +3 | Aoyijiang 2026-03-23 | 3/150 |
|
|
|
[考研] 329求調(diào)劑 +5 | 1() 2026-03-22 | 5/250 |
|
|
|
[考研]
 085602化學(xué)工程求調(diào)劑。
+4
|
平樂樂樂 2026-03-26 | 4/200 |
|
|
|
[考研] 材料考研求調(diào)劑 +3 | Dendel 2026-03-23 | 6/300 |
|
|
|
[考研] 一志愿北京化工大學(xué)材料與化工（085600）296求調(diào)劑 +9 | 稻妻小編 2026-03-26 | 9/450 |
|
|
|
[考研] 【2026考研調(diào)劑】制藥工程 284分 求相關(guān)專業(yè)調(diào)劑名額 +4 | 袁奐奐 2026-03-25 | 8/400 |
|
|
|
[考研] 318求調(diào)劑 +3 | plum李子 2026-03-23 | 3/150 |
|
|
|
[考研] 生物學(xué)學(xué)碩求調(diào)劑 +7 | 小羊睡著了? 2026-03-23 | 10/500 |
|
|
|
[考研] 接收2026碩士調(diào)劑(學(xué)碩+專碩) +4 | allen-yin 2026-03-23 | 6/300 |
|
|
|
[考研]
一志愿東華大學(xué)化學(xué)070300，求調(diào)劑
+7
|
2117205181 2026-03-21 | 8/400 |
|
|
|
[考研] 280分求調(diào)劑 一志愿085802 +4 | PUMPT 2026-03-22 | 7/350 |
|