怎样算高斯烟团模型的算法,编程的代码

解放军生活杂志 2024-11-17 17:41:34
最佳回答
% gauss.mfunction x = gauss(a,b)%the sizes of matrices a,b are supposed to be na x na and na x nb.%th** function solves ax = b by gauss elimination algorithm.na = size(a,2); [nb1,nb] = size(b);if nb1 ~= na, error('a and b must h**e compatible dimensions'); endn = na + nb; ab = [a(1:na,1:na) b(1:na,1:nb)]; % augmented matrixepss = eps*ones(na,1);for k = 1:na %scaled partial pivoting at ab(k,k) by eq.(2.2.20) [akx,kx] = max(abs(ab(k:na,k))./max(abs([ab(k:na,k + 1:na) epss(1:na - k + 1)]'))'); if akx < eps, error('singular matrix and no unique solution'); end mx = k + kx - 1; if kx > 1 % row change if necessary tmp_row = ab(k,k:n); ab(k,k:n) = ab(mx,k:n); ab(mx,k:n) = tmp_row; end % gauss forward elimination ab(k,k + 1:n) = ab(k,k+1:n)/ab(k,k); ab(k,k) = 1; %make each diagonal element one for m = k + 1: na ab(m,k+1:n) = ab(m,k+1:n) - ab(m,k)*ab(k,k+1:n); %eq.(2.2.5) ab(m,k) = 0; endend%backward substitution for a upper-triangular matrix eqation% h**ing all the diagonal elements equal to onex(na,:) = ab(na,na+1:n);for m = na-1: -1:1 x(m,:) = ab(m,na + 1:n)-ab(m,m + 1:na)*x(m + 1:na,:); %eq.(2.2.7)end例:>> a = [0 1 1;2 -1 -1;1 1 -1]; b = [2 0 1]';>> x = gauss(a,b)x = 1.00000000000000 1.00000000000000 1.00000000000000 baiduyixiajiukeyia 20210311
汇率兑换计算器

类似问答
  • 科学计算的程序编写和通常所说的码农的编程有多大的区别?
    • 2024-11-17 17:43:50
    • 提问者: 未知
    科学计算程序的编写通常需要利用数学函数或运算库,除此之外它和通常意义上的码农的工作有何差异?平常仅…
  • 关于capm模型的一道计算题
    • 2024-11-17 04:55:59
    • 提问者: 未知
    基金收益率=10%*7%90%*15%14.2%
  • 跪求简单计算器的j**a代码,要入门级的编程~~
    • 2024-11-17 01:30:30
    • 提问者: 未知
    import j**a.awt.*;import j**a.awt.event.*;import j**ax.swing.*;import j**a.util.*;public class counter { public static void main(string[] args) { counterframe frame = new counterframe(); ...
  • 一箭穿心3b编程代码
    • 2024-11-17 05:07:54
    • 提问者: 未知
    //将有序数组a[]和b[]合并到c[]中void memeryarray(int a[], int n, int b[], int m, int c[]){int i, j, k;i = j = k = 0;while (i < n && j < m){if (a[i] < b[j])c[k++] = a[i++];elsec[k++] = b[j++]; }while (i ...
  • 网络团购模式的会计核算流程是怎么样的
    • 2024-11-17 15:08:09
    • 提问者: 未知
    华律网小编拟从团购网站和商户的角度,结合会计核算的相关原则,举例分析网络团购业务的具体会计核算流程,希望这些内容能对大家有所帮助。一、网络团购的一般交易流程及...
  • 单位工程概算的编制
    • 2024-11-17 15:03:25
    • 提问者: 未知
    1.单位建筑工程概算的编制依据 (1)设计任务书内容随着建设项目的性质而不同,一般包括建设目的、建设规模、建设根据、建设布局、建设内容、建设进度、投资估算、新产品方案和原料来源等。(2)设计文件设计文件的内容包括初步设计图纸、设计说明书、总平面布置图、工程项目一览表、设备表、...
  • 科学计算的程序编写和通常所说的码农的编程有多大的区别?
    • 2024-11-17 14:45:46
    • 提问者: 未知
    一切以物理研究维准绳,你能提出一个新公式,甚好;能在计算中发现新现象,不错;把程序用面向对象…
  • 计算机编程群
    • 2024-11-17 18:29:01
    • 提问者: 未知
    你刚上大学?还是想自学编程?如果刚上大学,跟着课程安排学就ok。如果想自学编程,就先从c语言开始吧。网上很多c语言的学习资料,也有很多c语言学习的书籍。
  • 计算机编程的魅力在哪?
    • 2024-11-17 08:03:16
    • 提问者: 未知
    帮你干活。当时是学按键精灵,去打游戏。后来是学python,去处理文件…
  • 美国李维斯的尺码怎么算啊
    • 2024-11-17 05:20:29
    • 提问者: 未知
    尺码 腰围(cm)你的腰围2,4尺=80 应该穿w32的裤子 w28-74 w30-78 w32-83 w34-90 w36-96
汇率兑换计算器

热门推荐
热门问答
最新问答
推荐问答
新手帮助
常见问题
房贷计算器-九子财经 | 备案号: 桂ICP备19010581号-1 商务联系 企鹅:2790-680461

特别声明:本网为公益网站,人人都可发布,所有内容为会员自行上传发布",本站不承担任何法律责任,如内容有该作者著作权或违规内容,请联系我们清空删除。