MATLAB课程设计0920862班40号林战波.doc

太原工业学院 MATLAB 课程设计班级学号 092086240 姓 名 林战波 成 绩 指导老师 李莉 基于MATLAB的四阶龙格-库塔法1、 实践性质、目的与任务（一）性质MATLAB的名称源自Matrix Laboratory,1984年由美国Mathworks公司推向市场。它是一种科学计算软件，专门以矩阵的形式处理数据。MATLAB将高性能的数值计算和可视化集成在一起，并提供了大量的内置函数，从而被广泛的应用于科学计算、控制系统和信息处理等领域的分析、仿真和设计工作。MATLAB软件包括五大通用功能，数值计算功能（Nemeric）、符号运算功能（Symbolic）、数据可视化功能（Graphic）、数字图形文字统一处理功能（Notebook）和建模仿真可视化功能（Simulink）。其中，符号运算功能的实现是通过请求MAPLE内核计算并将结果返回到MATLAB命令窗口。该软件有三大特点，一是功能强大；二是界面友善、语言自然；三是开放性强。目前，Mathworks公司已推出30多个应用工具箱。MATLAB是一种以数值计算和数据图示为主的计算机软件，并包含适应多个学科的专业软件包，以及完善程序开发功能。本课程设计旨在要求学生熟悉Matlab的功能，掌握Matlab程序设计方法，培养学生的应用开发能力，为以后的毕业设计奠定一定的基础。（二）目的1、熟练掌握Matlab的基本编程方法，及其编程风格；2、熟练掌握Matlab常用函数的使用；3、与数值分析结合编程解决实际问题。（三）任务1、熟悉Matlab工作环境的使用，及其基本功能的设置，以及基本指令的使用；2、掌握Matlab数值数组及其矩阵操作的使用，以及基本绘图函数的使用规则；3、掌握Matlab字符串数组使用方法，掌握简单脚本文件和函数的编写规则；二、实践内容一金属片的初始温度为，其内部以的功率电加热，温度方程为： 其中温度单位为Kelvin，并有：（热传导率）（Stefan-Boltzmann常数）（表面积）（体积）（质量热容）（密度） （热交换系数）（发射率）用四阶龙格-库塔方法求解上的温度，取，作出图形，并给出相应的Matlab程序。三、实践步骤:首先使用ode45：function dT=lzb(t,T)dT=(3000-0.8*5.67*1.0e-8*0.25*(T4-2984)-30*0.25*(T-298)/(0.001*3000*900); t,T=ode45(lzb,0 600,298),grid,plot(t,T), grid然后使用四阶龙格-库塔方法求解：将程序命名为：lzb.mfunction f=f_x10(T)f=(3000-0.8*5.67*1.0e-8*0.25*(T4-2984)-30*0.25*(T-298)/(0.001*3000*900);y=298;h=0.1e-3;n=0;t=0;y_rec(1)=y;t_rec(1)=0;while t 10 n=n+1; k1=h*f_x10(y); k2=h*f_x10(y+k1/2); k3=h*f_x10(y+k2/2); k4=h*f_x10(y+k3/2); y=y+1/6*(k1+2*k2+2*k3+k4); t=(n-1)*h;y_rec(n+1)=y; t_rec(n+1)=t;endplot(t_rec,y_rec)四、心得体会 MATLAB是个编程工具，利用它你可以把理论上的数学知识很精确地表示出来，学习MATLAB要有数学基础，又要具备编程基础。 我在市里面学习编程，但通过这次课程设计，我发现自己存在一个问题，我考虑事情太单调，有点计算器那种感觉，一根筋，就不懂得融会变通，我想这也是我做编程一直存在的问题，轻易的不会去创新。但作为编程人员，要想在软硬件方便都很娇娇，这一点很重要，我会引以为戒的。参考文献：1. 张圣勤.MATLAB7.0实用教程.北京:机械工业出版社，2006.32. 张智星.MATLAB程序设计与应用.北京：清华大学出版社，20023.