- 博客(4)
- 收藏
- 关注
RSS订阅原创 解线性方程组的高斯消去法
clc;clear all;disp('此为列主元消去法');A=input('矩阵A=');b=input('右端项b=');eps=input('控制精度eps=');[n,n]=size(A);B=[A, b]; %%矩阵增广ra=rank(A); %求系数矩阵的秩rb=rank(B); %求增广矩阵的秩%判断方程解的情况if ra~=rb disp('系数矩阵...
2020-03-04 10:12:26
190
原创 matlab Newton-Cotes数值积分公式
function c=cotescoeff(n) %cotes系数syms tfor i=1:n+1 k=i-1; for j=1:n Ck(i)=(-1)^(n-k)/factorial(k)/factorial(n-k)/n*quad1(@(t)intfun(t,n,k),0,n) endendfunction [y,ck,ak]=newto...
2020-03-04 10:01:24
11571
2
原创 函数的最佳平方逼近
1、 试修改参考程序，使用正交多项式作为基底考虑使用勒让德多项式，修改程序如下：function S=square(a,b,n) %定义逼近函global i jif nargin<3 n=1;endX=zeros(n+1);for i=0:n for j=0:n X(i+1,j+1)=quad(@rho_phi,a,b); %求rho_phi...
2020-03-04 09:55:28
1428
原创 对函数做拉格朗日插值/牛顿插值
函数：函数：function y=f(x)y=5./(1+x.^2);end拉格朗日插值法函数：function [c,L]=lagran(x,y)w=length(x);n=w-1;l=zeros(w,w);for k=1:n+1 v=1; for j=1:n+1 if k~=j v=conv(v,poly(x(j)...
2020-03-04 09:19:23
1450
空空如也
空空如也
TA创建的收藏夹 TA关注的收藏夹
TA关注的人