你的参数没有定义,我定义为:a=1;h=2;w=pi;
syms t
f1=int('cos(atan(sin(2*t)/2))')
f2=int('sin(atan(sin(2*t)/2))')
对f1,f2求积分,
f1 =1/2*(-(4+sin(2*t)^2)*(-1+sin(2*t)^2))^(1/2)*(1-sin(2*t)^2)^(1/2)/(4-sin(2*t)^4-3*sin(2*t)^2)^(1/2)*EllipticF(sin(2*t),1/2*i)/cos(2*t)
f2 =
1/4*(-(4+sin(2*t)^2)*(-1+sin(2*t)^2))^(1/2)*asin(2/5*sin(2*t)^2+3/5)/cos(2*t)/(4+sin(2*t)^2)^(1/2)
将其定义为函数文件:
function y=test11(t)
y=1/2.*(-(4+sin(2.*t).^2).*(-1+sin(2.*t).^2)).^(1/2).*(1-sin(2.*t).^2).^(1/2)./(4-sin(2.*t).^4-3.*sin(2.*t).^2).^(1/2)./cos(2.*t);
function y=test22(t)
y=1/4.*(-(4+sin(2.*t).^2).*(-1+sin(2.*t).^2)).^(1/2).*asin(2/5*sin(2.*t).^2+3/5)./cos(2.*t)./(4+sin(2.*t).^2).^(1/2);
命令窗口调用:
t=-2*pi:0.001:2*pi;
>> plot(t,test11(t),t,test22(t))
注意:
f1中的EllipticF(sin(2*t),1/2*i) 这个函数我的matlab说没有定义,可能是少了某个工具包。。
我把它删掉后再画图的,
但是大概意思你应该懂了。
结果:这个图像可能就是删掉那函数变成这样的,不太完整