电脑桌面
添加小米粒文库到电脑桌面
安装后可以在桌面快捷访问

北航数值分析计算实习1

北航数值分析计算实习1_第1页
1/6
北航数值分析计算实习1_第2页
2/6
北航数值分析计算实习1_第3页
3/6
《 数 值 分 析 》 计 算 实 习 题 目 1 10091013 劳 云 杰 一 、 算 法 设 计 方 案 根 据 提 示 的 算 法 , 首 先 使 用 幂 法 求 出 按 模 最 大 的 特 征 值 λ t1, 再 根 据 已 求 出 的 λ t1 用 带原 点 平 移 的 幂 法 求 出 另 一 个 特 征 值 λ t2, 比 较 两 个 λ 的 大 小 , 根 据 已 知 条 件 , 可 以 得 出 λ 1和 λ 501. 至 于 λ s, 由 于 是 按 模 最 小 的 特 征 值 , 使 用 反 幂 法 求 之 , 由 于 反 幂 法 需 要 解 线 性 方 程 组 ,故 对 矩 阵 进 行 Doolittle 分 解 。 再 通 过 带 原 点 平 移 的 反 幂 法 求 跟 矩 阵 的 与 数 最 接 近 的 特 征 值 。 对 非 奇 异 的 矩 阵 A, 根 据 条 件 数 定 义 , 取 λ t1/λ s 的 绝 对 值 , 两 个 特 征 值 在 之 前 步 骤 中均 以 求 得 。 由 于 对 矩 阵 进 行 了 Doolittle 分 解 , 所 以 矩 阵 的 行 列 式 detA 可 由 分 解 得 出 的 上 三 角 阵 U的 对 角 线 上 元 素 相 乘 求 得 。 为 了 使 A 的 所 有 零 元 素 都 不 存 储 , 使 用 书 本 25 页 的 压 缩 存 储 法 对 A 进 行 存 储 , 在 计 算 时通 过 函 数 在 数 组 C 中 检 索 A 中 元 素 即 可 。 由 于 A 是 501*501 矩 阵 , C 应取 为 5*501 矩 阵 。 由于 数 据 不 大 ,为 了 方 便起见,在 程 序中 取 502*502 矩 阵 或者502 向量,C 也取 为 6*502 矩 阵 。 程 序编写参考《 数 值 分 析 》 颜庆津著 和 [C 数 值 算 法 ].(美国)WILLIAM.H.PRESS.扫描版。 二、 全部源程 序 #include #include #define XS 1.0e-12//精度水平 void fz_a();//对 矩 阵 A 赋值 double js(int,int);//在 压 缩 矩 阵 中 检 索 A 的 元 素 double mf(double);//幂 法 double fmf(double);//反 幂 法 int lu(double);//Doolittle 分 解 int jfc(double[],double[]);//解 方 程 int max(int,int); int min(int,int); double (*u)[502]=n...

1、当您付费下载文档后,您只拥有了使用权限,并不意味着购买了版权,文档只能用于自身使用,不得用于其他商业用途(如 [转卖]进行直接盈利或[编辑后售卖]进行间接盈利)。
2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。
3、如文档内容存在违规,或者侵犯商业秘密、侵犯著作权等,请点击“违规举报”。

碎片内容

北航数值分析计算实习1

确认删除?
VIP
微信客服
  • 扫码咨询
会员Q群
  • 会员专属群点击这里加入QQ群
客服邮箱
回到顶部