task:19&20 微分方程(8.6-8.9)

语法知识

补充：（8.5 极限的计算）

8.6 常微分方程的符号解

8.7 平面曲线族的包络线

实例演练

8.6

♪解二阶微分方程：

>> dsolve('D2y+4*y=sin(x)','x')
 
ans =
 
sin(5*x)/24 - sin(3*x)/8 + sin(x)/6 - sin(2*x)*(cos(2*x)/6 - cos(x)/3 + cos(x)*(cos(2*x)/6 - cos(x)/3)) + C1*cos(2*x) + C2*sin(2*x)
 
>> simplify(ans)
 
ans =
 
sin(2*x)/6 + sin(x)/3 + C1*cos(2*x) + C2*sin(2*x)
 

♪求回路电流：

>> dsolve('DI+I/T=I0*sin(w*t)','I(0)=0,DI(0)=0')
 
ans =
 
(I0*w*exp(-t/T))/(1/T^2 + w^2) - (I0*(w*cos(t*w) - sin(t*w)/T))/(1/T^2 + w^2)
 
>> I=simplify(ans)
 
I =
 
(I0*w*exp(-t/T))/(1/T^2 + w^2) - (I0*(w*cos(t*w) - sin(t*w)/T))/(1/T^2 + w^2)
 

♪追线问题：

此题感悟…QAQ(一定要这么为难人嘛)：首先哈，这个图先是看不懂，本以为导弹曲线应该是一个开口朝下的抛物线。 其次前期工作这个流程大概看半懂，具体解题方法还是懵逼的，接下来的代码应该是前期工作的从后往前求解从而得到我们所需要的。（♞：之后有点脑子了可以回来康康这题）

8.7

♪直线滑动的包络线：

>> clear
>> syms x y a c theta
>> f=x/a+y/(c^2-a^2)^(1/2)-1;
>> eq=diff(f,'a')
 
eq =
 
(a*y)/(c^2 - a^2)^(3/2) - x/a^2
 
>> [x,y]=solve(f,eq)
 
x =
 
a^3/((c^2 - a^2)^(3/2)*(1/(c^2 - a^2)^(1/2) + a^2/(c^2 - a^2)^(3/2)))
 
 
y =
 
1/(1/(c^2 - a^2)^(1/2) + a^2/(c^2 - a^2)^(3/2))
 
>> simplify(x)
 
ans =
 
a^3/c^2
 
>> simplify(y)
 
ans =
 
(c^2 - a^2)^(3/2)/c^2
 
>> x=c*cos(theta)^3;
>> y=c*sin(theta)^3;
>> x^(2/3)+y^(2/3)
 
ans =
 
(c*cos(theta)^3)^(2/3) + (c*sin(theta)^3)^(2/3)

>> ans1=simplify(ans)
 
ans1 =
 
(c*cos(theta)^3)^(2/3) + (c*sin(theta)^3)^(2/3)
 
>> ans1=c^(2/3);
>> c=5;
>> for a=0.5:0.5:5
      x=0:0.05:5;
      y=(1-x/a)*(c^2-a^2)^(1/2);
    plot(x,y),hold on
   end
>> axis([0 5 -0.5 5.5])
>> theta=0:pi/50:pi/2;
>> x1=c*cos(theta).^3;y1=c*sin(theta).^3;
>> plot(x1,y1,'-b')
>> plot(x1,y1,'or')
>> grid on

没有意外的话我也许勉勉强强大概应该是基本都能理解了吧这题。biubiubiu~(づ￣3￣)づ╭❤～♥

♪炮弹抛物线包络线：

朕 乏 了
（又名：菜鸡不配拥有此题）

总结思考