一：
%自由落体问题的高阶偏微分方程 s=0.5*(s'')*t^2
function dy=myfun_1(t,y)
dy=zeros(2,1);
dy(1)=y(2);
dy(2)=2*y(1)/(t*t);
[x,y]=ode45(@myfun_1,[0,20],[0,0]);
plot(t,y(:,1),'-',t,y(:,2),'.')
问题1:初值t=0，而t在分母上；
问题2:或者什么其他算法可以实现。
难道只能把初值t=0改成一个很小的数值吗？
二：
类似问题dy(1)=y(2);   dy(2)=（右边是关于y(1)和y(2)的函数且初值为零不含非齐次项）;
用龙格库塔法y(1)始终为零怎么解决，例如我看的是黄怀纬博士论文第六章
确定m，n时，ξ系数皆为常数，图像应该一直为零，不应该有波峰。
闭关学习中，谢谢好心大侠们。


`UE@$[V9OG8Y86U]H6L{DDR.png



IC8M`SCVF][K{0CDZ0)_K_F.png 返回小木虫查看更多