- 1、本文档共20页,可阅读全部内容。
- 2、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
- 3、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 4、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
查看更多
精彩文档
实用标准文案
数值代数实验报告
Numerical Linear Algebra And Its Applications
学生所在学院:理 学 院
学生所在班级:计算数学10-1
学 生 姓 名:戈 东 潮
指 导 教 师:于 春 肖
教 务 处
2012年 12月
实验一
实验名称: Poisson方程边值问题的五点差分格式
实验时间: 2012年12月13日 星期四 实验成绩:
一、实验目的:
通过上机利用Matlab数学软件实现Poisson方程的边值问题的五点差分格式的线性代数方程组来认真解读Poisson方程边值问题的具体思想与方法,使我们掌握得更加深刻,并且做到学习与使用并存,增加学习的实际动手性,不再让学习局限于书本和纸上,而是利用计算机学习来增加我们的学习兴趣。
二、实验内容:
利用Poisson方程来解下列问题:
其中。
边值问题的解释是
用例题中模型问题的方法取N=5,10(也可以取N=20)列出五点差分格式的线性代数方程组
三、实验过程:
系数矩阵A的Matlab实现函数如下:
%文件名:poisson.m
function T=poisson(N)
for i=1:N
a(i)=4;
end
b=diag(a);
for j=1:N-1
b(j,j+1)=-1;
end
for j=2:N
b(j,j-1)=-1;
end
for i=1:N
I(i)=-1;
end
I=diag(I);
for i=1:N:N*N
for j=1:N
T(i+j-1,i:i+N-1)=b(j,:);
end
end
for i=1:N:N*N-N
for j=1:N
T(i+j-1,i+N:i+N-1+N)=I(j,:);
end
end
for i=1+N:N:N*N
for j=1:N
T(i+j-1,i-N:i-1)=I(j,:);
end
end
求解常数项b的Matlab实现函数如下:
%函数名:constant.m
function b=constant(N)
n=N+1;
h=1/n;
i=1;
for y=1/n:1/n:N/n
for x=1/n:1/n:N/n
f(i)=2*pi*pi*sin(pi*x)*sin(pi*y);
i=i+1;
end
end
b=h*h*f;
b=b;
四、实验结果(总结/方案)
在Matlab运行窗口输入
A=poisson(5) Enter
输出结果
A =
Columns 1 through 11
4 -1 0 0 0 -1 0 0 0 0 0
-1 4 -1 0 0 0 -1 0 0 0 0
0 -1 4 -1 0 0 0 -1 0 0 0
0 0 -1 4 -1 0 0 0 -1 0 0
0 0 0 -1 4 0 0 0 0 -1 0
-1 0 0 0 0 4 -1 0 0 0 -1
0 -1 0 0 0 -1 4 -1 0 0 0
0 0 -1 0 0 0 -1 4 -1 0 0
0 0 0 -1 0 0 0 -1 4 -1 0
0 0 0 0 -1 0 0 0 -1 4 0
0 0
您可能关注的文档
- 示范验收地总结材料.doc
- 市调网考的题目201620980.docx
- 试的题目及其问题详解--弹性力学与有限元分析报告.doc
- 试论《巴黎圣母院》地浪漫主义艺术特征.doc
- 试论1962年版越剧《红楼梦》地结构特征.doc
- 试析新课程实用标准指导下有效教学行为地基本特征.doc
- 试验大纲设计编写规定.doc
- 试验监理实施地研究细则.doc
- 试验仪操作规程指导书合集.doc
- 试验仪器操作规程82798.doc
- 2023届黑龙江哈尔滨第三中学高三二模英语试题 (解析版).docx
- 三校生对口升学考试计算机专题复习试题及答案解析.pdf
- 2023届江西省鄱阳县第一中学高三一模英语试题(原卷版).docx
- 2024届甘肃省高三下学期二模物理试题(解析版).pdf
- 2023届江苏省南京市等2地南京市教育科学研究所等2校一模英语试题(原卷版).docx
- 2023届江西省萍乡市芦溪中学高三一模英语试题(解析版).docx
- 2024届甘肃省高三二模语文试题(原卷版).pdf
- 2023届河南省郑州市高中毕业年级第二次质量预测英语试题(解析版).docx
- 2023届江西省抚州市金溪县第一中学等七校高三下学期第一次联考英语试卷(原卷版).docx
- 2023届吉林省吉林市普通中学高三上学期第二次调研测试英语试题(解析版).docx
最近下载
- 2024届高考语文写作分类训练-成功面面观.docx VIP
- 统一身份认证系统的设计与实现.docx VIP
- (高清版)-B-T 2659.1-2022 世界各国和地区及其行政区划名称代码 第1部分:国家和地区代码.pdf VIP
- ICN20531~32扫PWM恒流输出LED显示屏驱动芯片.pdf
- 2024世界献血者日20周年PPT课件(精选图文).pptx
- 妇科临床诊疗指南及操作规范.docx VIP
- 基于PLC的风力发电控制系统设计.doc VIP
- OH卡-亲密关系完整版.ppt
- 全套ISO27001信息安全文件一览表手册,程序文件,作业文件,表单.pdf VIP
- 2024第六届(2024年)“信用电力”知识竞赛活动总试题库-上(单选题汇总).docx
文档评论(0)