一：
%自由落體問題的高階偏微分方程 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:或者什么其他算法可以實現(xiàn)。
難道只能把初值t=0改成一個很小的數(shù)值嗎？
二：
類似問題dy(1)=y(2);   dy(2)=（右邊是關于y(1)和y(2)的函數(shù)且初值為零不含非齊次項）;
用龍格庫塔法y(1)始終為零怎么解決，例如我看的是黃懷緯博士論文第六章
確定m，n時，ξ系數(shù)皆為常數(shù)，圖像應該一直為零，不應該有波峰。
閉關學習中，謝謝好心大俠們。


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



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