| 查看: 732 | 回复: 1 | ||
[求助]
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程序那里出错了 |
回复此楼» 猜你喜欢
十年后又回来了，论文投稿求助
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O.54,科目齐全,可+急
已经有4人回复
-
fileCode有新解读？
已经有17人回复
-
这样的filecode谁见过
已经有8人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有6人回复
-
售SCI一区T0P文章，我:8O.55.1.O.54,科目全,可伽急
已经有6人回复
-
售SCI文章，我:8O5.5.1.O.54,科目齐全,可+急
已经有6人回复
-
售SCI一区T0P文章，我:8.O.55.1.O54,科目全,可伽急
已经有5人回复
-
售SCI一区文章，我:8.O.55.1.O.54,科目齐全,可伽急
已经有4人回复
天天进步啊
木虫 (著名写手)
- 应助: 87 (初中生)
- 金币: 6058.7
- 散金: 101
- 红花: 42
- 沙发: 9
- 帖子: 1793
- 在线: 1005.5小时
- 虫号: 3689648
- 注册: 2015-02-15
- 性别: GG
- 专业: 理论和计算化学
2楼2018-03-28 14:37:22
20