smart55427
2008-11-08, 23:05
这是我在命令编辑窗口中写的一个程序
x3=sym('exp(i*2*pi*t/3)+exp(i*pi*t)');
x33=subs(x3,'-i','i');
%E3=sym(int(x3*conj(x3),-3,3))
E3=sym(int(x3*x33,'-a','a'))
subs(E3,'a',3)
运行结果里出现
E3 =
2*(2*i*pi*a-3*exp(-1/3*i*pi*a)+3*exp(1/3*i*pi*a))/i/pi
ans =
12
如果把倒数第二行注释掉，而把倒数第三行取消注释会出现
Warning: Explicit integral could not be found.
> In D:\MATLAB6p5p1\toolbox\symbolic\@sym\int.m at line 58
In D:\MATLAB6p5p1\work\E1_E2_E3.m at line 3
E3 =
int((exp(2/3*i*pi*t)+exp(i*pi*t))*conj(exp(2/3*i*pi*t)+exp(i*pi*t)),t = -3 .. 3)
ans =
int((exp(2/3*i*pi*t)+exp(i*pi*t))*conj(exp(2/3*i*pi*t)+exp(i*pi*t)),t = -3 .. 3)
这是什么原因，而且倒数第二行和第三行有什么区别吗？希望高手指点指点
x3=sym('exp(i*2*pi*t/3)+exp(i*pi*t)');
x33=subs(x3,'-i','i');
%E3=sym(int(x3*conj(x3),-3,3))
E3=sym(int(x3*x33,'-a','a'))
subs(E3,'a',3)
运行结果里出现
E3 =
2*(2*i*pi*a-3*exp(-1/3*i*pi*a)+3*exp(1/3*i*pi*a))/i/pi
ans =
12
如果把倒数第二行注释掉，而把倒数第三行取消注释会出现
Warning: Explicit integral could not be found.
> In D:\MATLAB6p5p1\toolbox\symbolic\@sym\int.m at line 58
In D:\MATLAB6p5p1\work\E1_E2_E3.m at line 3
E3 =
int((exp(2/3*i*pi*t)+exp(i*pi*t))*conj(exp(2/3*i*pi*t)+exp(i*pi*t)),t = -3 .. 3)
ans =
int((exp(2/3*i*pi*t)+exp(i*pi*t))*conj(exp(2/3*i*pi*t)+exp(i*pi*t)),t = -3 .. 3)
这是什么原因，而且倒数第二行和第三行有什么区别吗？希望高手指点指点