计算物理与教学教改.pptVIP

概述  北京师范大学物理系 彭芳麟 一. 科学计算已经成为一种必不可少的科研方法 2004年诺贝尔物理学奖获得者David Gross 教授 在“The Future of Physics”的报告中将计算物理列为物理学面临的25个问题之一 定量计算越来越重要，计算机能否代替解析技术？应该如何调整对物理学家的训练？计算机何时能成为有创造力的物理学家？ 例1：用随机数计算π 例2：分子振动的半经典量子化　　　　　　　　　　　　 u=2.6:0.001:4; x=0.6; for j=1:150, x=u.*(x-x.^2); end for i=1:100 x=u.*(x-x.^2); plot(u,x,r.) hold on; End Mandelbrot集 Julia集 动画Mandelbrot分形.flv 分形艺术大奖\Tn0226a.tif 教学课件 二。 计算物理—科学计算—战略计算 计算物理: 计算替代实验 核武器研制 科学计算：高性能计算（硬件与软件）， 海量存储，计算科学 战略计算：信息高速公路，数字地球， 宽带网络 计算物理发展史.doc 国外名校都开《计算物理》 物理与天文教指委建议: 1998年,《计算物理》为研究生基础课 1999年, 重申这一建议 2004年,《计算物理基础》为本科必修课 计算物理基础规范　 计算物理基础是结合数值计算方法，利用计算机对物理问题进行数值计算、模拟和仿真的一门基础性学科。通过本课程的学习，使学生学会用计算机来描述物理图像、模拟物理问题，掌握用计算机处理物理问题的基本计算原理和必要的计算机知识。 　　　　　　　　　主 要 内 容 1．? 计算物理基本知识：计算物理发展史,科学计算基本原理,计算误差 2．? 计算基础：数值微分、积分和拟合,常微分方程与偏微分方程数值解,矩阵运算, 随机数和随机行走 3．? 基本物理问题计算：非谐问题（非谐振动的微扰）,非线性振动的付 里分析,混沌（钟摆问题）,多变量问题（静电、热传导、波动方程等偏微分方程数值解）4．? 数值模拟：蒙特卡罗模拟,分子动力学模拟 5．? 常用计算软件介绍：MATLAB等 教学时数：72-96学时，（其中课堂学时48-64，上机操作24-32学时） 教 改 要 求 如何学计算物理之一 计算物理 解题步骤 计算物理是综合训练，包括建模，算法，编程，物理解释。基础课要有基础性，实用性。 “进行计算的本领很难通过标准的大学生物理系课程安排来培养，因为它需要三方面训练（物理学，数值分析和计算机程序设计）的综合，很少有物理系学生毕业时知道如何进行计算；那些具有这种知识的学生通常是在独立工作的过程中，比如通过学习一项研究课题或者一篇毕业论文，而学到一些有限的一些技巧。” ——Steven E Konnin 如何学计算物理之二 软件编程优点：快捷，可视化 例：核磁共振图数据 物理场图示 如何学计算物理之三 浮点数运算与计算误差 为什么有：0.1+0.1+0.1+0.1+0.1+0.1+ 0.1+0.1+0.1+0.1+0.1+0.1 =9.999999999999999e-001 浮点数的表示 f 的取值限制数的精度，e的取值限制了数的范围 例： 在每个二进制区间， ，数按间隔 等距离排列 eps= 是两个浮点数的最大相对间距， eps/2是计算结果的最大相对误差。 注意： 计算机不是使用全部实数作计算，计算工具有缺陷！ 如何学之四 掌握编程技巧 两程序画相同图形,时间相差约600倍 u=[0,i]; subplot(3,3,1); plot(u) for k=2:8 m=u/3; uu=[m, i/3+m*(sqrt(3)*0.5+0.5i),m+i/3, ... 2i/3+m*(sqrt(3)*0.5-0.5i), m+2i/3]; axis([-0.5 0.5 0 1]); subplot(3,3,k); plot(uu) u=uu; end theta=pi/6; u =[0,0;0,1]; subplot(3,3,1) plot(u(:,1),u(:,2)) axis([-0.5,0.5,0,1]) for n=1:7 uuu=[]; for I=0:(lengt