|
|
EDA365欢迎您登录！
您需要 登录 才可以下载或查看，没有帐号？注册
x
本人菜鸟，刚入行学习matlab编程，参照温正编写的书籍《精通MATLAB科学计算》，在学习欧拉法中的向前欧拉法时，按照书中的例子，一直有问题产生，想请教一下大神们。
8 S+ D. v' c d9 B$ T软件：matlab R2017b，电脑：win10! o: c$ j6 \. G; ^
书中向前欧拉法函数文件如下：, G0 w2 U e/ T9 j5 ?! S
' g0 T8 e. u# G2 |* i: i4 r- function P=Eulerli1(x0,y0,b,h)
- %x0初值
- %y0为初值
- %b为x0取值区间的右端点
- %h为步长
- n=(b-x0)/h;
- X=zeros(int(n),1);
- Y=zeros(int(n),1);
- k=1;
- X(k)=x0;
- Y(k)=y0;
- for k=1:n
- X(k+1)=X(k)+h;
- Y(k+1)=Y(k)+h*(X(k)-Y(k));
- k=k+1;
- end
- y=X-1+2*exp(-X);
- plot(X,Y,'mp',X,y,'b-')
- grid
- xlabel('自变量 X'),ylabel('因变量')
- legend('h=0.075 数值解','精确解')
- jwY=y-Y;
- xwY=jwY./y;
- k1=1:n;
- k=[0,k1];
- P=[k',X,Y,y,jwY,xwY];
1 b9 H5 u/ i& y+ [! w
& G* y& M6 s0 {% r1 X( G
! b+ d1 g. E) t* P6 f* a2 L1 m+ U
红色部分就是一直出问题的地方，我上网看了一下，zeros函数里面必须是整数，可是我应该如何修改这个函数文件呢？, ?( Q7 P/ D: P
自己也看了其他人的向前欧拉法，比如：' z; F1 t* L) S0 N+ }
- % forward euler method
- function [ x, y ] = euleRF ( f, y0, a, b, n )
- y ( 1 ) = y0 ;
- h = ( b - a ) / n;
- x = a : h : b;
- for i = 1 : n
- y ( i +1 ) = y ( i ) + h * feval ( f, x ( i ), y ( i ) ) ;
- end;
+ O# q& t f7 ?0 O5 n
6 c( @+ U6 L' T, A; P关键第一种书上的还有图形出来，但是可能函数的一些编写理由自己不是很清楚，麻烦大神指点一下，非常感谢！* {/ Q2 l% ^& h% T
|
|
/1 
发表于 2020-6-1 15:10
收藏
淘帖
支持！
反对！
发表于 2020-6-1 15:58