24小时热门版块排行榜    

查看: 2095  |  回复: 4

vivienmama

铜虫 (小有名气)

[求助] matlab画热加工图遇到的小问题

最近在用matlab画热加工图,遇到如下问题,请大侠们来帮忙!
(1)功率耗散功率等高图上的值 怎么让它只显示小数点后两位?
(2)热加工失稳图怎样把小于0的设置为阴影?

附件:


clc;clear;
X=[900,1050,1200]        %X轴,温度
Y=[-3,-1,1]        %Y轴,log应变速率
Z=[90.49,39.2,12.44,
196.36,96.3,48,
241.92,154.71,85.18]                %X列,Y行
Zlog=log10(Z)        %log应力
i=50        %设置插值矩阵的行数
j=50        %设置插值矩阵的列数
Xi=(900: (1200-900)/(j-1):1200)        %赋值给X插值点
Ylogi=(-3: (1-(-3))/(i-1):1)        %赋值给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,(1-(-3))/(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=-3        %赋值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,20)           %η二维图,contour(Xi,Ylogi,eta,10) ,10表示10根等高线,不加10表示自动
clabel(C,h)

1.jpg
111.jpg



[ Last edited by vivienmama on 2013-1-22 at 22:20 ]
回复此楼

» 本帖附件资源列表

  • 欢迎监督和反馈:小木虫仅提供交流平台,不对该内容负责。
    本内容由用户自主发布,如果其内容涉及到知识产权问题,其责任在于用户本人,如对版权有异议,请联系邮箱:xiaomuchong@tal.com
  • 附件 1 : untitled1.fig
  • 2013-01-22 22:17:44, 41 K

» 猜你喜欢

» 本主题相关价值贴推荐,对您同样有帮助:

,,,
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

yuhuidavid

木虫 (著名写手)

【答案】应助回帖


Zpffpz: 金币+1, 专家考核, 好建议 2013-01-30 15:38:38
改用origin,很容易的说~!
狭路相逢勇者胜!
2楼2013-01-30 12:08:24
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaomuchong916

木虫 (著名写手)

画出等于零的等值线,再采用fill填充即可。把等值线的handle先通过str2num转化,再通过sprintf设置有效数字位数,再以string显示即可。

[ 发自手机版 http://muchong.com/3g ]

» 本帖已获得的红花(最新10朵)

3楼2013-01-30 21:39:13
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

xiaomuchong916

木虫 (著名写手)

【答案】应助回帖

画出等于零的等值线,再采用fill填充即可。把等值线的handle先通过str2num转化,再通过sprintf设置有效数字位数,再以string显示即可。

[ 发自手机版 http://muchong.com/3g ]
4楼2013-01-30 21:39:36
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖

vivienmama

铜虫 (小有名气)

送鲜花一朵
引用回帖:
3楼: Originally posted by xiaomuchong916 at 2013-01-30 21:39:13
画出等于零的等值线,再采用fill填充即可。把等值线的handle先通过str2num转化,再通过sprintf设置有效数字位数,再以string显示即可。

谢谢呀!
不过我matlab基础很差,上面的程序是参考别人论文的。
小女子时间有限,不能深入学习matlab
求大侠把这些改成具体的语句吧,相信这对你来说 很容易!
,,,
5楼2013-02-26 17:23:43
已阅   回复此楼   关注TA 给TA发消息 送TA红花 TA的回帖
相关版块跳转 我要订阅楼主 vivienmama 的主题更新
最具人气热帖推荐 [查看全部] 作者 回/看 最后发表
[基金申请] 面上函评意见出来了,像什么等级? 20+3 Tsingking1 2026-08-27 14/700 2026-09-01 08:29 by 尼古拉斯小虫
[基金申请] 怎么看青基中了没有啊 +6 叶九微 2026-08-26 6/300 2026-08-31 23:54 by yudaoqian88
[基金申请] 投票:  有多少人是今天查系统知道结果的? +17 爱看书的可乐 2026-08-26 19/950 2026-08-31 21:30 by xiangy672
[基金申请] 基金未中,这种答复是模板吗? +6 zhaosm1982 2026-08-27 7/350 2026-08-31 21:18 by qdxxmc
[基金申请] 哪位高人中了,把查询到的截图贴出来让我看看,让我长长见识 +6 yuleib84 2026-08-26 7/350 2026-08-31 19:46 by 鱼翔浅底1
[基金申请] 面上意见出来了 +12 黄鸟于飞Chao 2026-08-29 23/1150 2026-08-31 18:57 by 黄鸟于飞Chao
[基金申请] 能否申诉? +7 echo8914667 2026-08-30 8/400 2026-08-31 17:00 by yihongxu
[文学芳草园] 梦想 +5 myrtle 2026-08-26 5/250 2026-08-31 14:40 by hahaboy
[基金申请] 有没有仍没收到信息的 +7 德尚中行 2026-08-27 8/400 2026-08-30 20:52 by purplejack
[考博] 找导师 +6 yuanjiabao 2026-08-29 7/350 2026-08-30 14:40 by 生科新手
[基金申请] 我就是申请一个面上项目而已,这评审意见是按照杰青的条件评的吧? +6 gouxfjh 2026-08-28 11/550 2026-08-30 07:57 by gouxfjh
[基金申请] 国自然面上复盘~欢迎讨论 (金币+15) +15 晴天加油 2026-08-26 16/800 2026-08-29 18:28 by symmetry
[基金申请] 怎么查啊 +6 huang1991js 2026-08-26 6/300 2026-08-28 08:42 by winsaint
[基金申请] 面上合作单位盖章 +5 ssyjh 2026-08-27 5/250 2026-08-27 20:50 by gdfollow
[基金申请] 2026年8月25日国自然放榜前突然收到列入评审专家邮件,有关系吗? +25 木水思豆 2026-08-25 28/1400 2026-08-26 14:53 by draco1987
[基金申请] 为什么国自然不能直接公布 +4 bjdxyxy 2026-08-26 4/200 2026-08-26 13:12 by qingmu1201
[基金申请] 国际合作可查了,中了面上 (EPI+1)(金币+50) +18 Ldrop2023 2026-08-26 18/900 2026-08-26 11:15 by cmrandy
[基金申请] 牛来!米来!面来! +8 beefly 2026-08-26 8/400 2026-08-26 08:37 by xuzhipiao
[基金申请] 在坚冰还盖着北海的时候,我看到了怒放的梅花。 (金币+10) +6 ziyangfang 2026-08-25 9/450 2026-08-25 20:26 by huagongfeihu
[基金申请] 某些机构,以效率低为荣,以效率低作为存在感 +9 yuleib84 2026-08-25 10/500 2026-08-25 17:14 by alexon
信息提示
请填处理意见