早教吧作业答案频道 -->其他-->
求助MATLAB符号计算:%声明符号变量syms c w t b a D n T v Q x y;%给已知的变量赋值D=0.038;w=2*pi/T;v=n*w*D/2;%正切值b= (sin (w*t))/(c-cos (w*t));%余弦值a=sqrt (1/ (1+b^2));%正弦值d=sqrt ((b^2)/ (1+b^2));x=int (v*sqrt (1/ (1
题目详情
求助MATLAB符号计算:
%声明符号变量
syms c w t b a D n T v Q x y;
%给已知的变量赋值
D=0.038;
w=2*pi/T;
v=n*w*D/2;
%正切值
b= (sin (w*t))/(c-cos (w*t));
%余弦值
a=sqrt (1/ (1+b^2));
%正弦值
d=sqrt ((b^2)/ (1+b^2));
x=int (v*sqrt (1/ (1+b^2)),t);
y=int (v*sqrt ((b^2)/ (1+b^2)),t);
x=simple(x)
y=simple(y)
Warning:Explicit integral could not be found.
x =
int((19*pi*n*(1/(sin((2*pi*t)/T)^2/(c - cos((2*pi*t)/T))^2 + 1))^(1/2))/(500*T),t)
y =
(19*n*(sin((2*pi*t)/T)^2/(2*c*(2*sin((pi*t)/T)^2 - 1) + c^2 + 1))^(1/2)*(2*c*(2*sin((pi*t)/T)^2 - 1) + c^2 + 1))/(1000*c*sin((2*pi*t)/T))
请问为什么y化得出,x化不出?
可以给出x的结果吗?
%声明符号变量
syms c w t b a D n T v Q x y;
%给已知的变量赋值
D=0.038;
w=2*pi/T;
v=n*w*D/2;
%正切值
b= (sin (w*t))/(c-cos (w*t));
%余弦值
a=sqrt (1/ (1+b^2));
%正弦值
d=sqrt ((b^2)/ (1+b^2));
x=int (v*sqrt (1/ (1+b^2)),t);
y=int (v*sqrt ((b^2)/ (1+b^2)),t);
x=simple(x)
y=simple(y)
Warning:Explicit integral could not be found.
x =
int((19*pi*n*(1/(sin((2*pi*t)/T)^2/(c - cos((2*pi*t)/T))^2 + 1))^(1/2))/(500*T),t)
y =
(19*n*(sin((2*pi*t)/T)^2/(2*c*(2*sin((pi*t)/T)^2 - 1) + c^2 + 1))^(1/2)*(2*c*(2*sin((pi*t)/T)^2 - 1) + c^2 + 1))/(1000*c*sin((2*pi*t)/T))
请问为什么y化得出,x化不出?
可以给出x的结果吗?
▼优质解答
答案和解析
奇怪,怎么我的MATLAB和你运行结果不一样?
>> %声明符号变量
syms c w t b a D n T v Q x y;
%给已知的变量赋值
D=0.038;
w=2*pi/T;
v=n*w*D/2;
%正切值
b= (sin (w*t))/(c-cos (w*t));
%余弦值
a=sqrt (1/ (1+b^2));
%正弦值
d=sqrt ((b^2)/ (1+b^2));
x=int (v*sqrt (1/ (1+b^2)),t);
y=int (v*sqrt ((b^2)/ (1+b^2)),t);
x=simple(x)
y=simple(y)
x =
-19/1000*n*(c^2*cos(2*pi/T*t)-c^2+2*c*cos(2*pi/T*t)-i*(c^2-2*c*cos(2*pi/T*t)+1)^(1/2)*c*EllipticE(i*sin(2*pi/T*t)/(1+cos(2*pi/T*t)),((c^2+2*c+1)/(c^2-2*c+1))^(1/2))*sin(2*pi/T*t)+2*i*(c^2-2*c*cos(2*pi/T*t)+1)^(1/2)*EllipticF(i*sin(2*pi/T*t)/(1+cos(2*pi/T*t)),((c^2+2*c+1)/(c^2-2*c+1))^(1/2))*c*sin(2*pi/T*t)-2*c*cos(2*pi/T*t)^2+cos(2*pi/T*t)-1+i*(c^2-2*c*cos(2*pi/T*t)+1)^(1/2)*EllipticE(i*sin(2*pi/T*t)/(1+cos(2*pi/T*t)),((c^2+2*c+1)/(c^2-2*c+1))^(1/2))*sin(2*pi/T*t))/sin(2*pi/T*t)/c/(c^2-2*c*cos(2*pi/T*t)+1)^(1/2)
y =
19/1000*i*n*(c^2-2*c*cos(2*pi/T*t)+1)^(1/2)*(-1+cos(2*pi/T*t)^2)^(1/2)/sin(2*pi/T*t)/c
>>
%%%%%%%%%%%%%%%%%%%%%%%%%
以上完全复制MATLAB界面.
>> %声明符号变量
syms c w t b a D n T v Q x y;
%给已知的变量赋值
D=0.038;
w=2*pi/T;
v=n*w*D/2;
%正切值
b= (sin (w*t))/(c-cos (w*t));
%余弦值
a=sqrt (1/ (1+b^2));
%正弦值
d=sqrt ((b^2)/ (1+b^2));
x=int (v*sqrt (1/ (1+b^2)),t);
y=int (v*sqrt ((b^2)/ (1+b^2)),t);
x=simple(x)
y=simple(y)
x =
-19/1000*n*(c^2*cos(2*pi/T*t)-c^2+2*c*cos(2*pi/T*t)-i*(c^2-2*c*cos(2*pi/T*t)+1)^(1/2)*c*EllipticE(i*sin(2*pi/T*t)/(1+cos(2*pi/T*t)),((c^2+2*c+1)/(c^2-2*c+1))^(1/2))*sin(2*pi/T*t)+2*i*(c^2-2*c*cos(2*pi/T*t)+1)^(1/2)*EllipticF(i*sin(2*pi/T*t)/(1+cos(2*pi/T*t)),((c^2+2*c+1)/(c^2-2*c+1))^(1/2))*c*sin(2*pi/T*t)-2*c*cos(2*pi/T*t)^2+cos(2*pi/T*t)-1+i*(c^2-2*c*cos(2*pi/T*t)+1)^(1/2)*EllipticE(i*sin(2*pi/T*t)/(1+cos(2*pi/T*t)),((c^2+2*c+1)/(c^2-2*c+1))^(1/2))*sin(2*pi/T*t))/sin(2*pi/T*t)/c/(c^2-2*c*cos(2*pi/T*t)+1)^(1/2)
y =
19/1000*i*n*(c^2-2*c*cos(2*pi/T*t)+1)^(1/2)*(-1+cos(2*pi/T*t)^2)^(1/2)/sin(2*pi/T*t)/c
>>
%%%%%%%%%%%%%%%%%%%%%%%%%
以上完全复制MATLAB界面.
看了 求助MATLAB符号计算:%...的网友还看了以下: