第6章在普通物理中的应用.pptVIP

第6章 在普通物理中的应用 【例6-1-1】温度单位转换 命题：写出一个程序，能把用户输入的摄氏温度转为华氏，也可反求。 解：◆建模 两种温度之间的转换公式为： 摄氏变华氏 华氏变摄氏 程序中要先考虑由用户选择转换的方向，再给数据。 程序exn611 k = input(选择1：摄氏变华氏；选择2：华氏变摄氏；… 键入数字1或2：); Tin = input(输入待变换的温度(允许输入数组):) ; if k==1 Tout = Tin*9/5 +32; % 摄氏转华氏 elseif k==2 Tout = (Tin-32)*5/9; % 华氏转摄氏 else disp(未给转换方向,转换无效), end s = [华氏;摄氏]; s1 = [转换后的温度为,s(k,:),num2str(Tout),度], % 注意此语句的编写方法 【例6-1-2】多种单位间的换算 写出一个程序，能把用户输入的长度单位在厘米、米、千米、英寸、英尺、英里、市尺、市里之间任意转换。 解：◆建模 这里采取的技巧是分成两步，先把输入量变换为米，第二步再把米变换为输出单位，另外，把变换常数直接表示为一个数组，选择单位的序号也就成了数组的下标；这样程序就比较简明易读。程序如下： 长度单位换算程序ex612.m clear all; disp( 长度单位换算程序) fprintf(长度单位: \n); % 选择输入输出的单位 fprintf( 1) 厘米 2) 米 3) 千米 4) 英寸 \n); fprintf( 5) 英尺 6) 英哩 7) 市尺 8) 市里 \n); InUnits = input(选择输入单位编号: ); OutUnits = input(选择输出单位编号: ); % 令各种单位对米的变换常数数组为ToMeter ToMeter = [0.01, 1.00, 1000.0, 0.0254, 0.3048, 1609.3, 1/3, 500 ]; 程序ex612.m(续) FrmMeter= 1./ ToMeter; Value = input(输入待变换的值(0为退出): ); while( Value ~= 0 ) ValueinM = Value*ToMeter(InUnits); % 把输入值变为米 NewValue = ValueinM*FrmMeter(OutUnits); % 把米变为输出单位 fprintf(变换后的值是 %g \n,NewValue); % 打印变换后的值 Value = input(输入待变换的值(0为退出): ); % 提问下个输入值 end 【例6-1-3】实验数据拟合 命题：设在某一实验中，给某元件加[1,2,3,4,5]v电压，测得的电流为[ 0.2339, 0.3812, 0.5759, 0.8153, 0.9742 ]ma。求此元件的电阻。 解：◆ 建模 模型:设直线的方程为 y= a(1)x +a(2), 待定的系数是a(1),a(2)。将上述数据分别代入x,y, a(1)+a(2)=0.2339 2a(1)+a(2)=0.3812 3a(1)+a(2)=0.5759 4a(1)+a(2)=0.8153 5a(1)+a(2)=0.9742 把这五个方程联立,用矩阵表述, 得 datax *a(1) + ones(N,1) * a(2) = datay 程序exn613 这是一个超定方程组，写成A*a = B，其最小二乘解可以用左除运算符a = A \ B来求得。因此程序如下： lear, datax = [1:5]; datay=[ 0.2339, 0.3812, 0.5759, 0.8153, 0.9742] A = [datax , ones(5,1)]; B = datay; a = A\B, r=1/a(1) plot(datax,datay,o),hold on 运行结果为： a(1) = 0.1905 a(2) = 0.0247 画出曲线如右图。 6.2 力学基础 【例6-2-1】目标相对于射点的高度为 yf,给定初速和射角,计算物体在真空中飞行的时间和距离。 ◆建模：这里目标和射点不在同一高度上，不好求封闭形式的解，用MATLAB使整个计算和绘图过程自动化。其好处是可快速地计算其在不同初速和射角下的飞行时间和距离。关键在求落点时间tf时需要解一个二次线性代数方程 由 解出t，她就是落点时间tf。它会有两个解，我们只取其中一个有效解。再求 程序exn621 clear; y0 = 0; x0 = 0; %