 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 2 | 1/1 | 返回列表 |
| 查看: 2353 | 回復(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ù)此樓» 猜你喜歡
307求調(diào)劑
已經(jīng)有9人回復(fù)
-
288求調(diào)劑
已經(jīng)有15人回復(fù)
-
085700資源與環(huán)境308求調(diào)劑
已經(jīng)有5人回復(fù)
-
294求調(diào)劑材料與化工專碩
已經(jīng)有14人回復(fù)
-
0817 化學(xué)工程 299分求調(diào)劑 有科研經(jīng)歷 有二區(qū)文章
已經(jīng)有18人回復(fù)
-
265求調(diào)劑
已經(jīng)有3人回復(fù)
-
296求調(diào)劑
已經(jīng)有6人回復(fù)
-
0703化學(xué)調(diào)劑
已經(jīng)有13人回復(fù)
-
085600材料與化工 求調(diào)劑
已經(jīng)有15人回復(fù)
-
生物學(xué)調(diào)劑招人�。�！
已經(jīng)有4人回復(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ā)表 | |
|---|---|---|---|---|
|
[考研] 307求調(diào)劑 +9 | 冷笙123 2026-03-17 | 9/450 |
|
|---|---|---|---|---|
|
|
[考研] 296求調(diào)劑 +3 | www_q 2026-03-18 | 6/300 |
|
|
[考研] 317求調(diào)劑 +3 | 申子申申 2026-03-19 | 6/300 |
|
|
|
[考研] 求調(diào)劑 +3 | Mqqqqqq 2026-03-19 | 3/150 |
|
|
|
[考研] 材料080500調(diào)劑求收留 +4 | 一顆meteor 2026-03-13 | 4/200 |
|
|
|
[考研] 274求調(diào)劑 +6 | S.H1 2026-03-18 | 6/300 |
|
|
|
[考研] 295求調(diào)劑 +3 | 一志愿京區(qū)211 2026-03-18 | 5/250 |
|
|
|
[考研] 085601專碩，總分342求調(diào)劑，地區(qū)不限 +5 | share_joy 2026-03-16 | 5/250 |
|
|
|
[考研] 311求調(diào)劑 +6 | 26研0 2026-03-15 | 6/300 |
|
|
|
[考研] 一志愿西南交大，求調(diào)劑 +4 | 材化逐夢(mèng)人 2026-03-18 | 4/200 |
|
|
|
[考研] 304求調(diào)劑 +12 | 小熊joy 2026-03-14 | 13/650 |
|
|
|
[考研]
 0854，計(jì)算機(jī)類招收調(diào)劑
+3
|
胡辣湯放糖 2026-03-15 | 6/300 |
|
|
|
[考研] 278求調(diào)劑 +5 | 煙火先于春 2026-03-17 | 5/250 |
|
|
|
[考博] 26博士申請(qǐng) +3 | 1042136743 2026-03-17 | 3/150 |
|
|
|
[考研] 308求調(diào)劑 +4 | 是Lupa啊 2026-03-16 | 4/200 |
|
|
|
[考研] 0703化學(xué)調(diào)劑 290分有科研經(jīng)歷，論文在投 +7 | 膩膩gk 2026-03-14 | 7/350 |
|
|
|
[考研] 326求調(diào)劑 +3 | mlpqaz03 2026-03-15 | 3/150 |
|
|
|
[考研] 本科南京大學(xué)一志愿川大藥學(xué)327 +3 | 麥田耕者 2026-03-14 | 3/150 |
|
|
|
[考研] 297一志愿上交085600求調(diào)劑 +5 | 指尖八千里 2026-03-14 | 5/250 |
|
|
|
[考研] 一志愿哈工大材料324分求調(diào)劑 +5 | 閆旭東 2026-03-14 | 5/250 |
|