数学模型实验-.docVIP

数学模型课程实验报告 课程名称 班级 实验日期 姓名 学号 实验成绩 实验名称 初等模型 实验目的及要求 一、上机时请按实验要求进行实验。 二、每次上机请完成实验报告。 实验环境 操作系统：处理器： 内存：系统类型： MATLAB 软件版本：R2015b（8.6.0.267246） 64-bit(win64) 实 验 内 容 与 步 骤 1. ①（验证）分别计算出LCLV, CCLV和TCLV三个3行1列的列向量，然后输出结果，并与P26的表2（教材）比较。 ②（编程）对于LCAV, CCAV和TCAV，编写类似①的程序，并运行，结果与P26的表3（教材）比较。 2. ①运行问题(1)中的程序。 ②编程解决问题(2)：实际值与计算值比较（数据比较和和拟合图形）。 3.1 ① 求g，画f(t)和g的图形（与P35图1比较）。 ② 求c(t), t=0, 1, 2, …, 23, c(0)=0，并求其中的最小值M（与P36表3比较）。 求c(t), t=0, 1, 2, …, 23, c(0)=-M（与P36表4比较）。画c(t)分别当c(0)和c(-M)时的图形（与P37图2比较）。 3.2 (1) 求g，画f(t)和g的图形（与P35图1比较）。 (2) 求c(t), 0≤t≤23, c(0)=0时的最小值M。画c(t)初值条件分别为c(0)=0和c(0)=-M时的图形（与P37图2比较）。 实 验 程 序 源 代 码 1.1%% clear;clc;format compact; R1=22.5; R2=58; d=10^(-3)*[1.6,0.74,0.32]; rho=[121,387,800]; LCLV=pi*(R2^2-R1^2)./d; CCLV=rho.*LCLV/10^6; TCLV=CCLV/(0.62*60); s= ;S=[s;s;s]; [num2str(round(LCLV)),S,... num2str(round(CCLV)),S,... num2str(round(TCLV))] 1.2 %% clear;clc;format compact; R1=22.5; R2=58; d=10^(-3)*[1.6,0.74,0.32]; rho=[121,387,800]; LCAV=pi*R2^2./(2*d); CCAV=rho.*LCAV/10^6; TCAV=CCAV/(0.62*60); s= ;S=[s;s;s]; [num2str(round(LCAV)),S,... num2str(round(CCAV)),S,... num2str(round(TCAV))] 2.1 %% clear; clc; n=[1 2 4 8]; t=[ 7.21 6.88 6.32 5.84]; logt=log(t); logn=log(n); p=polyfit(logn,logt,1); beta=p(1) alfa=exp(p(2)) 2.2 %% n=[1 2 4 8]; t=[7.21 6.88 6.32 5.84]; logt=log(t); logn=log(n); p=polyfit(logn,logt,1); beta=p(1); alfa=exp(p(2)); t2=alfa*n.^beta; [n,t,t2] a=0:0.01:10; t3=alfa*a.^beta; plot(n,t,x,a,t3); 3.1.1 %% t=0:23; f=[150.12 115.56 84.96 66.60 68.04 71.64,... 82.08 132.84 185.04 226.80 246.60 250.92,... 261.00 271.44 273.96 279.00 291.60 302.04,... 310.68 290.52 281.16 248.40 210.24 186.84]; s=0; for i=1:24 s=s+f(i); end g=s/24 t2=0:0.01:23; plot(t,f,t2,g,r-); text(10,213.67,g=203.67); grid on 3.1.2 %% t=0:23; f=[150.12 115.56 84.96 66.60 68.04 71.64,... 82.08 132.84 185.04 226.80 246.60 250.92,... 261.00 271.44 273.96 279.00 291.60 302.04,... 310.68 290.52 281.16 248.4