24小时热门版块排行榜

返回列表

我是真滴皮

新虫 (初入文坛)

应助: 0 (幼儿园)
金币: 10
帖子: 1
在线: 1.6小时
虫号: 6924526
注册: 2017-07-09
专业: 物理学I

[求助] MATLAB萌新求助大神帮忙

X=[360,380,400,420,440,460,480,500]       %X轴,温度
Y=[-5,-4,-3,-2,-1,0,1,2,3]       %Y轴,log应变速率
Z=[55.62,68.62,91.44,105.20,44.37,51.56,74.29,87.50,32.57,43.24,60.21,74.37,24.74,31.58,43.30,58.74]             %X列，Y行
Zlog=log10(Z)       %log应力
i=50       %设置插值矩阵的行数
j=50       %设置插值矩阵的列数
Xi=(360: (500-360)/(j-1):500)       %赋值给X插值点
Ylogi=(-5: (3-(-5))/(i-1):3)       %赋值给Y插值点
Yi=10.^Ylogi       %应变速率的插值向量，10的Ylogi元素次方
Zlogi=interp2(X,Y,Zlog,Xi,Ylogi','spline')       %用spline方法插值(log应力)矩阵
Zi=interp2(X,Y,Z,Xi,Ylogi','spline')       %用spline方法插值应力矩阵
[Fx,M]=gradient(Zlogi,1,(3-(-5))/(i-1))       %计算log力/log速率=M。如果F是二维矩阵，返回F的二维数值梯度。[FX,FY]=gradient(F,HX,HY)。HX,HY参数表示各方向相邻两点的距离。
S=cumtrapz(Yi,Zi)       %通过递归循环计算积分,计算累积梯形积分,应力在log应变速率上的点积分
P=zeros(i,j)       %设置全0矩阵
Es=-5       %赋值log应变速率最小值
Smin=Zi(1,: ).*10^Es./(M(1,: )+1)       %计算应变速率最小时的第一项定值,(3,: )代表第三行，所有的列。(:,3)代表所有行,第三列。(:,1,1)代表所有第一维,第二维和第三维的第一个元素。
for n=1:j
G(:,n)=S(:,n)+Smin(n)       %每项的累积积分值加上最小值-G值
P(:,n)=Zi(:,n).*Yi(: )       %每项应力和应变速率的乘积-P值
end
eta=2*(1-G./P)*100       %计算η值矩阵
ksai=2*M./(eta/100)-1       %计算ζ值矩阵
[C,h] =contour(Xi,Ylogi,eta,8)          %η二维图，contour(Xi,Ylogi,eta,10) ,10表示10根等高线，不加10表示自动
clabel(C,h)
求大佬看一下这个matlab程序那里出错了

回复此楼