- 博客(3)
- 收藏
- 关注
RSS订阅原创 基于MATLAB实现的龙贝格求积分算法
function [int] =integRob(xi,a,b)tk=zeros(1,100);sk=zeros(1,100);ck=zeros(1,100);rk=zeros(1,100);tk(1)=(b-a)*(0.5*polyval(xi,a)+0.5*polyval(xi,b));for n=2:20 sum=0; for m=1:2^(n-2) sum=sum+polyval(xi,a+(2*m-1)*(b-a)/2^(n-1)); end
2021-01-11 18:54:33
3249
4
原创 基于MATLAB实现最小二乘法的曲线拟合
function [a] = fitting(t,x,order)k=zeros(length(x),order+1);for i=1:length(x) for j=1:order+1 k(i,j)=t(i)^(j-1); endendy=k'*x';k=k'*k;a=sorm(k,y,1.8);%超松弛迭代法求解正规方程a=fliplr(a);end输出的a向量即为拟合方程曲线的多项式系数向量...
2021-01-11 18:51:10
1191
原创 基于MATLAB实现固定步长的四阶龙格库塔方法求解常微分方程初值问题
代码构造函数function [yn,xn] = RKorder4(f,a,b,h,y0)%a，b为上下界；h为步长；y0为初值。%返回yn、xn向量，并绘制曲线yn=zeros(1,((b-a)/h)+1);yn(1)=y0;xn=a:h:b;for i=1:((b-a)/h) k1=h*f(xn(i),yn(i)); k2=h*f(xn(i)+h/2,yn(i)+k1/2); k3=h*f(xn(i)+h/2,yn(i)+k2/2); k4=h*f(xn(
2021-01-11 18:50:52
3682
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人