MATLAB中出现 索引超出矩阵维度,老是提示索引超出矩阵维度,求大神指导

function main

clear

L=32;%基底尺寸的设定;

step1=500;%生长层数的设定,

niu=1;deltt=0.05;deltx=1;dd=16;aa=(2*dd/deltx)^(1/2);%基本参数的设定;

cs=100;            %重复计算次数

h=zeros(L+2,step1);%基底高度的初值;

for k=1:step1-1

ww(k)=0;

end

for cc=1:cs;%重复计算

for t=2:step1

for i=3:L+1

%以下是分数生长方程的离散形式

h(1,t)=h(L+1,t);%边界条件赋值

h(L+2,t)=h(2,t);%边界条件赋值

random=unifrnd(-1/2,1/2,1,1);%随机噪声;随机数;

h(i,t+1)=h(i,t)-niu*deltt*(h(i+2,t)-4*h(i+1,t)+5*h(i,t)-4*h(i-1,t)+h(i-2,t))+aa*(12*deltt)^(1/2)*random;

end

end

% h(1:L+1,1:step1);

for k=1:step1-1;

w(k)=std(h(2:L+1,k+2));%surface width;surface roughnenss;表面粗糙度

ww(k)=ww(k)+w(k);

t(k)=(k+2)*deltt; %生长时间=生长层数step1*基底尺寸L*时间间隔deltt

end

ccc=cc

end

save date

plot(log2(t),log2(ww/cs),'square')

p=polyfit(log2(t),log2(ww/cs),1)

% plot(h(2:L,100),'.-')

  • 0
    点赞
  • 0
    收藏
    觉得还不错? 一键收藏
  • 0
    评论
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值