jljhld 发表于 2008-12-11 14:17

关于用matlab求完积分画图的问题

我的程序如下:
syms x;
y=-1.5:0.01:1.5;
p1=int(x*exp((-0.25*x^4+x.*y)/1)./(int(exp((-0.25*x^4+y*x)/1),x,-inf,inf)),x,-inf,inf);
plot(y,double(p1))
但是画出图形后出现(Warning: Imaginary parts of complex X and/or Y arguments ignored.)我不知道什么问题,请高人指点。

sogooda 发表于 2008-12-11 14:23

p1中有虚部,绘图时已忽略。

jljhld 发表于 2008-12-12 11:04

p1中有虚部,绘图时已忽略。

我也知道啊,能不能改进一下。不出现(Warning: Imaginary parts of complex X and/or Y arguments ignored.)

sigma665 发表于 2008-12-12 16:05

回复 板凳 jljhld 的帖子

你可以忽略,但是得到的结果是复数域的,画图没有意义

ChaChing 发表于 2008-12-12 21:21

回复 楼主 jljhld 的帖子

symbolic这块个人接触少, 但几点小地方, 还是请楼主参考下! 个人认为写程序还是要养成好习惯
1. x.*y中间的点乘多馀!
2. (-0.25*x^4+x.*y)/1 除法多馀, 直接写成 -0.25*x^4+x.*y
3. -0.25*x^4+x*y 与 -0.25*x^4+y*x 意义相同, 干嘛写不同
所以 p1=int(x*exp((-0.25*x^4+x.*y)/1)./(int(exp((-0.25*x^4+y*x)/1),x,-inf,inf)),x,-inf,inf);
修改为 ee=exp(-0.25*x^4+x*y); p1=int( x*ee./int(ee,x,-inf,inf),x,-inf,inf);
参考参考!
页: [1]
查看完整版本: 关于用matlab求完积分画图的问题