 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 5 | 1/1 | 返回列表 |
| 查看: 938 | 回復(fù): 3 | ||
| 當(dāng)前只顯示滿足指定條件的回帖，點(diǎn)擊這里查看本話題的所有回帖 | ||
夢(mèng)魘拾柒木蟲 (著名寫手)
|
[求助]
Matlab程序求助 已有1人參與
|
|
|
function multifractal(A) format long g L=length(A); i=1; modify=1; tmin=2; % 邊框間距，“※” tmax=10; ttmin=-10; ttmax=10; % 自定義 q 的范圍 for r=tmin:1:tmax c(i,1)=mod(L,r); i=i+1; end c'; % 計(jì)算不能被邊長(zhǎng)r整除的余數(shù) a=L-c'; % 計(jì)算并剔除掉不能被邊長(zhǎng)r整除的原始數(shù)據(jù) n=length(a); % 求解格網(wǎng)化邊長(zhǎng)的個(gè)數(shù)，即為 n TT=[]; j=1; r=tmin; % 自定義項(xiàng)，“※-2” for i=1:1:n % 即n=25-10+1，自定義的結(jié)果 B=A(1:a(i),1); U=reshape(B,r,length(B)/r); T=mean(U); T=T'.*r^3; TT(1: length(T),i)=[T]; modifying(modify,1)=length(T); r=r+1; if r>= tmin+n % 自定義項(xiàng)，“※-3” break; % 邊長(zhǎng)超過10+n，超過初始限制，則程序自動(dòng)終止 end modify=modify+1; end modifying; TT= nthroot(TT,1); % 或者不縮小 % TT; TT=nonzeros(TT); for cugb=1:1:n modifying_modifying(cugb,1)=sum(modifying(1:cugb)); % 有影響的新加卷 % & * % ￥ # @ �。� …… end modifying_modifying; j=1; % q 為任意數(shù)，這里取1到n，為n，與 k取值保持一致，q過大，計(jì)算機(jī)無(wú)法 %識(shí)別，默認(rèn)為無(wú)窮大，q過小，結(jié)果接近0，則意義不明確 for q=ttmin:ttmax %這里取 q=-10:1:10 for k=1:1:n X=TT(1:modifying_modifying(k,1),1).^q; if k>1 X=TT(modifying_modifying(k-1,1)+1:modifying_modifying(k,1),1).^q; end t=sum(X); XX(k,j)=[t]; % 這里用到兩個(gè)循環(huán)，即考慮到了冪函數(shù)，又需考慮求和 end j=j+1; end XX; % 得到質(zhì)量分配函數(shù)，Xq(ξ)， X=log(tmin:1:tmax); % X=log(tmax:-1:tmin); % 此系以前的自定義輸入結(jié)果，“※—4” Y=log(XX); % figure(1) % plot(X',Y,'o-k') % 至此，計(jì)算多重分形譜的第一步，分配函數(shù)構(gòu)建完畢 side_length= tmin:1:tmax; % 自定義網(wǎng)格邊長(zhǎng)，“※—5” side_length=side_length'; q=ttmin:ttmax; % q=-5:1:n-10，q=-10:1:10 m=1; [ha,hb]=size(XX); for i=1:1:hb % XX=XX’; % or not s=XX(:,i); % XX b=polyfit(log(side_length),log(s),1); % 在對(duì)數(shù)尺度下計(jì)算斜率 slope(m,1)=b(1,1); m=m+1; end slope; % 這里的Slope即為質(zhì)量指數(shù)，τ(q) N=polyfit(q', slope,1); plot(q', slope) % 此步是考察τ(q)-q 之間的關(guān)系， a=diff(slope)./diff(q'); % 第三步計(jì)算，diff函數(shù)求偏導(dǎo)確實(shí)少一列 q=q'; f_a=a.*q(1:end-1)-slope(1:end-1); % figure(2) % a=sort(a,'ascend'); plot(a,f_a,'o-k') xlabel('α','FontSize',12); ylabel('f(α)','FontSize',12); % polyfit(a,f_a,3) a=sort(a,'ascend'); % a+1 % f_a+1 運(yùn)行后出現(xiàn)??? Attempted to access modifying_modifying(0,1); index must be a positive integer or logical. Error in ==> multifractal at 47 X=TT(1:modifying_modifying(k-1,1),1).^q;哪位童鞋能幫我看看嘛 |
回復(fù)此樓» 猜你喜歡
伙伴們，祝我生日快樂吧
已經(jīng)有5人回復(fù)
-
274求調(diào)劑0856材料化工
已經(jīng)有10人回復(fù)
-
337一志愿華南理工0805材料求調(diào)劑
已經(jīng)有6人回復(fù)
-
（081700）化學(xué)工程與技術(shù)-298分求調(diào)劑
已經(jīng)有7人回復(fù)
-
材料與化工304求B區(qū)調(diào)劑
已經(jīng)有5人回復(fù)
-
NSFC申報(bào)書里申請(qǐng)人簡(jiǎn)歷中代表性論著還需要在申報(bào)書最后的附件里面再上傳一遍嗎
已經(jīng)有9人回復(fù)
-
一志愿北京化工大學(xué)材料與化工296分求調(diào)劑
已經(jīng)有17人回復(fù)
-
0703化學(xué)一志愿211 總分320求調(diào)劑
已經(jīng)有3人回復(fù)
-
0703求調(diào)劑
已經(jīng)有6人回復(fù)
-
一志愿天津大學(xué)材料與化工275求調(diào)劑
已經(jīng)有22人回復(fù)
【答案】應(yīng)助回帖
感謝參與，應(yīng)助指數(shù) +1
|
or k=1:1:n X=TT(1:modifying_modifying(k,1),1).^q; if k>1 X=TT(modifying_modifying(k-1,1)+1:modifying_modifying(k,1),1).^q; end 這個(gè)中的TT時(shí)二維數(shù)組嗎？1：modifying_modifying(k,1)是想說(shuō)從一到modifying_modifying(k,1)這個(gè)數(shù)，1^q,2^q,3^q......modifying_modifying(k,1)^q嗎？ |
3樓2019-05-29 20:59:31
1107040177
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 386.4
- 散金: 10
- 帖子: 110
- 在線: 7.5小時(shí)
- 蟲號(hào): 1888288
- 注冊(cè): 2012-07-11
- 專業(yè): 食品科學(xué)基礎(chǔ)
2樓2019-05-29 17:37:25
夢(mèng)魘拾柒
木蟲 (著名寫手)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 2225.7
- 散金: 6355
- 紅花: 4
- 帖子: 1027
- 在線: 240.3小時(shí)
- 蟲號(hào): 987900
- 注冊(cè): 2010-04-02
- 性別: GG
- 專業(yè): 地質(zhì)學(xué)
送紅花一朵 |
k=1:1:n X=TT(1:modifying_modifying(k,1),1).^q; if k>1 X=TT(modifying_modifying(k，1)+1:modifying_modifying(k,1),1).^q; 有點(diǎn)小錯(cuò)誤。應(yīng)為（k,1）,是一維數(shù)組。。兄臺(tái)能幫我調(diào)試一下，A為一維列向量即可。。 發(fā)自小木蟲Android客戶端 |
4樓2019-05-29 23:32:26
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] （081700）化學(xué)工程與技術(shù)-298分求調(diào)劑 +7 | 11啦啦啦 2026-03-11 | 7/350 |
|
|---|---|---|---|---|
|
|
[考研] 材料與化工304求B區(qū)調(diào)劑 +4 | 邱gl 2026-03-11 | 5/250 |
|
|
|
[考研] 341求調(diào)劑 +3 | 搗蛋豬豬 2026-03-11 | 3/150 |
|
|
|
[考研] 288求調(diào)劑 +13 | 王曉陽(yáng)- 2026-03-09 | 18/900 |
|
|
[考研] 化工0817調(diào)劑 +8 | 燦若星晨 2026-03-10 | 8/400 |
|
|
|
[考研] 298求調(diào)劑 +3 | Vv呀！ 2026-03-10 | 3/150 |
|
|
|
[考研] 327求調(diào)劑 +3 | Ffff03 2026-03-10 | 3/150 |
|
|
|
[考研] 313分生物學(xué)求調(diào)劑 +5 | Yyt楊1 2026-03-09 | 7/350 |
|
|
|
[考研] 2026考研求調(diào)劑-材料類-本科211一志愿985-初試301分 +10 | 蟲友233 2026-03-07 | 10/500 |
|
|
|
[考研] 085701環(huán)境工程295求調(diào)劑，一志愿A區(qū)211 +8 | 等春來(lái)， 2026-03-06 | 9/450 |
|
|
|
[碩博家園] 木蟲好像不熱鬧了，是不是？ +4 | 偏振片 2026-03-10 | 4/200 |
|
|
|
[基金申請(qǐng)] 面上項(xiàng)目還需要AI說(shuō)明嗎？ +3 | liyundong 2026-03-08 | 3/150 |
|
|
|
[考研]
 295復(fù)試調(diào)劑
+4
|
簡(jiǎn)木ChuFront 2026-03-09 | 4/200 |
|
|
|
[考研] 296求調(diào)劑 +4 | Xinyu Wu311 2026-03-09 | 4/200 |
|
|
|
[考研] 0856求調(diào)劑 +3 | squirtle11 2026-03-07 | 3/150 |
|
|
|
[考研] 290 材料與化工求調(diào)劑 +7 | Nebulala 2026-03-08 | 7/350 |
|
|
|
[考博] 26申博自薦 +9 | 烷基八氮- 2026-03-06 | 11/550 |
|
|
|
[考研] 材料專碩323求調(diào)劑 +5 | 李白26 2026-03-07 | 5/250 |
|
|
|
[考研] 求調(diào)劑推薦 +4 | 微辣不吃 2026-03-06 | 4/200 |
|
|
|
[考研] 求調(diào)劑 +4 | 呼呼？~+123456 2026-03-06 | 4/200 |
|