关于matlab积分的使用,syms w x m oj b r q ob o h s0 b0 t p M kf1=x*{exp(w*m*i-0.5*w^2*oj^2)-1-w*m*i};f2=b - r -q +0.5*o^2 - h*o^2 + subs(f1,'w',-i*h-i)-subs(f1,'w',-i*h);f3=x * exp( m * h + 0.5 * h^2 * oj^2) * ( exp( i* w * (m + h * oj^2 ) - 0.
来源:学生作业帮助网 编辑:作业帮 时间:2024/06/30 07:34:02
![关于matlab积分的使用,syms w x m oj b r q ob o h s0 b0 t p M kf1=x*{exp(w*m*i-0.5*w^2*oj^2)-1-w*m*i};f2=b - r -q +0.5*o^2 - h*o^2 + subs(f1,'w',-i*h-i)-subs(f1,'w',-i*h);f3=x * exp( m * h + 0.5 * h^2 * oj^2) * ( exp( i* w * (m + h * oj^2 ) - 0.](/uploads/image/z/14960965-13-5.jpg?t=%E5%85%B3%E4%BA%8Ematlab%E7%A7%AF%E5%88%86%E7%9A%84%E4%BD%BF%E7%94%A8%2Csyms+w+x+m+oj+b+r+q+ob+o+h+s0+b0+t+p+M+kf1%3Dx%2A%7Bexp%28w%2Am%2Ai-0.5%2Aw%5E2%2Aoj%5E2%29-1-w%2Am%2Ai%7D%3Bf2%3Db+-+r+-q+%2B0.5%2Ao%5E2+-+h%2Ao%5E2+%2B+subs%28f1%2C%27w%27%2C-i%2Ah-i%29-subs%28f1%2C%27w%27%2C-i%2Ah%29%3Bf3%3Dx+%2A+exp%28+m+%2A+h+%2B+0.5+%2A+h%5E2+%2A+oj%5E2%29+%2A+%28+exp%28+i%2A+w+%2A+%28m+%2B+h+%2A+oj%5E2+%29+-+0.)
关于matlab积分的使用,syms w x m oj b r q ob o h s0 b0 t p M kf1=x*{exp(w*m*i-0.5*w^2*oj^2)-1-w*m*i};f2=b - r -q +0.5*o^2 - h*o^2 + subs(f1,'w',-i*h-i)-subs(f1,'w',-i*h);f3=x * exp( m * h + 0.5 * h^2 * oj^2) * ( exp( i* w * (m + h * oj^2 ) - 0.
关于matlab积分的使用,
syms w x m oj b r q ob o h s0 b0 t p M k
f1=x*{exp(w*m*i-0.5*w^2*oj^2)-1-w*m*i};
f2=b - r -q +0.5*o^2 - h*o^2 + subs(f1,'w',-i*h-i)-subs(f1,'w',-i*h);
f3=x * exp( m * h + 0.5 * h^2 * oj^2) * ( exp( i* w * (m + h * oj^2 ) - 0.5 * w^2 * oj^2) -1 - i * w *(m + h * oj^2) );
cf1=exp( i * w * log ( s0 / b0 ) + t * (0.5 * i * w * o^2) - 0.5 * w^2 * o^2 + f3 + subs(f3 ,'w' ,i) );
cf2=exp( i * w * log ( s0 / b0 ) + t * (-0.5 * i * w * o^2) - 0.5 * w^2 * o^2 + subs(f3,'h',h+1 ) - i * w * subs( subs(f3,'h',h+1 ) ,'w' ,-i) );
f4=exp(-i * w * k * p / M )/(i*w) * cf1
f5=int(f4,w,-10000,10000)
这里计算总是报错
Operation terminated by user during ==> mupadengine.mupadengine>mupadengine.evalin at 119
In ==> mupadengine.mupadengine>mupadengine.feval at 157
[S,err] = evalin(engine,stmt);
In ==> mupadfeval at 24
[S,err] = feval(symengine,cmd,varargin{:}); % returns sym
In ==> map at 14
[result,status] = mupadfeval(fcn,obj,varargin{:});
In ==> sym.int at 56
r = reshape(map(f(:),'mlint',x,a,b),size(f));
In ==> test at 8
f5=int(f4,w,-10000,10000)
关于matlab积分的使用,syms w x m oj b r q ob o h s0 b0 t p M kf1=x*{exp(w*m*i-0.5*w^2*oj^2)-1-w*m*i};f2=b - r -q +0.5*o^2 - h*o^2 + subs(f1,'w',-i*h-i)-subs(f1,'w',-i*h);f3=x * exp( m * h + 0.5 * h^2 * oj^2) * ( exp( i* w * (m + h * oj^2 ) - 0.
我算了一下,没有解析解
syms w x m oj b r q ob o h s0 b0 t p M k
f1=x*{exp(w*m*i-0.5*w^2*oj^2)-1-w*m*i};
f2=b - r -q +0.5*o^2 - h*o^2 + subs(f1,'w',-i*h-i)-subs(f1,'w',-i*h);
f3=x * exp( m * h + 0.5 * h^2 * oj^2) * ( exp( i* w * (m + h * oj^2 ) - 0.5 * w^2 * oj^2) -1 - i * w *(m + h * oj^2) );
cf1=exp( i * w * log ( s0 / b0 ) + t * (0.5 * i * w * o^2) - 0.5 * w^2 * o^2 + f3 + subs(f3 , 'w' ,i) );
cf2=exp( i * w * log ( s0 / b0 ) + t * (-0.5 * i * w * o^2) - 0.5 * w^2 * o^2 + subs(f3,'h',h+1 ) - i * w * subs( subs(f3,'h',h+1 ) , 'w' ,-i) );
f4=exp(-i * w * k * p / M )/(i*w) * cf1
f5=int(f4,w,-10000,10000)
f4 =
-(i*exp(i*w*log(s0/b0) - (o^2*w^2)/2 + x*exp((h^2*oj^2)/2 + m*h)*(m + exp(oj^2/2 - h*oj^2 - m) + h*oj^2 - 1) - x*exp((h^2*oj^2)/2 + m*h)*(i*w*(h*oj^2 + m) - exp(i*w*(h*oj^2 + m) - 1/2*oj^2*w^2) + 1) + (i*o^2*t*w)/2))/(w*exp((i*k*p*w)/M))
Warning: Explicit integral could not be found.
> In sym.int at 64
f5 =
int(-(i*exp(i*w*log(s0/b0) - (o^2*w^2)/2 + x*exp((h^2*oj^2)/2 + m*h)*(m + exp(oj^2/2 - h*oj^2 - m) + h*oj^2 - 1) - x*exp((h^2*oj^2)/2 + m*h)*(i*w*(h*oj^2 + m) - exp(i*w*(h*oj^2 + m) - 1/2*oj^2*w^2) + 1) + (i*o^2*t*w)/2))/(w*exp((i*k*p*w)/M)), w = -10000..10000)
>>