资源描述：
《计算机要化学中应用.ppt》由会员上传分享，免费在线阅读，更多相关内容在教育资源-天天文库。
1、计算机要化学中应用1.数学工具：MatLab,Maple,Mathematic,SciLab,Octave2.数据绘图Origin,Sigmaplot,surf3D3.化学结构绘图(2D)ChemDraw,ISIS(3D)GaussianView,DSVisualizer4.公式编辑MathType5.Linux系统初步6.热力统计学第1章MATLAB语言程序设计基础薛定宇、陈阳泉著《高等应用数学问题的MATLAB求解》清华大学出版社，20042021/8/42高等应用数学问题的MATLAB求解东北大学信息学院MATLAB语言的简洁高效性MATLAB
2、语言的科学运算功能MATLAB语言的绘图功能MATLAB庞大的工具箱与模块集MATLAB强大的动态系统仿真功能MATLAB语言是当前国际上自动控制领域的首选计算机语言，也是很多理工科专业最适合的计算机数学语言，本课程以MATLAB为主要计算机语言。本章主要内容MATLAB程序设计语言基础基本数学运算MATLAB语言流程控制二维图形绘制三维图形绘制Chapter1MATLAB程序设计语言基础MATLAB语言的变量名规则由一个字母引导，后面可以为其他字符区分大小写AbcABc有效MYvar12,MY_Var12和MyVar12_错误的变量名12MyVar
3、,_MyVar12MATLAB的保留常量eps,i,j,pi,NaN,Inf,i=sqrt(-1)lastwarn,lasterr数值型数据结构双精度数值变量IEEE标准，64位(8字节)，11指数位，53数值位和一个符号位double()函数的转换其他数据类型uint8，常用于图像表示和处理，8位int8(),int16(),int32(),uint16(),uint32()符号型变量数据类型符号型，sym(A),常用于公式推导变量声明采用变精度函数求值【例】求出p的300位有效数字MATLAB的基本语句结构【例】赋值语句【例】矩阵赋值【例】复数矩
4、阵基本矩阵运算基本向量运算基本数论运算【例】对下面的数据进行取整运算-0.2765,0.5772,1.4597,2.1091,1.191,-1.6187【例】1856，12最大公约数、最小公倍数，因数分解【例】1-1000间质数1.MATLAB语言流程控制循环结构for结构while结构MatlabC++【例】用循环求解MatlabC++【例】求tics=0;forj=1:10000s=s+1/2^j+1/3^j;endtocticj=1:10000s=sum(1./2.^j+1./3.^j);toc转移结构【例】用循环求解求最大的ms=0;forj
5、=1:10000s=s+j;if(s>10000)break;endend;j开关结构2.二维图形绘制二维图形绘制基本语句构造向量:【例】直线，二次曲线x=[-10:0.5:10]y1=8*x+8y2=x.*x+2*x-13;plot(x,y1,x,y2);gridon【例】三角函数x=[0:0.05:2*pi];y1=sin(x);y2=cos(x);plot(x,y1,x,y2);gridonxlabel('angle');ylabel('value');title('sinandcos');【例】l=1.0x=0:0.01:l;n=2;f=sq
6、rt(2/l)*sin(n*pi*x)plot(x,f)gridonx=[-pi:0.05:pi]y=sin(tan(x))-tan(sin(x));plot(x,y);gridon【例】【例】H的1s径向分布函数r=0:0.1:5.0;Z=1.0;a0=1.0;R=1/(sqrt(pi))*(Z/a0)^1.5*exp(-1.0/a0*r);R2=R.*R;r2=4*pi*r.*r;D=R2.*r2;plot(r,D)gridonr=0:0.1:18.0;Z=1.0;a0=1.0;sigma=Z*r/a0;N=1/(2*sqrt(2));R=N*(
7、Z/a0)^1.5*(2-sigma).*exp(-sigma/2)R2=R.*R;r2=r.*r;D=R2.*r2;plot(r,D);gridon;【例】H的2s径向分布函数【例】绘制极坐标曲线theta=0:0.01:6*pi;rho=5*sin(4*theta/3);polar(theta,rho);theta=0:0.01:6*pi;rho1=5*sin(theta/3);polar(theta,rho1);其他二维图形绘制语句【例】用不同曲线绘制函数表示正弦曲线t=0:0.2:2*pi;y=sin(t);subplot(2,2,1),st
8、airs(t,y);subplot(2,2,2),stem(t,y);subplot(2,2,3),bar(
