Matlab在物理上的应用举例
MATLAB在大学物理课程及实验教学中的应用(1)
MATLAB在大学物理中的应用形式十分相似,故用MATLAB来解算问题要比用C,FORTRAN等语言完成相同的事情简捷得多,并且MATLAB也吸收了像Maple等软件的优点,使MATLAB成为一个强大的数学软件。
在新的版本中也加入了对C,FORTRAN,C++,JAVA的支持。
可以直接调用,用户也可以将自己编写的实用程序导入到MATLAB函数库中方便自己以后调用,此外许多的MATLAB爱好者都编写了一些经典的程序,用户可以直接进行下载就可以用。
20世纪70年代,美国新墨西哥大学计算机科学系主任Cleve Moler为了减轻学生编程的负担,用FORTRAN编写了最早的MATLAB。
1984年由Little、Moler、Steve Bangert 合作成立了的MathWorks公司正式把MATLAB推向市场。
到20世纪90年代,MATLAB已成为国际控制界的标准计算软件。
布朗运动表述的是悬浮在流体中的微粒受到流体分子与粒子的碰撞而发生的不停息的随机运动,其运动的剧烈程度受温度影响.这一运动比较抽象,我们可以用MATLAB 仿真,将其可视化.模拟效果如图 1 所示.图1 布朗运动的模拟效果图部分程序如下:axes(handles.axes1);global k;n = 50;x = rand(n,1)-0.5;y = rand(n,1)-0.5;h = plot(x,y,'.');axis([-1 1 -1 1]);axis square;grid off;set(h,'EraseMode','Xor','MarkerSize',24);if k >= 0;for i = linspace(1,10,10000)drawnowx = x + k /10000.* randn(n,1);y = y + k /10000.* randn(n,1);set(h,'xdata',x,'ydata',y);endelseinput('Error');end力学中的简谐振动和阻尼振动也可以用MATLAB 将其化为可视图形.如图 2 所示:2. 2 光学现象如干涉和衍射是光学中的重要物理现象,我们通过MATLAB 来实现衍射现象的可视化.单缝衍射程序:·40·李仲,董松,炊万年,沈武鹏:MATLAB 在大学物理课程及实验教学中的应用图2简谐振动规律和阻尼振动规律模拟效果图图3单缝衍射效果图functionpushbutton1_Callback(hObject,eventdata,handles)a = -2* pi:0.0001* pi:2* pi;P = (1 -sinc(a)).^2;plot(a,P)lgray = zeros(256,3);for i = 0:255lgray(i + 1,:)= (255 -i)/255;endimagesc(P)colormap(lgray)可视化效果如图 3 所示.还可根据单缝宽度来获得不同的衍射图像.2. 3 电磁学现象等量异种电荷电场线的[x ,y ]= meshgrid ( - 2: 0. 1: 2,- 2: 0. 1: 2) ; % 以 0. 1 为步长建立平面数据网格 z = 1. / sqrt ( ( x - 1) . ^2 + y. ^2 + 0. 01)- 1. / sqrt ( ( x + 1) . ^2 + y. ^2 + 0. 01) ; % 写电势表达式 [px ,py ]= gradient ( z ) ; % 求电势在 x ,y 方向的梯度即电场强度 contour ( x ,y ,z ,[- 12,- 8,- 5,- 3,- 1,- 0. 5,- 0. 1,0. 1,0. 5,1,3,5,8,12]) % 画等势线·41·青海民族大学学报(教育科学版)hold on % 作图控制quiver(x,y,px,py,'k')% 画出各点上电场的大小和方向等量异种电荷电势线描绘:[x,y]= meshgrid(-5:0. 2:5,-4:0. 2:4);% 建立数据网格z = 1. / sqrt((x -2). ^2 + y. ^2 + 0. 01)-1. / sqrt((x + 2). ^2 + y. ^2 + 0. 01);% 电势的表达式mesh(x,y,z)% 三维曲面绘图图4 等量异种电荷电场线及电势线描绘3 物理实验数据处理在物理实验中,通常对数据采用的是手工处理方法,常见的有列表法,作图法等,这些方法往往速度慢,效率低,而且过于繁琐.如最小二乘法是采用数理统计的方法来处理实验数据的,相比于图解法,用该方法处理实验数据更科学、更可信.但由于该方法复杂,计算量大,学生们很难顺利完成,而运用MATLAB 可以精确实现.它克服了最小二乘法计算量大的缺点,又简便易懂,很容易为学生所掌握,同时能拟合出较准确的曲线,轻松实现数据可视化.如在电阻的伏安特性曲线的绘制中,已知测得的电流、电压值分别为0. 662,0. 712,0. 782,0. 841,0. 931,0. 988A0. 1720,0. 1846,0. 2024,0. 2182,0. 2364,0. 2560V;可以用MATLAB 所提供的数据拟合多项式函数polyfit 和评估多项式函数polyval 来进行曲线拟合,并且可以计算误差平方和,方便的得到伏安特性曲线、电阻值和计算误差. 程序如下:x = [0. 1720,0. 1846,0. 2024,0. 2182,0. 2364,0. 2560];y = [0. 662,0. 712,0. 782,0. 841,0. 931,0. 988];·42·李仲,董松,炊万年,沈武鹏:MATLAB 在大学物理课程及实验教学中的应用p1 = polyfit(x,y,1);% 一次多项式拟合p2 = polyfit(x,y,9);% 九次多项式拟合q1 = polyval(p1,x);q2 = polyval(p2,x);s1 = sum((y -q1). ^2);% 一次多项式误差平方和s2 = sum((y -q2). ^2);% 九次多项式误差平方和plot(x,y,'* ')hold onplot(x,q1,'r')hold onplot(x,q2,'b:o')grid onp1,q1,s1,p2,q2,s2MATLAB 绘制的伏安特性曲线如图5 所示.拟合次数越高,曲线越精确,可以看出,一次拟合结果为:R = 3.9653 ±(2.5345e -004)Ω4 结语应用MATLAB 设计程序对物理现象规律进行模拟仿真,实现了物理现象、规律描述的可视化.通过物理实验数据处理程序的设计,可推动大学物理实验教学现代化的改革.这项实践活动不仅有助于增强学生对物理课程学习及MATLAB 软件应用开发的兴趣,还可实际应用于大学物理课程及实验教学的活动中.学以致用,是提高大学课程教学效果的有效途径和手段.参考文献:[1]苏金明,张莲花,刘波.MATLAB 工具箱应用[M].北京:电子工业出版社,2002.[2]彭芳麟.数学物理方程的Matlab 解法与可视化[M].北京:清华大学出版社,2004.[3]熊万杰.MATLAB 用于物理教学[J].物理通报,2004,(2)16 -19.·43·感谢您试用AnyBizSoft PDF to Word。
MATLAB在物理教学中的应用
分别画出 t = 1 S , t =1 . 2 5 s , t = 1 . 5 s 时刻的波形. 分析 : 将 题 中波 动 方 程表 达 式 与 平 面简 谐 波 标
作者 简介 : 田
澍( 1 9 8 5~) , 女, 山西兴县人 , 助教 , 研 究方向为发光材料
对 于 很 多学 生来 说 物理 是 一 门较 难 的学 科 , 一 个 主要 的原 因是物 理 中包 含 大 量 复 杂 的计 算 . 枯 燥 的理论 计 算 、 繁琐 的 数 据 处 理很 难 激 发学 生 的 学 习 兴趣 , 同时过 度依 赖 书 本 会 导 致 想 象 力 不太 好 的 学
田 澍 , 白 心 爱
( 吕梁学院 物理系 , 山西 离石 0 3 3 0 0 0 )
摘
要 :M A T L A B具有 强大的计算和绘 图功 能, 语言 简洁容 易掌握. 通过 对麦克 斯 韦速 率分布 函数 、 不 同时刻
平面 简谐波 的波动方程和 电偶极子 电 势函数 的作 图和 分析 , 表 明 MA T L A B可以极 大地 简化物 理计 算和绘 图工作
第 3卷 第 2期
V o 1 . 3 No . 2
吕梁 学 院学 报
J o u r n a l o f L t i l i a n g U n i v e r s i t y
2 0 1 3年 4月
Apr . 2 01 3
・
物理 学研 究 ・
MA T L A B在物理教学 中的应用
y =0 . 0 6 c o s ( 1 0  ̄ r t 一 8  ̄ r x )
数 k= 1 . 3 8 × 1 0 。 J / K . 麦克斯韦速率分布函数是 比 较复杂的 , 从式 中不能够直接观察 出其与温度 的关 系. 所 以笔 者用 MA T L A B进行 了 编程 作 图 , 取T 1=
适合用matlab解决的经典物理例题
适合用Matlab解决的经典物理例题在物理学领域,经典物理例题一直是学习和研究的重要内容。
而Matlab作为一种强大的数学软件,非常适合解决各种物理问题。
本文将从力学、电磁学和热力学等多个方面,选取一些经典的物理例题,通过Matlab进行分析和求解,展示Matlab在解决物理问题时的强大用途。
1. 简谐振动问题简谐振动是物理学中一个重要的模型,涉及到弹簧振子、单摆等问题。
通过Matlab可以很方便地求解简谐振动的运动规律。
对于弹簧振子的运动方程,可以通过Matlab进行数值模拟,得到振动的周期、频率、位移等参数,从而更好地理解简谐振动的特性。
2. 电场问题在电磁学中,电场是一个重要的研究对象。
通过Matlab可以很容易地分析不同形状的电荷分布所产生的电场分布。
可以通过Matlab计算出点电荷、均匀带电细棒等情况下的电场分布,并绘制出电场线图,直观地展现电场的分布规律。
这样的分析对于理解电场的性质和相互作用具有重要意义。
3. 热传导问题热传导是热力学研究的一个重要方面,涉及到导热方程的求解和热量分布的分析。
通过Matlab可以对不同材料和形状的热传导问题进行数值模拟和求解。
可以通过Matlab计算出棒状材料中的温度分布随时间的演化,从而得到材料的热传导性能。
这样的分析对于工程实践中的热设计和材料选型具有重要指导意义。
4. 万有引力问题在力学中,万有引力是一个经典的例题,涉及到行星轨道、卫星运动等问题。
通过Matlab可以很方便地进行万有引力场下的物体运动模拟。
可以通过Matlab计算地球和月球的引力作用下的月球轨道,从而揭示天体运动的规律和特性。
这样的模拟对于探索宇宙中天体运动规律具有重要帮助。
总结回顾:通过以上例题的分析,我们不仅了解了Matlab在经典物理例题中的应用,也可以发现Matlab在解决物理问题时的便捷和高效。
当然,实际物理问题可能具有更多的复杂性和多样性,需要结合理论分析和实验数据进行综合研究。
Matlab在物理学中的应用
Matlab在物理学中的应用Matlab是一种广泛应用于工程计算及数值分析领域的功能强大的计算机高级语言,它集科学计算、图象处理于一身,并提供了丰富的图形界面设计方法。
它的特点是语法结构简单,数值计算高效,图形功能完备,特别适合于非计算机专业的编程人员完成日常数值计算、科学实验数据处理、图形图象生成等通用性任务时使用,因而在统计、信号处理、自动控制、图象处理、人工智能及计算机通信等领域得到了广泛应用。
现在,我们将介绍Matlab在物理学中的一些应用。
首先,我们来介绍Matlab在动力学中的应用。
动力学的最基本的公式为F=m*a,其中m是物体的质量,a是物体的加速度。
如何验证上述公式的正确性呢?我们设计了一些实验,得到了在质量一定的情况下a与F的一些实验数据。
而根据实验数据来确定a与F的关系,可以用Matlab中的拟合曲线来完成。
实例如下:我们测得当m=0.45kg时,a与F的对应关系如下表所示:我们在Matlab中输入以下指令:F=[0.3726 0.7453 1.1157 1.4994 1.8566 2.2271];a=[0.0845 0.169 0.253 0.34 0.421 0.505];p=polyfit(F,a,1);i=0:0.1:2.5;q=polyval(p,i);plot(F,a,'r.',i,q,'b-')axis([0 2.5 0 0.6])xlabel('F 单位:N');ylabel('a 单位:m/s^2');我们可以得到如下视图:上图中蓝色的直线是F与a线性拟合曲线,我们看到实验数据点(F,a)几乎全分布在拟合曲线上,由此我们可以大致得到结论,F与a呈线性关系。
延长拟合曲线可发现其大致过原点,由此我们可以推断F与a成正比关系,直线的斜率即为物体质量的倒数。
再次,我们来看看Matlab在绘制简谐运动图像上的应用。
MATLAB在大学物理课程教学中的应用-以多普勒效应为例
MATLAB在大学物理课程教学中的应用-以多普勒效应为例摘要:运用MATLAB仿真大学物理实验,能够帮助学生更好的对物理概念和规律进行理解和掌握,同时有效提升学生运用科学计算的能力,极大的提高物理教学的效果。
本文以多普勒效应相关内容为例,进行MATLAB仿真模拟分析。
关键词:多普勒效应 MATLAB 仿真分析在计算机仿真日益盛行的今天,作为一种重要的科学工具,计算机已经广泛渗透到人们生活的方方面面。
随着计算机仿真技术的发展,利用仿真虚拟技术展示客观物理现象,在各行各业均得到了广泛应用,逐渐成为社会发展进步不可或缺的手段。
在高校物理教学领域里,信息技术与教学相结合所带来的教育信息化已经成为当前高校物理教育改革的热门研究课题和教育发展的必然趋势。
一方面,利用计算机仿真技术优势,将枯燥难懂的物理问题变成一幅幅生动的画面,增强了教学内容的直观性,生动性,扩展了教学内容,开拓了学生视野,促进了学生对知识的理解和掌握。
另一方面,利用MATLAB仿真大学物理实验可以大大激发学生对物理知识的好奇心和求知欲,强化了学生继续探索的动机,激发了学生的创新意识,同时也极大地提高大学物理课程教学效果。
另外,用 MATLAB 制作的软件有极其丰富的表现内涵和巨大的表现能力,能够具体形象地再现各种实验概念,有效地揭示事物的本质和内在联系,将它应用于课堂教学,极大地扩展教学空间,化繁为简,变难为易,使学生对教学内容更容易理解和掌握。
本文就以物理课程中的多普勒效应为例进行仿真模拟分析,研究接收者接收到的频率变化规律。
我们知道,当一辆汽车在我们的身旁疾驰而过的时候,车上喇叭的音调呈现出从高到低的突然变化过程。
同样的,当我们在铁路旁听列车的汽笛声也能够发现,列车迅速迎面而来时音调较静止时高,而列车迅速离去时则音调较静止时低。
这种由于波源和接收者相对运动而出现接收者接收频率变化的现象,称之为多普勒效应。
多普勒效应最早由奥地利物理学家多普勒在1842年首先发现。
Matlab技术在物理建模和数值计算中的应用
Matlab技术在物理建模和数值计算中的应用近年来,Matlab成为了物理学研究和数值计算的重要工具。
其强大的数据处理和可视化能力以及灵活的编程语言,为物理建模和数值计算提供了极大的便利。
本文将探讨Matlab技术在物理建模和数值计算中的应用,并展示其在不同领域中的实际案例。
一、物理建模物理建模是物理学研究的关键步骤,它通过建立数学模型来描述和解释自然现象。
Matlab作为一种高效的数学计算工具,可以帮助研究人员快速地进行物理建模,并对模型进行分析和优化。
例如,在光学领域,研究人员经常使用Matlab来建立光传输模型。
他们可以基于光的传输方程和介质的光学性质,编写Matlab程序来计算和可视化光在不同介质中的传播路径和强度分布。
这种建模和仿真的方法可以帮助研究人员理解光的传输机制,并为光学器件设计提供指导。
在力学领域,物理建模可以用于模拟和分析复杂的物体运动。
例如,研究人员可以使用Matlab的刚体运动模型来研究机器人的动力学特性。
他们可以通过在Matlab中定义机器人的几何结构和关节运动,以及描述力和力矩的方程,来模拟机器人在不同任务下的运动轨迹和力学响应。
这种建模方法可以为机器人控制算法的开发和优化提供测试平台。
二、数值计算数值计算是物理学研究中的重要工具,它通过数值方法对物理问题进行近似求解。
Matlab提供了丰富的数值计算库和函数,使得研究人员能够高效地进行数值计算和分析。
在量子力学领域,研究人员常常使用Matlab进行量子系统的数值模拟和求解。
他们可以编写Matlab程序来求解薛定谔方程,从而获得量子系统的波函数和能谱。
这种数值求解的方法在研究纳米尺度物理系统以及量子信息领域具有重要的应用价值。
在流体力学领域,数值计算在流体流动的模拟和分析中起到关键作用。
研究人员可以使用Matlab的流体力学工具箱来建立流体模型,并通过数值方法求解雷诺方程等流体动力学方程。
这种基于Matlab的数值模拟方法可以帮助研究人员研究流体流动的行为、优化流体系统的设计,并解决与流体力学相关的实际问题。
matlab物理案例分析
Matlab及其应用
24
线型和标记
Plot命令中,在每一对数组数据后面,给plot一个 附加参量,就可以指定所要的颜色和线型。 >>plot(x,y,’rd’,x,z,’k--’,x,y,’b-.’,x,z,’co’) b蓝色,c青色,g绿色,k黑色,m紫红色,r红色,w白 色,y黄色 加格线和标注,grid on命令在当前图形的单位标 记处加格线。grid off命令取消格线。 用text命令可以在图形的任意指定位置增加标记和 其他文本信息,格式是text(x,y,’string’).
Matlab及其应用
2
Matlab的学习方法
必须做大量的练习,熟悉其中的函数 联系和自己的课题相关,深刻的体会 多看帮助文件,又一本好的参考书 记住:Matlab可以做很多事情
Matlab及其应用
3
1.1 Matlab的出现
70年代中期,Cleve Moler和他的同事开发了 LINPACK和EISPACK的Fortran子程序库 70年代末期,Cleve Moler 在新墨西哥大学给学 生开线性代数,为学生编写了接口程序,这程序取 名为MATLAB,即MATrix LABoratory 1983年春天,工程师John Little与Moler、 Steve Bangert一起开发了第二代专业版MATLAB 1984年,MathWorks公司成立,MATLAB正是推 向市场。
22
Matlab及其应用
五、数据的可视化
二维图形
在二维曲线的绘制命令中,最重要、最基本的命令 时plot,其他许多特殊绘图命令都以它为基础形成。 1 plot(x,’s’) x为实向量时,以该向量元素的下标为横坐标、元素值 为纵坐标画出一条连续曲线。X为实矩阵时,则按列绘 制每列元素值相对其下标的曲线,图中曲线数等于X阵 列数。s用来指定线型。(正方形)
MATLAB在高中物理学习中的应用
教育研究课程教育研究110学法教法研究1.MATLAB 在力学分析中的应用如图1,质量为m 的小球,从静止开始沿半径为R 的1/4圆弧光滑无摩擦轨道MN 滑下,重力的瞬时功率如何变化?图1 小球下降过程示意图1.1 建模小球下降过程中的瞬时功率:(θ为F 和v 的夹角),通过受力分析,最后可解得重力瞬时功率的表达式为:(1)通过式(1),学生的一般分析是:重力在圆弧轨道上的分离逐渐减小,但小球下滑的速度会逐渐加大,但两者乘积的结果如何难以判断。
通过复杂数学计算可以得出P 的极值,进而可推断出瞬时功率先增大后减小。
但功率随角度θ具体变化的情况学生根据已有知识却无法很快做出判断。
1.2 程序及结果运行结果如图2所示,从图中很容易看出随着角度变化,瞬时功率先增加后减小,在35度左右发生转折。
图2 功率随角度变化曲线2.MATLAB 在电学分析中的应用如图3所示,倾角为30度的直角三角形底边长2l ,且处于水平位置,斜面为光滑绝缘导轨,现在底边中点O 处固定一正电荷Q,让一质量为m 的带负电的点电荷q 从斜面顶端A 处释放沿斜面滑下(不脱离斜面)。
现测得它滑到B 点在斜边上的垂足D 处的速度为v ,问质点的运动过程如何?图3 电荷受力分析图2.1 建模我们可以分析得:点电荷q 在斜面上任一点处受力情况,β,则两电荷间的距离r 可用:(2))其中,β在20度到150度之间。
,所以加速度的变化主要随的变化而变化。
2.2程序及结果为了绘图方便,这里设,β在20度到150度之间变化,A 对角度变化的曲线如图4所示,能很直观的看出,加速度先增加再减小。
图4 质点加速度曲线3.MATLAB 在电学实验中的应用MATLAB/Simulink 仿真实验模拟真实实验环境,让学生在虚拟的实验环境中进行实验,能便捷地熟悉物理实验并学习物理知识。
下面通过MATLAB/Simulink 建模仿真对伏安法测电阻实验进行分析。
3.1 建模根据欧姆定律,在同一电路中,通过某一导体R 的电流I U 成正比,即:4)(4)进行变形,可以导出测量电阻的公式为:5)建模用MATLAB/Simulink 中的Sim-powerSystems 的电力系统模型库,选取相应的电阻、直流电源、电压表、电流表及示波器等基本功能模块,模块按照图5建立仿真模型,设置电压值为24V ,电阻为48欧姆,然后进行仿真。
matlab计算物理
matlab计算物理【原创版】目录1.MATLAB 简介2.MATLAB 在物理学中的应用3.MATLAB 的优点4.如何学习和使用 MATLAB5.总结正文一、MATLAB 简介MATLAB 是一种广泛应用于科学计算和工程设计的编程语言。
它被广泛应用于各个领域,包括物理学、数学、工程学等。
MATLAB 的最大特点是它的矩阵计算能力,使得它在处理大量的数据时具有非常高的效率。
二、MATLAB 在物理学中的应用在物理学中,MATLAB 可以用于许多物理现象的模拟和计算。
例如,它可以用于模拟电磁场、计算流体力学、模拟量子物理现象等。
通过MATLAB,物理学家可以更方便地进行数据分析和模型构建。
三、MATLAB 的优点MATLAB 有以下几个优点:1.高效的矩阵计算能力:MATLAB 的矩阵计算能力非常强大,可以处理大量的数据。
2.丰富的函数库:MATLAB 有丰富的函数库,包括数学、统计、图像处理等多个领域。
3.易于学习和使用:MATLAB 的语法简单,易于上手。
四、如何学习和使用 MATLAB如果你想学习和使用 MATLAB,可以按照以下步骤进行:1.下载并安装 MATLAB:你可以在 MATLAB 的官方网站上下载并安装MATLAB。
2.学习 MATLAB 的基本语法:你可以通过阅读 MATLAB 的官方文档或参加 MATLAB 的培训课程来学习 MATLAB 的基本语法。
3.练习 MATLAB 编程:你可以通过编写一些简单的程序来熟悉MATLAB 的编程环境。
4.应用 MATLAB 进行物理计算:你可以使用 MATLAB 的函数库进行物理计算,例如使用 MATLAB 进行数据分析、模拟物理现象等。
五、总结总的来说,MATLAB 是一种非常强大的科学计算工具,它在物理学中的应用非常广泛。
适合用matlab解决的经典物理例题
适合用matlab解决的经典物理例题(最新版)目录1.引言:介绍 MATLAB 以及经典物理例题2.MATLAB 在物理学中的应用3.经典物理例题解析4.总结:MATLAB 在解决经典物理例题中的优势正文一、引言MATLAB 是一种广泛应用于科学计算和工程设计的语言,特别是在物理学领域,它为学生们提供了解决复杂数学问题的强大工具。
经典物理例题是物理学教学中的重要组成部分,通过解决这些例题,学生可以深入理解物理学的基本原理和定律。
本文将介绍如何利用 MATLAB 解决一些经典的物理例题。
二、MATLAB 在物理学中的应用MATLAB 在物理学中有着广泛的应用,例如在数据分析、图像处理、数值计算等方面。
它提供了丰富的函数库,如数学运算、矩阵操作、绘图功能等,使得物理学家和学生可以更加高效地完成各种计算任务。
此外,MATLAB 还具有高级的编程功能,可以实现自定义的算法和模型,为解决复杂的物理问题提供了可能。
三、经典物理例题解析1.牛顿运动定律利用 MATLAB 可以方便地解决与牛顿运动定律相关的问题,如物体在给定力的作用下的运动轨迹、速度和加速度等。
通过绘制物体运动的速度- 时间图像,可以直观地观察到物体在不同时间内的运动状态。
2.库仑定律MATLAB 可以用于解决库仑定律相关的问题,例如两个点电荷之间的电力大小和方向、电场强度分布等。
通过 MATLAB 的数值计算功能,可以得到精确的解。
3.波动方程MATLAB 可以方便地解决波动方程相关的问题,如声波、光波等。
通过绘制波动图像,可以直观地观察到波的传播过程和特性。
四、总结MATLAB 作为一种功能强大的科学计算工具,在解决经典物理例题中具有明显的优势。
它不仅可以完成复杂的数学运算,还可以进行数据分析和绘图,使得物理问题的解决更加直观和形象。
同时,MATLAB 的编程功能为解决复杂的物理问题提供了可能。
matlab在力学中的应用
Matlab在力学中有广泛的应用。
以下是一些示例:1.力学模拟和分析:Matlab提供了丰富的工具和函数,可以进行力学系统的建模、仿真和分析。
可以使用Matlab进行刚体力学、弹性力学、流体力学等各种力学问题的数值模拟和分析。
2.力学数据处理和可视化:在力学实验中,可以使用Matlab对实验数据进行处理和分析。
Matlab提供了强大的数据处理功能和绘图工具,可以用来处理和可视化力学实验数据,生成曲线、图像和动画等结果。
3.力学方程求解:Matlab具备解常微分方程和偏微分方程的能力,可以用来求解力学问题中的方程和模型。
可以使用Matlab解析解方程,数值求解微分方程,并应用于刚体力学、振动、弹性力学等领域。
4.优化和最优化问题:力学中经常需要进行优化和最优化问题的求解,例如寻找最佳设计、最优控制、力学系统的最佳参数等。
Matlab提供了强大的优化工具和算法,可以用来解决这些问题。
5.控制系统设计和分析:在力学中,控制系统的设计和分析是一个重要的方面。
Matlab 提供了广泛的控制系统工具箱,可以用于设计和分析力学系统的控制器,进行稳定性分析、频谱分析、响应分析等。
6.有限元分析:有限元分析是力学中一种常用的数值分析方法。
Matlab提供了用于有限元分析的工具和函数,可以进行结构强度分析、振动分析、流体力学分析等。
可以使用Matlab进行有限元模型的建立、求解和后处理。
总之,Matlab在力学中具有广泛的应用,可以用于力学建模、仿真、分析,数据处理和可视化,方程求解,优化问题,控制系统设计和分析,以及有限元分析等领域。
它为力学工程师和研究人员提供了强大的工具和资源,帮助他们解决力学问题并进行深入的研究和开发。
matlab在中学物理教学中的基本应用
matlab在中学物理教学中的基本应用文章标题:深入探讨MATLAB在中学物理教学中的基本应用随着科技的发展,计算机在教学中扮演着越来越重要的角色。
MATLAB作为一种强大的计算工具,在中学物理教学中也有着广泛的应用。
本文将从多个角度探讨MATLAB在中学物理教学中的基本应用,以便让读者更深入地理解这一话题。
一、MATLAB在物理实验中的应用在中学物理教学中,实验是非常重要的一环。
通过实验,学生可以亲身感受物理原理,加深对知识点的理解。
而MATLAB在物理实验中的应用,可以让实验过程更加直观、高效。
利用MATLAB可以更精确地绘制实验数据的曲线图,进行数据的拟合和分析。
还可以利用MATLAB对实验数据进行快速处理,帮助学生更好地理解实验结果。
二、MATLAB在物理模拟中的应用除了实验,物理模拟也是物理教学中常用的手段。
通过MATLAB,可以轻松实现物理现象的模拟,如自由落体运动、简谐振动等。
学生可以通过自行编写MATLAB代码,模拟各种物理现象,从而深入理解物理原理。
通过调整参数、改变条件,他们可以在模拟中观察到不同的现象,从而更好地理解物理规律。
三、MATLAB在物理问题的数学建模中的应用物理问题与数学紧密相连,而MATLAB恰恰是一个强大的数学工具。
在中学物理教学中,通过MATLAB可以进行物理问题的数学建模,从而将抽象的物理问题转化为具体的数学模型。
可以利用MATLAB进行势能、动能方程的求解,复杂物理系统的模拟等。
这一应用不仅培养了学生的数学建模能力,同时也加深了他们对物理问题的理解。
总结回顾MATLAB作为一种强大的计算工具,为中学物理教学提供了丰富的应用场景。
通过MATLAB,可以帮助学生更好地理解物理原理,提高他们的实验能力和数学建模能力。
在实验、模拟和数学建模中,MATLAB都发挥着不可替代的作用。
我们应该充分发挥MATLAB在中学物理教学中的作用,从而提高教学质量,培养更多对物理感兴趣的学生。
Matlab在大学物理中应用的几个实例
●
●
●
, . ● ● ,
●
●
,
【
,
一
。
一
等 去 瓦 l[ 丽 √
0 函
4M
¨
q
+
缸
x
:
R
。
]
I
,
、
果
,
利 用 函 数 命令 d i f f ( ) 很 容 易计 算 出结 其语 句 如 下 :
> > s y rn s > >
、
,
s
q q
s /
X r
R
‘ ‘
y
"
m s
X
Y
‘
一
,
。
、
、
。
’ 】 1 M a t l a b 7 0 有关功能简 介【
一
,
任 意 条 电场 线 应满 足 方程 :
。
,
,
等 型f
:
出
E
.
r
.
…
v Y
l
、 , (1 )
‘
,
。
一
,
毡 求 解式 ( 1 ) 可 得 【
、
”
,
、
圆环 轴 线 上
(
£
r
)
一
点 的 电势为0 ∥
。
‘
砥
焉亍意亍
垒! !! !
,+
一
,
!
x
。
一
,
(j √
“
,
+ x
。
√
+
y
。
。
MATLAB在大学物理中的应用共3篇
MATLAB在大学物理中的应用共3篇MATLAB在大学物理中的应用1MATLAB在大学物理中的应用MATLAB是一种数学软件,被广泛应用于大学物理的教学和研究中。
其功能强大,包含了许多求解数学和物理问题所需的工具和函数。
本文将探讨MATLAB在大学物理中的应用。
一、矢量和矩阵计算MATLAB中的矢量和矩阵计算功能可以方便地帮助学生学习和理解物理中的向量和矩阵。
例如,通过MATLAB可以进行向量叉乘、点乘等运算,帮助学生更深入了解向量的性质和运算规律。
在矩阵方面,MATLAB可以进行矩阵的转置、逆矩阵的计算、特征值和特征向量的计算等操作,这些在物理中常常遇到的矩阵计算可以大大简化学生的计算过程。
二、数值计算和绘图在物理中,我们经常需要对一些物理现象进行数值计算和绘图。
MATLAB中的数值计算和绘图功能可以方便地进行这些操作。
例如,使用MATLAB可以进行微积分的数值计算,包括求导、积分等。
同时,MATLAB还可以绘制函数图像、动画、示波器等,帮助学生更加直观地理解物理现象。
三、符号计算在大学物理中,符号计算也是重要的一部分。
MATLAB可以进行符号计算,包括求解方程、求解微分方程、求导、积分等。
这些功能可以帮助学生更加深入地理解物理中的数学公式和方程,同时也方便了学生在计算中的操作。
四、数值模拟MATLAB还可以进行数值模拟,模拟物理问题的数值计算和分析。
例如,可以使用MATLAB模拟机械振动、光学成像等。
数值模拟可以帮助学生更好地理解物理中的现象、规律和数学模型,同时也可以提高学生的实际操作能力。
五、数据分析最后,在大学物理中,数据分析也是一个重要的环节。
MATLAB 中可以进行数据分析,包括数据的导入、处理、分析和可视化等。
数据分析可以帮助学生更加准确地分析物理中的数据,进一步深入理解物理现象。
综上所述,MATLAB在大学物理中的应用非常广泛,涉及到矢量和矩阵计算、数值计算和绘图、符号计算、数值模拟和数据分析等多个方面。
MATLAB在大学物理课程及实验教学中的应用
J R L oF QI GHA AT O L T oU NA N IN I NA I ⅢS U I E S T NV R IY
( d ct nSineE io E uai ce c dt n) o i
2l 0 1年 第 5
关键词 : 大学物理 ; A L B 可视化 ; M TA ; 数据处理
M ATL AB —b sd Vi aiai n a d Da a P o esn o l g h sc a h n — a e s l t n t r c si g f rCol e P y isTe c i g u z o e LIZh n DONG o g , og, S n CHUIW a n—n a S in , HEN u —p n 4 W e g
易 懂的特点. 应用于物理课程教学中可使学生头脑中形成生动形象 的动态物理图像. 化抽象为直观, 帮助
.
理解大学物理的理论. 以下举例说明物理现象 、 规律的可视化 : 动表述的是悬浮在流体 中的微粒受到流体分子与粒子 的碰撞而发生 的不停息 的随机运动 , 其运 『 剧烈程度受温度影响. 这一运动比较抽象 , 我们可以用 M T A A L B仿真 , 将其可视化. 模拟效果如图 1 所
( .. . . col f hs s Eet ncIfr t nE g er g Qnhi ao at sU i r t, iigQnh i 107 C ia 12 3 4 Sh o o P yi & l r i no i ni ei , iga N t nli nv syX nn iga 80 0 , hn ) c c o ma o n n i ie ei
’咀 子 二^ 上l J u W " 1 子 J 7 鬯 ’ Y圳 ’ A 1 3 d T , 疋 同 千 上 , 付 旺 ^¨ 匕 ,丁 J 1 r 、 日匕/J 。 , 口 J
matlab编程例题自由落体运动
自由落体运动是物理学中的一个基础概念,它描述了在没有空气阻力作用下物体在重力场中运动的规律。
在实际工程和科学研究中,我们常常需要使用matlab编程来模拟自由落体运动,并分析物体的轨迹、速度、加速度等参数。
本文将通过实际例题,介绍如何利用matlab编程来模拟自由落体运动,并分析其运动规律。
1. 问题描述假设有一个质量为m的物体从高度h自由落体,重力加速度为g,求在不考虑空气阻力的情况下,物体落地时的速度、运动时间以及运动过程中的轨迹。
2. 解决思路为了解决上述问题,我们可以采用matlab编程,利用欧拉法对自由落体运动进行模拟。
欧拉法是一种常用的数值求解算法,其基本思想是根据微分方程的定义,通过离散化时间和空间,逐步迭代计算物体运动过程中的位置、速度等参数。
3. matlab代码实现下面是利用matlab编写的自由落体运动模拟程序代码:```matlabclearclch = input('请输入初始高度h:'); 初始高度g = 9.8; 重力加速度t = sqrt(2*h/g); 运动时间v = g*t; 落地时的速度disp(['物体落地时的速度为:', num2str(v)]);disp(['物体运动时间为:', num2str(t)]);模拟自由落体运动轨迹t = 0:0.01:t; 时间间隔y = h - 0.5*g*t.^2; 物体高度plot(t, y);xlabel('时间');ylabel('高度');title('自由落体运动轨迹');```4. 代码说明上述matlab代码首先输入初始高度h,并利用公式计算出运动时间t 和落地时的速度v,然后利用plot函数绘制自由落体运动的轨迹图。
在模拟运行程序后,我们可以得到落地时的速度、运动时间以及运动轨迹,并通过图形直观地观察物体的运动规律。
利用matlab进行仿真的案例
利用matlab进行仿真的案例利用Matlab进行仿真可以涉及多个领域的案例,下面列举10个案例:1. 汽车碰撞仿真:利用Matlab中的物理仿真库,可以模拟汽车碰撞的过程,分析碰撞时车辆的变形、撞击力等参数。
可以根据不同的碰撞角度和速度,评估不同碰撞条件下的安全性能。
2. 电力系统仿真:利用Matlab中的电力系统仿真工具,可以模拟电力系统的运行情况,包括电压、电流、功率等参数的变化。
可以用于分析电力系统的稳定性、短路故障等问题,并进行相应的优化设计。
3. 通信系统仿真:利用Matlab中的通信系统仿真工具箱,可以模拟无线通信系统的传输过程,包括信号的发送、接收、调制解调等环节。
可以用于评估不同调制方式、编码方式等对通信系统性能的影响。
4. 智能控制仿真:利用Matlab中的控制系统仿真工具,可以模拟各种控制系统的运行情况,包括PID控制、模糊控制、神经网络控制等。
可以用于设计、优化和评估各种控制算法的性能。
5. 雷达系统仿真:利用Matlab中的雷达仿真工具,可以模拟雷达系统的工作原理和性能,包括发射、接收、信号处理等过程。
可以用于评估雷达系统的探测能力、跟踪精度等指标,并进行系统参数的优化设计。
6. 气候变化模拟:利用Matlab中的气候模型,可以模拟气候系统的变化过程,包括温度、降水、风速等参数的变化。
可以用于研究气候变化对生态环境、农业生产等方面的影响,以及制定相应的应对策略。
7. 人体生理仿真:利用Matlab中的生理仿真工具箱,可以模拟人体的生理过程,包括心血管系统、呼吸系统、神经系统等。
可以用于研究不同疾病、药物对人体的影响,以及评估各种治疗方案的效果。
8. 金融市场仿真:利用Matlab中的金融工具箱,可以模拟金融市场的价格变化过程,包括股票、期货、汇率等。
可以用于研究不同投资策略、风险管理方法等对投资收益的影响,并进行相应的决策分析。
9. 电子器件仿真:利用Matlab中的电子器件仿真工具,可以模拟各种电子器件的工作原理和性能,包括二极管、晶体管、集成电路等。
巧用“MATLAB”解决物理中的疑难问题
作者: 徐富清
作者机构: 江苏省扬州市邗江区公道中学
出版物刊名: 实验教学与仪器
页码: 59-60页
主题词: 疑难问题;计算软件;可视化;最小速度;发射角;物理问题;计算结果;发射速度;分子势能;作用力
摘要: 'MATLAB'是一款非常强大的科学计算软件,在众多优点中比较突出的是以可视化的方式将计算结果快速、准确地用图形的方式形象、直观地描绘出来。
它可以帮助学生解决因数学困难而无法解决的物理问题。
下面笔者举例说明'MATLAB'在解决物理疑难问题中的应用。
例1在我方前沿阵地2000m远有一座高70m。
- 1、下载文档前请自行甄别文档内容的完整性,平台不提供额外的编辑、内容补充、找答案等附加服务。
- 2、"仅部分预览"的文档,不可在线预览部分如存在完整性等问题,可反馈申请退款(可完整预览的文档不适用该条件!)。
- 3、如文档侵犯您的权益,请联系客服反馈,我们会尽快为您处理(人工客服工作时间:9:00-18:30)。
1. 单列波%%单列波t=0:0.001:10;A=input('振幅A=');w=input('频率w=');a=input('a=');y=A.*sin(w.*t+a);plot(t,y);pause(1),sound(y);ylabel('y'),xlabel('t')2. %%光的单缝衍射现象Lambda=500e-9; %a=input('a='); % 可取0.2e-3,1e-3,2e-3三种情况z=1 %ymax=3*Lambda*z/a; %Ny=51; %ys=linspace(-ymax,ymax,Ny); %NPoints=51; %yPoint=linspace(-a/2,a/2,NPoints); %for j=1:Ny %L=sqrt((ys(j)-yPoint).^2+z^2); %Phi=2*pi.*(L-z)./Lambda; %SumCos=sum(cos(Phi)); %SumSin=sum(sin(Phi)); %B(j)=(SumCos^2+SumSin^2)/NPoints^2; %endclf,plot(ys,B,'*',ys,B);grid; %3. %%用毕奥-沙伐尔定律计算电流环产生的磁场mu0=4*pi*1e-7;I0=5.0;Rh=1;C0=mu0/(4*pi)*I0;NGx=21;NGy=21;x=linspace(-Rh,Rh,NGx);y=linspace(-3,3,20);y=x;Nh=20;theta0=linspace(0,2*pi,Nh+1);theta1=theta0(1:Nh);y1=Rh*cos(theta1);z1=Rh*sin(theta1);theta2=theta0(2:Nh+1);y2=Rh*cos(theta2);z2=Rh*sin(theta2);xc=0;yc=(y2+y1)/2;zc=(z2+z1)/2;for i=1:NGyfor j=1:NGxrx=x(j)-xc;ry=y(i)-yc;rz=0-zc;r3=sqrt(rx.^2+ry.^2+rz.^2).^3;dlXr_x=dly.*rz-dlz.*ry;dlXr_y=dly.*rx-dlx.*rz;Bx(i,j)=sum(C0*dlXr_x./r3);By(i,j)=sum(C0*dlXr_y./r3); endendclf;quiver(x,y,Bx,By);4. %%多普勒效应x0=500;v=50;y0=20;c=330;w=1000;t=0:0.001:30;r=sqrt((x0-v*t).^2+y0.^2);t1=t-r/c;u=sin(w*t)+sin(1.1*w*t);u1=sin(w*t1)+sin(1.1*w*t1); sound(u);pause(5);sound(u1);5.亥姆霍兹线圈clear allmu0=4*pi*1e-7;I0=5.0;Rh=1;C0=mu0/(4*pi)*I0;NGx=21;NGy=21;x=linspace(-Rh,Rh,NGx);y=linspace(-Rh,Rh,NGy);Nh=20;theta0=linspace(0,2*pi,Nh+1); theta1=theta0(1:Nh);y1=Rh*cos(theta1);z1=Rh*sin(theta1);theta2=theta0(2:Nh+1);y2=Rh*cos(theta2);z2=Rh*sin(theta2);xc=0;yc=(y2+y1)/2;zc=(z2+z1)/2;for i=1:NGyfor j=1:NGxrx=x(j)-xc;ry=y(i)-yc;rz=0-zc;r3=sqrt(rx.^2+ry.^2+rz.^2).^3;dlXr_x=dly.*rz-dlz.*ry;dlXr_y=dly.*rx-dlx.*rz;Bx(i,j)=sum(C0*dlXr_x./r3);By(i,j)=sum(C0*dlXr_y./r3);endendBax=Bx(:,11:21)+Bx(:,1:11);Bay=By(:,11:21)+By(:,1:11);subplot(1,2,1)mesh(x(11:21),y,Bax);xlabel('x');ylabel('y'); subplot(1,2,2),plot(y,Bax),grid,xlabel('y');ylabel('Bx'); 6.库仑引力clear all;N=input('电荷数目N:');for ic=1:N %fprintf('-----\n 对电荷#%g\n',ic);rc=input('电荷位置[x y](米):');x(ic)=rc(1); %y(ic)=rc(2); %q(ic)=input('输入电荷量(库仑)'); endE0=8.85e-12; %C0=1/(4*pi*E0); %for ic=1:NFx=0.0; Fy=0:0; %for jc=1:N %if(ic~=jc) %xij=x(ic)-x(jc); yij=y(ic)-y(jc);Rij=sqrt(xij^2+yij^2); %Fx=Fx+C0*q(ic)*xij/Rij^3;Fy=Fy+C0*q(ic)*yij/Rij^3;endendfprintf('其它电荷作用在电荷#%g上的合力为:\n',ic); %fprintf('x-分量:%gN\n',Fx);fprintf('y-分量:%gN\n',Fy)end7.李萨如图形% lisaru.msyms t a1 a2 w1 w2x=cos(w1.*t+a1);y=sin(w2.*t+a2);a1=input('a1=');a2=input('a2=');w1=input('w1=');w2=input('w2=');tf=10;Ns=1000;t=linspace(0,tf,Ns);dt=tf/(Ns-1); %分Ns个点,求出时间增量dtxplot=eval(x);yplot=eval(y); %计算Ns个点的位置x(t),y(t)figure(gcf);subplot(1,2,1),for i=1:750plot(yplot(1:i),xplot(1:i)); %画点的轨迹图axis('equal'); grid ; %使两轴比例相同pause(0.01)end8.耦合振子m1=2;m2=2;K1=16;K2=4; % x0=[1;0];xd0=[6;6];tf=10; %M=[m1,0;0,m2];K=[K1+K2,-K2;-K2,K1+K2] ; %u(:,s)�[u,L]=eig(K,M) ; %t=linspace(0,tf,101);x=zeros(2,101); % for s=1:2 %alfa=sqrt(u(:,s)'*M*u(:,s)) ; %u(:,s)= u(:,s)/alfa; %w(j)=sqrt(L(j,j)); % xt=u(:,j)*(u(:,j)'*M*x0*cos(w(j)*t)+u(:,j)'*M*xd0/w(j)*sin(w(j)*t));x=x+xt; %endfor r=1:2 %x1,x1?? subplot(2,1,r)plot(t.x(r,:)),grid;xlabel('xxx');ylabel(['yyy',num2str(r)]);end9.拍频%%t=0:0.001:10;a1=input('??1=');w1=input('??1=');a2=input('??2=');w2=input('??2=');y1=a1*sin(w1*t);y2=a2*sin(w2*t);y=y1+y2;subplot(3,1,1),plot(t,y1),ylabel('y1')subplot(3,1,2),plot(t,y2),ylabel('y2')subplot(3,1,3),plot(t,y),ylabel('y'),xlabel('t')pause,sound(y1);pause(5),sound(y2);pause(5),sound(y),pausesubplot(1,1,1)10。
物块下滑m1=input('m1=');m2=input('m2=');theta=input('theta(度)=');theta=theta*pi/180; g=9.81;A=[m1*cos(theta),-m1,-sin(theta),0;...m1*sin(theta),0,cos(theta),0;...0,m2, -sin(theta),0;...0,0,-cos(theta),1];B=[0,m1*g,0,m2*g]'; X=A\B;a1=X(1), a2=X(2), N1=X(3),N2=X(4)11。
循环%%xunhuanR=8.31;gama=1.4; %(注:由于在MATLAB中无拉丁文,所以用gama代替?)nMoles=0.5;P(1)=1e5;V(1)=0.012;WTotal=0;QTotal=0;iPoint=1;NCurve=100;PPlot=P(1);VPlot=V(1);%变量和图形初始化,输入气体的摩尔数nMoles,初始压力P(1),初始容积V(1),%气体常数R=8.314;给定起始总功WTotal=0;点序号iPoint=1;画等温线用%的点数NCurve=100;P-V图第一点坐标PPlot=P(1);VPlot=V(1)T(1)=P(1)*V(1)/(nMoles*R); %算出初始温度%为了进入循环,先要设两个不相等的PathType和QuitType值QuitType=5;PathType=0;while(PathType~=QuitType) %在菜单上选择‘退出’之前不断循环,%选择路径类型或退出iPoint=iPoint+1; %下一点fprintf('对过程#%g\n',iPoint-1);PathType=menu(sprintf('过程%g:选择下一路径',iPoint-1),...'等压','等容','等温','绝热','退出'); %图形界面菜单生成语句switch PathTypecase 1 %等压路径V(iPoint)=input('输入新容积:');P(iPoint)=P(iPoint-1); %压力不变T(iPoint)=P(iPoint)*V(iPoint)/(nMoles*R); %按新容积算出温度W=P(iPoint)*(V(iPoint)-V(iPoint-1)); %计算等压过程所做的功Q=(gama*nMoles*R/(gama-1))*(T(iPoint)-T(iPoint-1));VPlot=[VPlot,V(iPoint)]; %加上新的容积和压力点,用以绘图PPlot=[PPlot,P(iPoint)];case 2 %等容路径P(iPoint)=input('输入新压力:');V(iPoint)=V(iPoint-1); %容积不变T(iPoint)=P(iPoint)*V(iPoint)/(nMoles*R); %按新压力算出温度W=0; %等容路径上所做的功为零Q=(nMoles*R/(gama-1))*(T(iPoint)-T(iPoint-1));VPlot=[VPlot,V(iPoint)]; %加上绘图用的新容积和压力点PPlot=[PPlot,P(iPoint)];case 3 %等温路径V(iPoint)=input('输入新容积:');T(iPoint)=T(iPoint-1); %温度不变P(iPoint)=nMoles*R*T(iPoint)/V(iPoint); %按新容积求新压力W=nMoles*R*T(iPoint)*log(V(iPoint)/V(iPoint-1)); %求所做的功Q=W;%用元素群运算求等温路径上的P和V,加进绘图数据中VNew=linspace(V(iPoint-1),V(iPoint),NCurve);PNew=nMoles*R*T(iPoint)./VNew;VPlot=[VPlot,VNew]; %将新的V,P点加入绘图数据中PPlot=[PPlot,PNew];case 4 %绝热路径V(iPoint)=input('输入新容积:');P(iPoint)=P(iPoint-1)*V(iPoint-1)^gama/V(iPoint)^gama; %按新容积求新压力T(iPoint)=T(iPoint-1)*V(iPoint-1)^(gama-1)/V(iPoint)^(gama-1); %按新容积求新温度W=-(nMoles*R/(gama-1))*(T(iPoint)-T(iPoint-1)); %求所做的功Q=0;VNew=linspace(V(iPoint-1),V(iPoint),NCurve);PNew=P(iPoint-1)*V(iPoint-1)^gama./VNew.^gama;VPlot=[VPlot,VNew]; %将新的V,P点加入绘图数据中PPlot=[PPlot,PNew];otherwiseend%画出到目前为止的PV图if(PathType~=QuitType)WTotal=WTotal+W; %将新做的功加进总功if(Q>0)QTotal=QTotal+Q;endfigure(gcf);plot(V,P,'o',VPlot,PPlot,'-')%图形窗移前,绘图size=axis;axis([size(1)*0.9,size(2)*1.1,size(3)*0.9,size(4)*1.1]);%上两句用于美观图象,使曲线不要紧靠边框%标注语句略endendWTotalQTotaly=WTotal/QTotal12。