Saturday, 2009-3- 21, 22:00:11
第10页/共39页
双精度数据结构
• 矩阵的表示
• 矩阵简单分析 F=inv(A) • 复数矩阵的表示
Saturday, 2009-3- 21, 22:00:11
第11页/共39页
特殊矩阵的输入
• 单位矩阵 A=eye(m,n) • 零矩阵 A=zeros(n,m) • 对角矩阵 A=diag([1,3,5,7,2]) • 均匀分布随机数矩阵 A=rand(n,m) • 多项式
Saturday, 2009-3- 21, 22:00:11
第4页/共39页
其他“高级”语言(如C)的局限性
• 语句繁杂 • 没有基本数学运算求解工具 • 程序冗长,容易出错,得出结果难以令人相信 • 各种各样的隐患 • 结果可视化效果不佳
Saturday, 2009-3- 21, 22:00:11
第5页/共39页
例1 矩阵特征值求解
• 用数学库EISPACK
• 用MATL AB: eig(A)
Saturday, 2009-3- 21, 22:00:11
第6页/共39页
例2 可能的隐患 Fibonacci数列
• 数学描述 • C语言实现
• 存在问题:数据结构的预定义 • MATL AB实现
第28页/共39页
• 函数 • MATL AB命令
• 检验 • 数据存储 • 导入
Saturday, 2009-3- 21, 22:00:11
第29页/共39页
极坐标图形绘制
• 绘制极坐标函数 • 仍然首先生成q 向量,再用点运算计算r • 调用polar函数绘图
• 绘制曲线
• 注意周期检验