 |
 |
 |
 |
 |
 |
||
 |
 |
 |
 |
 |
| 4 | 1/1 | 返回列表 |
| 查看: 930 | 回復(fù): 3 | ||
夢魘拾柒木蟲 (著名寫手)
|
[求助]
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'; % 計算不能被邊長r整除的余數(shù) a=L-c'; % 計算并剔除掉不能被邊長r整除的原始數(shù)據(jù) n=length(a); % 求解格網(wǎng)化邊長的個數(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; % 邊長超過10+n，超過初始限制，則程序自動終止 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ī)無法 %識別，默認(rèn)為無窮大，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]; % 這里用到兩個循環(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') % 至此，計算多重分形譜的第一步，分配函數(shù)構(gòu)建完畢 side_length= tmin:1:tmax; % 自定義網(wǎ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); % 在對數(shù)尺度下計算斜率 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'); % 第三步計算，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ù)此樓» 猜你喜歡
0817化學(xué)工程319求調(diào)劑
已經(jīng)有7人回復(fù)
-
083000環(huán)境科學(xué)與工程調(diào)劑
已經(jīng)有5人回復(fù)
-
求調(diào)劑 一志愿蘇州大學(xué)，0856化工323分 | 本科應(yīng)化 | 有專利/競賽/科研助手經(jīng)歷 |
已經(jīng)有8人回復(fù)
-
0856材料與化工353分求調(diào)劑
已經(jīng)有5人回復(fù)
-
282分材料專業(yè)求調(diào)劑院校
已經(jīng)有15人回復(fù)
-
化工284求調(diào)劑
已經(jīng)有12人回復(fù)
-
288求調(diào)劑
已經(jīng)有13人回復(fù)
-
一志愿安徽大學(xué)材料工程專碩313分，求調(diào)劑的學(xué)校
已經(jīng)有3人回復(fù)
-
考研一志愿長安大學(xué)材料與化工309分請求調(diào)劑
已經(jīng)有8人回復(fù)
-
085600材料與化工 326 求調(diào)劑
已經(jīng)有4人回復(fù)
1107040177
新蟲 (小有名氣)
- 應(yīng)助: 0 (幼兒園)
- 金幣: 386.4
- 散金: 10
- 帖子: 110
- 在線: 7.5小時
- 蟲號: 1888288
- 注冊: 2012-07-11
- 專業(yè): 食品科學(xué)基礎(chǔ)
2樓2019-05-29 17:37:25
【答案】應(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 這個中的TT時二維數(shù)組嗎？1：modifying_modifying(k,1)是想說從一到modifying_modifying(k,1)這個數(shù)，1^q,2^q,3^q......modifying_modifying(k,1)^q嗎？ |
» 本帖已獲得的紅花（最新10朵）
3樓2019-05-29 20:59:31
夢魘拾柒
木蟲 (著名寫手)
- 應(yīng)助: 13 (小學(xué)生)
- 金幣: 2225.7
- 散金: 6355
- 紅花: 4
- 帖子: 1027
- 在線: 240.3小時
- 蟲號: 987900
- 注冊: 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)小錯誤。應(yīng)為（k,1）,是一維數(shù)組。。兄臺能幫我調(diào)試一下，A為一維列向量即可。。 發(fā)自小木蟲Android客戶端 |
4樓2019-05-29 23:32:26
相關(guān)版塊跳轉(zhuǎn)
我要訂閱樓主 夢魘拾柒 的主題更新
| 4 | 1/1 | 返回列表 |
| 最具人氣熱帖推薦 [查看全部] | 作者 | 回/看 | 最后發(fā)表 | |
|---|---|---|---|---|
|
[考研] 求調(diào)劑 一志愿蘇州大學(xué)，0856化工323分 | 本科應(yīng)化 | 有專利/競賽/科研助手經(jīng)歷 | +6 | 橙子cyx 2026-03-06 | 8/400 |
|
|---|---|---|---|---|
|
|
[考研] 085600材料與化工 326 求調(diào)劑 +4 | 熱愛生活ing 2026-03-09 | 4/200 |
|
|
|
[考研] 考研材料與化工，求調(diào)劑 +3 | 戲精丹丹丹 2026-03-09 | 3/150 |
|
|
|
[考研] A區(qū)一本交叉課題組，低分調(diào)劑，招收機(jī)械電子信息通信等交叉方向 +54 | lisimayy 2026-03-04 | 67/3350 |
|
|
[考研] 334求調(diào)劑 +5 | Trying] 2026-03-07 | 5/250 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑 +4 | 很老實(shí)人 2026-03-09 | 4/200 |
|
|
|
[考研] 中科大材料299求調(diào)劑 +10 | DAIjiayo 2026-03-05 | 16/800 |
|
|
|
[考研] 求0856調(diào)劑 +5 | 小力氣珂珂 2026-03-08 | 5/250 |
|
|
|
[考研] 0703化學(xué)求調(diào)劑，總分309分，一志愿華南師范 +3 | 花與葉@ 2026-03-08 | 3/150 |
|
|
|
[考研] 理學(xué)287分求調(diào)劑 +4 | 看看我. 2026-03-05 | 4/200 |
|
|
|
[考研] 083000環(huán)境科學(xué)與工程調(diào)劑 +5 | 加油呀fxy 2026-03-07 | 6/300 |
|
|
|
[考研]
 調(diào)劑的同學(xué)，走過路過，不要錯過
+10
|
likeihood 2026-03-06 | 15/750 |
|
|
|
[考研] 334求調(diào)劑 +8 | Trying] 2026-03-06 | 8/400 |
|
|
|
[考博] 2026年博士名額撿漏 +4 | 科研ya 2026-03-04 | 7/350 |
|
|
|
[考研] 材料與化工304求調(diào)劑 +7 | 邱gl 2026-03-05 | 10/500 |
|
|
|
[考研] 求調(diào)劑，學(xué)校研究所都可以，材料與化工267分 +6 | wmx1 2026-03-05 | 6/300 |
|
|
|
[考研] 293求調(diào)劑 +3 | 是樂渝哇 2026-03-04 | 3/150 |
|
|
|
[考研] 復(fù)試調(diào)劑 +4 | 杰尼龜aaa 2026-03-03 | 4/200 |
|
|
|
[考研] 293求調(diào)劑 +4 | 是樂渝哇 2026-03-03 | 4/200 |
|
|
|
[考研] 306求調(diào)劑 +7 | 張張張張oo 2026-03-03 | 7/350 |
|