| 查看: 1007 | 回复: 12 | |||
| 当前主题已经存档。 | |||
| 当前只显示满足指定条件的回帖，点击这里查看本话题的所有回帖 | |||
[交流]
【求助】matlab对符号方程组如何画出可视化图形【已解决】
|
|||
|
我想画出缺陷浓度随样品中某一元素含量变化的曲线，原来都是用matlab做一些简单的 计算，现在碰到符号方程组，不知道怎么处理啦，盼前辈能指点帮忙： 问题是：如何画出x1,x2,x3,x4随Xa变化的图像？（ x1,x2,x3,x4为样品中四种缺陷的浓度，Xa为样品中一种元素的含量，Xa取值区间为【0.70：0.75】） 下面的程序没有完成，现在不知道怎么下手啦：） k=1.380658e-23; %k 为玻尔兹曼常数 ee=1.60217733e-19; %ee为单位电荷电量 K=1000*k/ee; %为1000°时对应的能量 E1=1.47;E2=1.92;E3=0.56;E4=0.59; %E1,E2,E3,E4为空位形成能，单位eV A=4.5; % x1,x2,x3,x4为四种缺陷浓度 %Xa为样品中一种元素的含量，Xa取值区间为0.5：0.75 tic; syms x1 x2 x3 x4 Xa x1=0.75*A*(4*x4)/(1-x4)^(Xa-1)*exp((-E1+(Xa-1)*E4)/K); %x1缺陷浓度方程 x2=0.25*A*(4*x4/(1-x4))^Xa*exp((-E2+Xa*E4)/K); %x2缺陷浓度方程 x3=0.75*(4*x4/(1-x4))^(-1)*exp((-E3-E4)/K); %x3缺陷浓度方程 x4=x3+(1-Xa)*x1-Xa*x2+Xa-0.75; %x4缺陷浓度方程 solve(x1,x2,x3,x4) toc; [ Last edited by wangen994 on 2009-8-15 at 22:50 ] |
回复此楼» 猜你喜欢
一志愿北京理工大学本科211材料工程294求调剂
已经有8人回复
-
349求调剂
已经有6人回复
-
329求调剂
已经有9人回复
-
求调剂一志愿武汉理工大学材料工程（085601）
已经有9人回复
-
327求调剂
已经有5人回复
-
求收留
已经有5人回复
-
375求调剂
已经有4人回复
-
340求调剂
已经有6人回复
-
298求调剂
已经有3人回复
-
085601材料工程找调剂
已经有9人回复
hitzhang
木虫 (正式写手)
- 应助: 0 (幼儿园)
- 贵宾: 2.15
- 金币: 1376.7
- 散金: 969
- 红花: 8
- 帖子: 862
- 在线: 226.4小时
- 虫号: 390575
- 注册: 2007-06-02
- 性别: GG
- 专业: 无机非金属类电介质与电解
★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★ ★
positron(金币+20,VIP+0):谢谢，大概x1,x3,x4的形状是对的，x2的形状是相反的，呵呵，不过很不错啦 7-27 13:06
woshilsh(金币+8,VIP+0):辛苦了，张兄！ 7-27 17:02
positron(金币+5,VIP+0): 7-30 00:03
positron(金币+20,VIP+0):谢谢，大概x1,x3,x4的形状是对的，x2的形状是相反的，呵呵，不过很不错啦 7-27 13:06
woshilsh(金币+8,VIP+0):辛苦了，张兄！ 7-27 17:02
positron(金币+5,VIP+0): 7-30 00:03
|
function code: \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ function y=myfunction(x,Xa) k=1.380658e-23; %k 为玻尔兹曼常数 ee=1.60217733e-19; %ee为单位电荷电量 K=1000*k/ee; %为1000°时对应的能量 E1=1.47;E2=1.92;E3=0.56;E4=0.59; %E1,E2,E3,E4为空位形成能，单位eV A=4.5; x1=x(1); x2=x(2); x3=x(3); x4=x(4); y=[x1-0.75*A*(4*x4)/(1-x4)^(Xa-1)*exp((-E1+(Xa-1)*E4)/K);x2-0.25*A*(4*x4/(1-x4))^Xa*exp((-E2+Xa*E4)/K);x3-0.75*(4*x4/(1-x4))^(-1)*exp((-E3-E4)/K);x4-x3+(1-Xa)*x1-Xa*x2+Xa-0.75]; \\\\\\\\\\\\\\\\\\\\\\\\\\\\\ command code \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ >> Xa=0.7:0.001:0.75; >> Xaa=0.7:0.001:0.75; for i=1:51 Xa=Xaa(i); x0=1e-4*ones(1,4); f=@(x)myfunction(x,Xa); x(i,  = fsolve(f,x0);end >> plot(Xaa',x) >> xlabel('Xa') >> ylabel('x') >> legend({'x1','x2','x3','x4'}) \\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\\ |
8楼2009-07-27 12:55:56
2楼2009-07-27 00:08:39
woshilsh
荣誉版主 (职业作家)
- 应助: 43 (小学生)
- 贵宾: 10.33
- 金币: 21115.1
- 散金: 737
- 红花: 22
- 帖子: 3777
- 在线: 639.7小时
- 虫号: 293038
- 注册: 2006-11-04
- 专业: 计算机应用技术
- 管辖: 仿真模拟
3楼2009-07-27 02:06:46
5楼2009-07-27 09:42:10
20