 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 2354 | 回復(fù): 1 | ||
天微涼61鐵蟲 (初入文壇)
|
[求助]
matlab求函數(shù)所有極小值點(diǎn)
|
|
function prl() %考慮第二種情況widely spaced islands a=sort(80*rand(1,5)); %生成一組隨機(jī)數(shù)然后從小到大排序 v(1)=a(1)+0.5*(a(2)-a(1)); %由Voronoi polygon算第一個(gè)island的體積 for i=2:4, v(i)=(a(i+1)-a(i-1))*0.5; %由Voronoi polygon算2到n-1個(gè)位置體積 end v(5)=80-a(5)+0.5*(a(5)-a(4)); %最后一個(gè)island的體積 y=[zeros(1,length(a));v]; subplot(2,1,1); plot([a;a],y,'k'); axis([0,80,0,100]); title('layer 1'); syms x f=((x^2+1)^(-3/2))*(1-3/(x^2+1)); for i=1:5, f=f+(((x-a(i))^2+1)^(-3/2))*(1-3/((x-a(i))^2+1)); %由張力公式疊加算第二層的島的位置 end h=diff(f,x,1); %對(duì)總張力公式求一階導(dǎo) g=diff(f,x,2); %求二階導(dǎo) h=inline(vectorize(f),'x'); options=optimset('Display','off'); k=1 q=zeros(1,80); for j=0:300 q(k)=fzero(h,j,options); %找出一階導(dǎo)數(shù)的零點(diǎn) k=k+1; end p=zeros(1,80); p(1)=q(1); n=1; for m=1:k % if q(m)~=p(n) if g(q(m))>0 %判斷二階導(dǎo)數(shù)是否大于零 n=n+1; p(n)=q(m); %得到極值點(diǎn) end %end end v=zeros(1,80); %求第二層體積 v(1)=p(1)+0.5*(p(2)-p(1)); for i=2:m-1, v(i)=(p(i+1)-p(i-1))*0.5; end v(length(p))=80-p(length(p))+0.5*(p(length(p))-p(length(p)-1)); subplot(2,1,2) plot([p;p],[zeros(1,length(p));v],'k') axis([0,80,0,100]) title('layer 2'); 請(qǐng)問哪錯(cuò)了，求二至n層核的位置及體積，核的位置為f極小值點(diǎn),f為張力公式，根據(jù)核的位置疊加求極小值，極小值就是下一層核位置，最后根據(jù)泰森邊形球體積，求大神指導(dǎo)，小弟感激不盡[ Last edited by 天微涼61 on 2013-4-29 at 10:47 ] |
回復(fù)此樓» 猜你喜歡
求調(diào)劑
已經(jīng)有3人回復(fù)
-
265求調(diào)劑
已經(jīng)有4人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿吉林大學(xué)材料學(xué)碩321求調(diào)劑
已經(jīng)有12人回復(fù)
-
286分人工智能專業(yè)請(qǐng)求調(diào)劑愿意跨考！
已經(jīng)有3人回復(fù)
-
329求調(diào)劑
已經(jīng)有5人回復(fù)
-
申請(qǐng)回稿延期一個(gè)月，編輯同意了。但系統(tǒng)上的時(shí)間沒變，給編輯又寫郵件了，沒回復(fù)
已經(jīng)有4人回復(fù)
-
材料學(xué)碩318求調(diào)劑
已經(jīng)有5人回復(fù)
-
一志愿中國(guó)海洋大學(xué)，生物學(xué)，301分，求調(diào)劑
已經(jīng)有6人回復(fù)
-
081700化工學(xué)碩調(diào)劑
已經(jīng)有3人回復(fù)
» 本主題相關(guān)價(jià)值貼推薦，對(duì)您同樣有幫助:
- 用matlab里的nlinfit如何擬合分段函數(shù)
已經(jīng)有8人回復(fù)
- Matlab中如何自己定義一個(gè)新函數(shù)，然后計(jì)算時(shí)可以調(diào)用？
已經(jīng)有12人回復(fù)
- 求matlab畫一個(gè)函數(shù)的圖
已經(jīng)有6人回復(fù)
- matlab的函數(shù)調(diào)用問題
已經(jīng)有3人回復(fù)
- MATLAB中產(chǎn)生高斯白噪聲，涉及到awgn和wgn函數(shù)
已經(jīng)有3人回復(fù)
- 問個(gè)關(guān)于matlab里find函數(shù)的問題！
已經(jīng)有4人回復(fù)
- 【分享】Matlab常用函數(shù)及常用數(shù)學(xué)函數(shù)【已搜索無重復(fù)】
已經(jīng)有80人回復(fù)
- matlab做函數(shù)圖象的問題
已經(jīng)有6人回復(fù)
- 【分享】MATLAB常用函數(shù)參考
已經(jīng)有36人回復(fù)
天微涼61
鐵蟲 (初入文壇)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 13
- 帖子: 3
- 在線: 8.5小時(shí)
- 蟲號(hào): 2426567
- 注冊(cè): 2013-04-19
- 專業(yè): 電路與系統(tǒng)
2樓2013-05-19 12:35:13
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 天微涼61 的主題更新
| 2 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 296求調(diào)劑 +3 | www_q 2026-03-18 | 6/300 |
|
|---|---|---|---|---|
|
|
[考研] 304求調(diào)劑 +3 | 曼殊2266 2026-03-18 | 3/150 |
|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
|
[考研] 本科鄭州大學(xué)物理學(xué)院，一志愿華科070200學(xué)碩，346求調(diào)劑 +4 | 我不是一根蔥 2026-03-18 | 4/200 |
|
|
|
[考研] 材料專碩英一數(shù)二306 +5 | z1z2z3879 2026-03-18 | 5/250 |
|
|
|
[考研] 354求調(diào)劑 +4 | Tyoumou 2026-03-18 | 7/350 |
|
|
|
[考研] 328求調(diào)劑，英語六級(jí)551，有科研經(jīng)歷 +3 | 生物工程調(diào)劑 2026-03-17 | 7/350 |
|
|
|
[考研] 一志愿武理材料305分求調(diào)劑 +5 | 想上岸的鯉魚 2026-03-18 | 6/300 |
|
|
|
[考研] 26調(diào)劑/材料/英一數(shù)二/總分289/已過A區(qū)線 +7 | 步川酷紫123 2026-03-13 | 7/350 |
|
|
|
[考研] 材料專碩306英一數(shù)二 +10 | z1z2z3879 2026-03-16 | 13/650 |
|
|
|
[考研] 收復(fù)試調(diào)劑生 +4 | 雨后秋荷 2026-03-18 | 4/200 |
|
|
|
[考研] 生物學(xué)071000 329分求調(diào)劑 +3 | 我愛生物生物愛?/a> 2026-03-17 | 3/150 |
|
|
|
[考研] 265求調(diào)劑 +3 | 梁梁校校 2026-03-17 | 3/150 |
|
|
|
[考研] 301求調(diào)劑 +9 | yy要上岸呀 2026-03-17 | 9/450 |
|
|
|
[考研] 326求調(diào)劑 +5 | 上岸的小葡 2026-03-15 | 6/300 |
|
|
|
[考研] 278求調(diào)劑 +3 | Yy7400 2026-03-13 | 3/150 |
|
|
|
[考研] 一志愿，福州大學(xué)材料專碩339分求調(diào)劑 +3 | 木子momo青爭(zhēng) 2026-03-15 | 3/150 |
|
|
|
[考研] 070303 總分349求調(diào)劑 +3 | LJY9966 2026-03-15 | 5/250 |
|
|
|
[考研] 288求調(diào)劑 +4 | 奇點(diǎn)0314 2026-03-14 | 4/200 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|