声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 1763|回复: 1

[综合讨论] 怎么求出数值解并画出图形

[复制链接]
发表于 2010-7-13 18:43 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
程序为:
syms w1 m
xj=9; q=1.6*10^(-19); delta=10^(-4);
d1=1.5; d2=0.19; d3=0.4; d4=0.4;
N1=2.5*10^15; N2=2*10^16; N3=6*10^16; N4=2*10^15;
v=0.12; x=9:0.01:11.45;
e01=(q*N4*((xj+w1)^2-xj^2)/(2*delta))/xj
g1=int(e01*xj/m-q*N4*(m^2-xj^2)/(m*2*delta),m,xj,xj+w1)+v
w1=solve(g1,'w1');
e01=(q*N4*((xj+w1)^2-xj^2)/(2*delta))/xj;
y=(e01*xj./x-q*N4.*(x.^2-xj^2)./(x.*2*delta)).*(xj<=x&x<=xj+w1);
plot(x,y)
错误提示:

Warning: Explicit solution could not be found.
> In solve at 81
  In f1 at 17
??? Error using ==> mupadmex
Error in MuPAD command: Either base or exponent must be a scalar.
Error in ==> sym.sym>sym.mpower at 198
            B = mupadmex('mllib::mpower',A.s,p.s);
Error in ==> f1 at 18
e01=(q*N4*((xj+w1)^2-xj^2)/(2*delta))/xj;

[ 本帖最后由 ChaChing 于 2010-7-13 23:51 编辑 ]
回复
分享到:

使用道具 举报

发表于 2010-7-14 00:32 | 显示全部楼层
个人水平专业有限, LZ没交代下相关资讯, 真的有点难懂!
但好奇, 既然要求数值解, 为何还使用sym?
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-11 11:05 , Processed in 0.071187 second(s), 17 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表