24小时热门版块排行榜

返回列表

当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖

百里道

木虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1797
散金: 5
帖子: 46
在线: 11.1小时
虫号: 1070800
注册: 2010-08-07
性别: GG
专业: 工程热物理相关交叉领域

[求助] matlab中二维图作图

已经作了温度关于时间和坐标的三维图，怎么做某一具体时刻的温度-位移图？
没有搜到这个，只好求救了。。

回复此楼

» 猜你喜欢

博士读完未来一定会好吗已经有15人回复
心脉受损已经有4人回复
Springer期刊投稿求助已经有4人回复
读博已经有3人回复
小论文投稿已经有3人回复
Bioresource Technology期刊，第一次返修的时候被退回好几次了已经有9人回复
到新单位后，换了新的研究方向，没有团队，持续积累2区以上论文，能申请到面上吗已经有8人回复
申请2026年博士已经有6人回复
请问哪里可以有青B申请的本子可以借鉴一下。已经有5人回复

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

关于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人回复

1楼2012-03-30 08:03:21

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

百里道

木虫 (初入文坛)

应助: 0 (幼儿园)
金币: 1797
散金: 5
帖子: 46
在线: 11.1小时
虫号: 1070800
注册: 2010-08-07
性别: GG
专业: 工程热物理相关交叉领域

引用回帖:

4楼: Originally posted by jinniu82 at 2012-03-30 11:24:03:
图片只是输出的一种方式，把代码中的输出部分的相关数据画出来不就行了，不行的话把代码放上来，研究研究

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