数学模型试验

本文格式为Word版，下载可任意编辑——数学模型试验

数学模型课程试验报告

课程名称姓名试验名称试验目的及要求试验环境班级学号初等模型一、上机时请按试验要求进行试验。二、每次上机请完成试验报告。操作系统：处理器：内存：系统类型：MATLAB软件版本：R2023b（8.6.0.267246）64-bit(win64)试验日期试验成绩1

实验内容与步骤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比较）。2

实验程序源代码1.1%%clear;clc;formatcompact;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;formatcompact;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=[1248];t=[7.216.886.325.84];logt=log(t);logn=log(n);p=polyfit(logn,logt,1);beta=p(1)alfa=exp(p(2))3

2.2%%n=[1248]';t=[7.216.886.325.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.12115.5684.9666.6068.0471.64,...82.08132.84185.04226.80246.60250.92,...261.00271.44273.96279.00291.60302.04,...310.68290.52281.16248.40210.24186.84];s=0;fori=1:24s=s+f(i);endg=s/24t2=0:0.01:23;plot(t,f,t2,g,'r-');text(10,213.67,'g=203.67');gridon

4

3.1.2%%t=0:23;f=[150.12115.5684.9666.6068.0471.64,...82.08132.84185.04226.80246.60250.92,...261.00271.44273.96279.00291.60302.04,...310.68290.52281.16248.40210.24186.84];s=0;fori=1:24s=s+f(i);endg=s/24;c1(1)=0;forj=1:23c1(j+1)=c1(j)+f(j)-g;endc1M=min(c1)c2(1)=-M;fork=1:23c2(k+1)=c2(k)+f(k)-g;endc2plot(t,c1,t,c2);text(7,160,'c(0)=876.15');text(7,-760,'c(0)=0');xlabel('t');ylabel('c');gridon3.2.1%%functiony=f(t)tt=0:23;ft=[150.12115.5684.9666.6068.0471.64,...8