2020年1月30日
摘要: 在下面的这段代码中,包含了高斯-勒让德、高斯-切比雪夫、以及拉盖尔和埃尔米特型求积公式,它们分别对应了不同的被积积分型 1.代码 %%高斯型求积公式 %%Y是函数表达式,interval是求积区间,n是求积阶数 %%对于求一般形式的非反常积分,可用勒让德型, %%对于求形如f(x)/sqrt(1-x 阅读全文
posted @ 2020-01-30 17:15 谷梁天 阅读(4576) 评论(0) 推荐(0) 编辑
摘要: 此段代码是基于辛普森公式的积分计算方法 1.代码 %%复合辛普森求积公式 %%Y是数值(attribute=0)或具体表达式(attribute=1),interval是求积区间,n是精度(如果是数值,则为数值长度-1) function CSQF = Compound_Simpson_quadra 阅读全文
posted @ 2020-01-30 16:34 谷梁天 阅读(3702) 评论(0) 推荐(0) 编辑
摘要: 这段代码实现的是最一般的数值积分法——梯形求积法,积分值的准确依赖于所取精度大小 1.代码 %%复合梯形求积公式 %%Y是数值(attribute=0)或具体表达式(attribute=1),interval是求积区间,n是精度(如果是数值,则为数值长度-1) function CTQF = Com 阅读全文
posted @ 2020-01-30 16:18 谷梁天 阅读(1737) 评论(0) 推荐(1) 编辑
摘要: 此段代码是牛顿- 柯特斯数值积分法,代码如下: 1.代码 %%牛顿-柯特斯公式(此方法对于8阶以下是有效的,8阶以上误差将非常大) %%interva为求积区间,Y随attribute变化(0或1)而对应不同选项(已知X对应的数值 或 表达式),n为步数 function NCF = Newton_ 阅读全文
posted @ 2020-01-30 16:11 谷梁天 阅读(2066) 评论(0) 推荐(0) 编辑
2020年1月16日
摘要: 本段代码可实现OLS法的线性回归分析,并可对回归系数做出分析 1.代码 %%OLS法下的线性回归 function prodict = Linear_Regression(X,Y) x = sym('x'); n = max(size(X)); %%定义画图窗格属性 h = figure; set( 阅读全文
posted @ 2020-01-16 18:08 谷梁天 阅读(1894) 评论(0) 推荐(0) 编辑
2019年12月30日
摘要: 这段代码可以给出常用的4类正交多项式的具体表达式,后续将给出可自定义的任意正交多项式代码 %%正交多项式 %%此函数包括勒让德正交多项式,切比雪夫正交多项式(两类),拉盖尔正交多项式,埃尔米特正交多项式,输入项数应从1开始 %%n是多项式的项数,n>=0,type是类型,分为Legendre、Che 阅读全文
posted @ 2019-12-30 17:23 谷梁天 阅读(1670) 评论(0) 推荐(0) 编辑
摘要: 1.代码 %%列主元消去法 function ECPE = Elimination_of_column_pivot_entries(M,b) global n; [n,n] = size(M); B =[M,b]; R_A = rank(M);R_B = rank(B); if R_A ~= R_B 阅读全文
posted @ 2019-12-30 14:00 谷梁天 阅读(1151) 评论(0) 推荐(0) 编辑
摘要: 1.代码 %%LU分解法 function LUDM = LU_Decomposition_method(A,b) global n;global B;global U;global L;global M; [n,n] = size(A); B = [A,b]; R_A = rank(A);R_B 阅读全文
posted @ 2019-12-30 13:58 谷梁天 阅读(1801) 评论(0) 推荐(0) 编辑
摘要: 1.代码 %%雅可比迭代法(此迭代法对于病态矩阵的解不理想) %%线性方程组M*X = b,M是方阵,X0是初始解向量,epsilon是控制精度 function JIM = Jacobian_iteration_method(M,b,X0,epsilon) [m,n] = size(M); d = 阅读全文
posted @ 2019-12-30 13:52 谷梁天 阅读(3589) 评论(0) 推荐(0) 编辑
摘要: 1.代码 %%高斯-塞得勒迭代法 %%线性方程组M*X = b,M是方阵,X0是初始解向量,epsilon是控制精度 function GSIM = Gauss_Seidel_iterative_method(M,b,X0,epsilon) [m,n] = size(M); d = diag(M); 阅读全文
posted @ 2019-12-30 13:48 谷梁天 阅读(975) 评论(0) 推荐(0) 编辑