| 查看: 2811 | 回复: 8 | |||
百里道木虫 (初入文坛)
|
[求助]
matlab中二维图作图
|
|
已经作了温度关于时间和坐标的三维图，怎么做某一具体时刻的温度-位移图？ 没有搜到这个，只好求救了。。 |
回复此楼» 猜你喜欢
070300化学 求调剂
已经有12人回复
-
化工求调剂！
已经有32人回复
-
291 求调剂
已经有6人回复
-
297求调剂
已经有23人回复
-
283求调剂，工科！
已经有6人回复
-
266求调剂
已经有22人回复
-
287求调剂
已经有18人回复
-
300求调剂
已经有6人回复
-
282，求调剂
已经有3人回复
-
071000生物学调剂求助
已经有8人回复
» 本主题相关价值贴推荐，对您同样有帮助:
- 关于Matlab作图的问题
已经有14人回复
- Matlab画图求助? (有人说, 标题一定要长------------------------->)
已经有5人回复
- Matlab中如何调用M文件作图
已经有9人回复
- 求助matlab三维画图中：Inner matrix dimensions must agree.
已经有4人回复
- MATLAB 画图小问题
已经有11人回复
- 求助：谁会用matlab画图？
已经有13人回复
- Matlab 画图后怎么获得图上的数据点
已经有6人回复
- 【求助】matlab 画出类似的二维图
已经有3人回复
- 科研作图时候matlab中使用latex
已经有3人回复
- 【求助】MATLAB多个mat中 数据画图问题
已经有7人回复
- 【求助】matlab做二维彩图【已解决】
已经有7人回复
- 【求助】如何采用Matlab绘制电子云密度图
已经有8人回复
jinniu82
铁杆木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 7338.8
- 红花: 1
- 帖子: 352
- 在线: 49.9小时
- 虫号: 711197
- 注册: 2009-02-28
- 专业: 光谱技术
2楼2012-03-30 08:12:45
百里道
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1797
- 散金: 5
- 帖子: 46
- 在线: 11.1小时
- 虫号: 1070800
- 注册: 2010-08-07
- 性别: GG
- 专业: 工程热物理相关交叉领域
3楼2012-03-30 08:50:18
jinniu82
铁杆木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 7338.8
- 红花: 1
- 帖子: 352
- 在线: 49.9小时
- 虫号: 711197
- 注册: 2009-02-28
- 专业: 光谱技术
4楼2012-03-30 11:24:03
dbb627
荣誉版主 (著名写手)
-
专家经验: +4 - 仿真EPI: 6
- 应助: 289 (大学生)
- 贵宾: 0.589
- 金币: 24640.4
- 散金: 551
- 红花: 61
- 沙发: 1
- 帖子: 1246
- 在线: 1794.8小时
- 虫号: 149791
- 注册: 2005-12-29
- 性别: GG
- 专业: 污染控制化学
- 管辖: 计算模拟
5楼2012-03-30 11:27:34
百里道
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1797
- 散金: 5
- 帖子: 46
- 在线: 11.1小时
- 虫号: 1070800
- 注册: 2010-08-07
- 性别: GG
- 专业: 工程热物理相关交叉领域
|
function yiweifourier c=1; %Phi=10; Lambda=10; qw=-10; %d=Phi/Lambda; d=@(x)x^2;g=@(t)exp(t); e=qw/Lambda; xspan=[0 1]; tspan=[0 0.5]; ngrid=[100 10]; f=@(x)4*x; [T,x,t]=rechuandao(c,d,e,g,f,xspan,tspan,ngrid); [x,t]=meshgrid(x,t); mesh(x,t,T); xlabel('x') ylabel('t') zlabel('T') function [U,x,t]=rechuandao(c,d,e,g,f,xspan,tspan,ngrid) n=ngrid(1); m=ngrid(2); h=range(xspan)/(m-1); x=linspace(xspan(1),xspan(2),m); k=range(tspan)/(n-1); t=linspace(tspan(1),tspan(2),n); r=c^2*k/h^2; if r>0.5 error('<为了保证算法的收敛，请增大步长h或减小步长k!') end s=1-2*r; s2=e*h; U=zeros(ngrid); for i=1:m U(1,i)=f(x(i)); end for j=2:n for i=2:m-1 U1(j,i)=d(x(i));U2(j,i)=g(t(j)); U(j,i)=s*U(j-1,i)+r*(U(j-1,i-1)+U(j-1,i+1))+50*U1(j,i)*U2(j,i); end %%%%%%%下面是边界条件 U(j,1)=s2+U(j,2);%定热流边界条件 U(j,m)=U(j,m-1);%绝热边界条件 end |
6楼2012-04-04 21:04:20
jinniu82
铁杆木虫 (正式写手)
- 应助: 2 (幼儿园)
- 金币: 7338.8
- 红花: 1
- 帖子: 352
- 在线: 49.9小时
- 虫号: 711197
- 注册: 2009-02-28
- 专业: 光谱技术
7楼2012-04-05 09:25:13
百里道
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1797
- 散金: 5
- 帖子: 46
- 在线: 11.1小时
- 虫号: 1070800
- 注册: 2010-08-07
- 性别: GG
- 专业: 工程热物理相关交叉领域
8楼2012-04-05 15:23:56
百里道
木虫 (初入文坛)
- 应助: 0 (幼儿园)
- 金币: 1797
- 散金: 5
- 帖子: 46
- 在线: 11.1小时
- 虫号: 1070800
- 注册: 2010-08-07
- 性别: GG
- 专业: 工程热物理相关交叉领域
9楼2012-04-10 08:36:49
