数学实验第三版（主编：李继成 赵小艳）课后练习答案（十四）（1）

实验十四：水塔水流量估计模型

练习一

1.海水温度随着深度的变化而变化,海面温度较高,随着深度的增加,海水温度越来越低.通过验观测得一组海水温度t与深度h的数据如下:

h/m

0

1.5

2.5

4.6

8.2

12.5

16.5

26.5

t/℃

23.5

22.9

20.1

19.1

15.4

11.5

9.5

8.2

要求:

(1)分别用多种数据插值方法找出温度t与深度h之间的近似函数关系;

(2)找出温度变化最快的深度位置,通过查询相关资料,了解这个特殊位置的实际应用价值.

(1)

clc;clear;
format long
h=[0,1.5,2.5,4.6,8.2,12.5,16.5,26.5];
t=[23.5,22.9,20.1,19.1,15.4,11.5,9.5,8.2];
x=0:0.1:26.5;
%拉格朗日插值
y=lglrcz(h,t,x);
figure(1)
plot(h,t,'.',x,y);
x0=ones(length(h),1);xx=[x0];
for i=1:length(t)-1x0=x0.*h';xx=[xx,x0];
end
p=inv(xx)*t'%多项式系数（从低到高排列）
%三次样条插值
yy=interp1(h,t,x,'spline');
figure(2)
plot(x,yy,h,t,'.');
pp=spline(h,t)
function y=lglrcz(x0,y0,x)
n=length(x0);
m=length(x);
for i=1:mz=x(i);s=0.0;for k=1:np=1.0;for j=1:nif j~=kp=p*(z-x0(j))/(x0(k)-x0(j));endends=p*y0(k)+s;endy(i)=s;
end
end

p =

  23.500000000000000

   5.519770260409278

  -6.784280947302284

   2.409485615066580

  -0.388966596330604

   0.030972950891118

  -0.001174710590966

   0.000016740427015

pp =

  包含以下字段的 struct:

      form: 'pp'

    breaks: [0 1.500000000000000 … ]

     coefs: [7×4 double]

    pieces: 7

     order: 4

       dim: 1

>> pp.coefs%三次样条插值的各区间的系数如下

ans =

   0.666350830669392  -3.625403322677568   3.538815615010219  23.500000000000000

   0.666350830669392  -0.626824584665302  -2.839526246004088  22.899999999999999

  -0.286563754712924   1.372227907342874  -2.094122923326514  20.100000000000001

   0.050422214507364  -0.433123747348549  -0.122004187338434  19.100000000000001

  -0.005736766041512   0.111436169330978  -1.280079468201690  15.400000000000000

  -0.000611295331494   0.037431887395469  -0.639946824277967  11.500000000000000

  -0.000611295331494   0.030096343417537  -0.369833901025942   9.500000000000000

（2）通过观察图像得知，变化最快的点应该在[0,8.2]内。

由于拉格朗日插值的龙格现象，我们不妨把x缩小范围，在小范围内求解。

x1=0:0.1:16.5;%由于龙格现象，我们把x缩小范围进行求解
%拉格朗日插值
y=lglrcz(h,t,x1);
figure(1)
plot(h,t,'.',x1,y);
x0=ones(length(h),1);xx=[x0];
for i=1:length(t)-1x0=x0.*h';xx=[xx,x0];
end
p=inv(xx)*t'%多项式系数（从低到高排列）
function y=lglrcz(x0,y0,x)
n=length(x0);
m=length(x);
for i=1:mz=x(i);s=0.0;for k=1:np=1.0;for j=1:nif j~=kp=p*(z-x0(j))/(x0(k)-x0(j));endends=p*y0(k)+s;endy(i)=s;
end
end


pd=polyder(p);
xl=abs(polyval(p,x1));
x1(find(xl==max(xl)))

ans =

   16.5000

感觉龙格现象的影响，拉格朗日插值不太适合此题。

对于三次样条插值，我们在前四个区间内求解。

syms t
f1=0.6664*(t)^3-3.6254*(t)^2 +3.5388*(t) + 23.5;
f2=0.6664*(t-1.5)^3 -0.6268*(t-1.5)^2  -2.8395*(t-1.5) + 22.9;
f3=-0.2866*(t-2.5)^3+1.3722*(t-2.5)^2 -2.0941*(t-2.5)+ 20.1;
f4=0.0504*(t-4.6)^3 -0.4331*(t-4.6)^2  -0.1220*(t-4.6)+ 19.1;
f11=matlabFunction(diff(f1));
f22=matlabFunction(diff(f2));
f33=matlabFunction(diff(f3));
f44=matlabFunction(diff(f4));
t1=0:0.1:1.5;t2=1.5:0.1:2.5;t3=2.5:0.1:4.6;t4=4.6:0.1:8.2;
f111=f11(t0)
f222=f22(t0);
f333=f33(t0);
f444=f44(t0);
max(abs(f111))
max(abs(f222))
max(abs(f333))
max(abs(f444))
%显然最大值在第一个区间内
t1(find(f111==3.5388))

ans =

    3.5388

 ans =

    3.0357

ans =

    2.0941

ans =

    1.3624

ans =

     0

这两种方法得到的结果差别很大，估计是由于选择的插值方法不是很合适，个人疑问：这个数据所给的为何会插值得到奇怪的曲线，刚开始为何深度变大而温度升高呢？

个人建议还是改成拟合较好。

温度变化快的地方可以用来海水温差发电，且发电效率更高。

2.表14.8给出了在低潮时某一平面区域内若干点(x,y)处的水深z值(单位:ft).已知船的吃水深度为5ft.试画出海底的地貌图,并在平面矩形区域(80,196)x(-70,145)内标注哪些地方船要避免进人·

表14.8 水域坐标数据

x

129.0

140.5

103.5

88.0

185.5

195.0

105.5

157.5

107.5

77.0

81.0

162.0

162.0

117.5

y

7.5

141.5

23.0

147.0

22.5

137.5

85.5

-6.5

-81.0

3.0

56.5

-66.5

84.0

-33.5

z

4

8

6

8

6

8

8

9

9

8

8

9

4

9

clc;clear;
x=[129.0,140.5,103.5,88.0,185.5,195.0,105.5,157.5,107.5,77.0,81.0,162.0,162.0,117.5];
y=[7.5,141.5,23.0,147.0,22.5,137.5,85.5,-6.5,-81.0,3.0,56.5,-66.5,84.0,-33.5];
z=[4,8,6,8,6,8,8,9,9,8,8,9,4,9];
xx=77:196;
yy=-81:145;
[xxx,yyy]=meshgrid(xx,yy);
zzz=griddata(x,y,z,xxx,yyy,'cubic');
figure(1)
mesh(xxx,yyy,zzz);
figure(2)
contourf(xxx,yyy,zzz);%等高线
figure(3)
rectangle('Position',[80,-70,116,215]);%画矩形，[左下横，左下纵，长，宽]
for xi=80:196for yi=-70:145if zzz(find(yy==yi),find(xx==xi))>5hold onplot(xi,yi,'r.');endend
end

3.估计煤矿的储量.表14.9给出了某露天煤矿在平面矩形区域(1100mx700m)内,纵横均匀的网格交点处测得的煤层厚度(单位:m).由于客观原因,有些点无法测量煤层厚度,用‘—’标出,其中每一网格均为100mx100m的小矩形,试根据这些数据,用不同的方法估算该矩形区域煤矿的储藏量(体积)。

表14.9 煤层厚度

A

B

C

D

E

F

G

H

I

J

K

1

—

—

12.5

13.5

17.2

—

8.8

14.7

8.0

13.0

—

2

—

—

—

15.6

18.2

13

6.4

8.9

9.2

11.7

—

3

—

12

13.5

13.5

17.8

16.9

13.2

—

—

—

—

4

7.5

12.6

14.9

18.7

17.7

17.5

14.7

13

—

—

6.5

5

8.9

7.8

12.4

13.5

15.7

17.6

11.7

9.6

9.2

9.5

8.6

6

—

—

—

13.7

13.6

16.5

12.5

8.7

9.7

—

—

7

—

—

8.6

11.8

12.5

11.3

13.4

—

—

—

—

首先，我对该题有一个疑问：图中只给出了7*11个数据（姑且将‘—’也称为数据），每个数据是网格交点处的煤层厚度，而题中又说该平面区域为1100mx700m，可是表中的数据怎么画出1100mx700m的区域呢（十一个点只有十个线段啊），那该怎么办呢？

姑且每行都拟合让其拥有7*12个数据，把平面区域改为1100mx600m，这样才能有效地估计煤矿的仓储量。

clc;clear;
format long
x=0:100:1100;
y=0:100:600;
z=xlsread("C:\Users\dell\Desktop\煤层厚度.xlsx");
zz=[];
for i=1:7l=length(find(isnan(z(i,:))==0));ll=find(isnan(z(i,:))==0);lll=z(i,:);p=polyfit(1:l,lll(ll),2);m=polyval(p,1:12);if length(find(m<0))>0p=polyfit(1:l,lll(ll),1);m=polyval(p,1:12);end      zz=[zz;m];
end
[x,y]=meshgrid(x,y);
x0=0:5:1100;
y0=0:5:600;
[xx0,yy0]=meshgrid(x0,y0);
zz0=interp2(x,y,zz,xx0,yy0,'cubic');
mesh(xx0,yy0,-zz0);
v=0;
for k=1:length(x0)-1for i=1:length(y0)-1h=(zz0(i,k)+zz0(i+1,k)+zz0(i,k+1)+zz0(i+1,k+1))/4;v=v+25*h;end
end
v

v =

     8.597807745430814e+06

本文由作者自创，由于时间原因，难免出现些许错误，还请大家多多指正。创作不易，请大家多多支持。

本文来自互联网用户投稿，该文观点仅代表作者本人，不代表本站立场。本站仅提供信息存储空间服务，不拥有所有权，不承担相关法律责任。如若转载，请注明出处：http://www.mzph.cn/news/689151.shtml

如若内容造成侵权/违法违规/事实不符，请联系多彩编程网进行投诉反馈email:809451989@qq.com，一经查实，立即删除！

相关文章

数据在内存中的存储以及百度笔试题

数据在内存中的存储以及百度笔试题

目录 一.整型家族 什么是大小端存储&#xff08;百度笔试题&#xff09; 大端字节序存储 小端字节序存储 为什么要讨论大小端字节序存储 写一个程序判断是大端还是小端存储&#xff08;百度笔试题&#xff09; 思路&#xff1a;用1去判断&#xff0c;如果返回1则是小端&a…
阅读更多...
腾讯云4核8G12M服务器支持多少人在线？

腾讯云4核8G12M服务器支持多少人在线？

4核8G服务器支持多少人同时在线访问&#xff1f;阿腾云的4核8G服务器可以支持20个访客同时访问&#xff0c;关于4核8G服务器承载量并发数qps计算测评&#xff0c;云服务器上运行程序效率不同支持人数在线人数不同&#xff0c;公网带宽也是影响4核8G服务器并发数的一大因素&…
阅读更多...
Open CASCADE学习|曲线向曲面投影

Open CASCADE学习|曲线向曲面投影

在三维空间中&#xff0c;将曲线向曲面投影通常涉及复杂的几何计算。这个过程可以通过多种方法实现&#xff0c;但最常见的是使用数学和几何库&#xff0c;如OpenCASCADE&#xff0c;来处理这些计算。 在OpenCASCADE中&#xff0c;投影曲线到曲面通常涉及以下步骤&#xff1a;…
阅读更多...
elementui 中 el-date-picker 控制选择当前年之前或者之后的年份

elementui 中 el-date-picker 控制选择当前年之前或者之后的年份

文章目录 需求分析 需求 对 el-date-picker控件做出判断控制 分析 给 el-date-picker 组件添加 picker-options 属性&#xff0c;并绑定对应数据 pickerOptions html <el-form-item label"雨量年份&#xff1a;" prop"date"><el-date-picker …
阅读更多...
开源模型应用落地-工具使用篇-SLB（二）

开源模型应用落地-工具使用篇-SLB（二）

一、前言 通过学习"开源模型应用落地"系列文章&#xff0c;我们成功地建立了一个完整可实施的AI交付流程。现在&#xff0c;我们要引入负载均衡&#xff0c;以提高我们的AI服务的性能和故障转移能力。本文将详细介绍如何使用腾讯云的负载均衡技术来将我们的AI服务部署…
阅读更多...
洛谷 P1150 Peter 的烟

洛谷 P1150 Peter 的烟

参考代码and代码解读 #include<iostream> using namespace std; int main() { int n,k,nonu; //n烟的数量&#xff0c;k需要多少根烟头换一支烟&#xff0c;nonu记录烟头的个数 cin>>n>>k; int sumn; //一开始就能吸n支烟 nonusum; …
阅读更多...
深度学习之pytorch实现线性回归

深度学习之pytorch实现线性回归

度学习之pytorch实现线性回归 pytorch用到的函数torch.nn.Linearn()函数torch.nn.MSELoss()函数torch.optim.SGD() 代码实现结果分析 pytorch用到的函数 torch.nn.Linearn()函数 torch.nn.Linear(in_features, # 输入的神经元个数out_features, # 输出神经元个数biasTrue # 是…
阅读更多...
Jlink+OpenOCD+STM32 Vscode 下载和调试环境搭建

Jlink+OpenOCD+STM32 Vscode 下载和调试环境搭建

对于 Mingw 的安装比较困难&#xff0c;国内的网无法正常在线下载组件&#xff0c; 需要手动下载 x86_64-8.1.0-release-posix-seh-rt_v6-rev0.7z 版本的软件包&#xff0c;添加环境变量&#xff0c;并将 mingw32-make.exe 名字改成 make.exe。 对于 OpenOCD&#xff0c;需要…
阅读更多...
mysql 2-18

mysql 2-18

加密与解密函数 其他函数 聚合函数 三者效率 GROUP BY HAVING WHERE和HAVING的区别 子查询 单行子查询和多行子查询 单行比较操作符 多行比较操作符 把平均工资生成的结果当成一个新表 相关子查询 EXISTS 一条数据的存储过程 标识符命名规则 创建数据库 MYSQL的数据类型 创建表…
阅读更多...
利用ChatGPT进行数据分析——如何提出一个好的prompt

利用ChatGPT进行数据分析——如何提出一个好的prompt

利用ChatGPT进行数据分析——如何提出一个好的prompt ​ 使用ChatGPT时&#xff0c;能否得到一个好的输出结果&#xff0c;关键在于能否提出好的prompt。 1.1 prompt是什么 ​ 现在大家把向ChatGPT输入的内容称作prompt&#xff08;提示&#xff09;&#xff0c;它的作用是引…
阅读更多...
Docker部署Java项目

Docker部署Java项目

打包 添加 <!-- 打包插件--> <build><plugins><plugin><groupId>org.springframework.boot</groupId><artifactId>spring-boot-maven-plugin</artifactId></plugin></plugins> </build> 通过执行以下命令进…
阅读更多...
MySQL数据库基础（九）：SQL约束

MySQL数据库基础（九）：SQL约束

文章目录 SQL约束 一、主键约束 二、非空约束 三、唯一约束 四、默认值约束 五、外键约束&#xff08;了解&#xff09; 六、总结 SQL约束 一、主键约束 PRIMARY KEY 约束唯一标识数据库表中的每条记录。主键必须包含唯一的值。主键列不能包含 NULL 值。每个表都应该有…
阅读更多...
String为什么是不可变的？

String为什么是不可变的？

原因一、String字符串类型的数据结构 字符串在String类的内部是用一个char[]数组表示的,而这个数组使用final关键字修饰的&#xff0c;所以不能修改。 举例说明&#xff1a; String ip"127.0.0.1"; String retip.replace(".","#"); System.out.…
阅读更多...
BUGKU-WEB 留言板1

BUGKU-WEB 留言板1

题目描述 题目截图如下&#xff1a; 进入场景看看&#xff1a; 解题思路 之间写过一题类似的&#xff0c;所以这题应该是有什么不同的那就按照之前的思路进行测试试试提示说&#xff1a;需要xss平台接收flag&#xff0c;这个和之前说的提示一样 相关工具 xss平台&#xf…
阅读更多...
银河麒麟操作系统自动同步时间更新

银河麒麟操作系统自动同步时间更新

1、银河麒麟操作系统基于Centos8的&#xff0c;因centos8取消了ntp服务器&#xff0c;所以导致之前使用ntpdate命令无法同步时间 2、centos默认使用chrony模块来进行同步时间 3、修改chrony配置同步时间服务器 vim /etc/chrony.conf 4、目前使用的是阿里云的时间服务器&…
阅读更多...
深入浅出熟悉OpenAI最新大作Sora文生视频大模型

深入浅出熟悉OpenAI最新大作Sora文生视频大模型

蠢蠢欲动&#xff0c;惴惴不安&#xff0c;朋友们我又来了&#xff0c;这个春节真的过的是像过山车&#xff0c;Gemini1.5 PRO还没过劲&#xff0c;OpenAI又放大招&#xff0c;人类真的要认输了吗&#xff0c;让我忍不住想要再探究竟&#xff0c;到底是什么让文生视频发生了质的…
阅读更多...
流星蝴蝶剑之七夜听雪中文版下载

流星蝴蝶剑之七夜听雪中文版下载

软件介绍&#xff1a; 中文名称: 流星蝴蝶剑七夜听雪 英文名称: Meteor 游戏类型: 3D武侠格斗 发行时间: 2002年08月 制作发行: 流星江湖悠悠客栈 语言 :中文 配置要求: 操作系统&#xff1a;Windows 95 / 98 / Me / 2000 / XP 最低配置 CPU&#xff1a;Pentium II 450MHz 以上…
阅读更多...
单向/双向V2G环境下分布式电源与电动汽车充电站联合配置方法(matlab代码)

单向/双向V2G环境下分布式电源与电动汽车充电站联合配置方法(matlab代码)

目录 1 主要内容 目标函数 电动汽车负荷建模 算例系统图 程序亮点 2 部分代码 3 程序结果 4 下载链接 1 主要内容 该程序复现博士文章《互动环境下分布式电源与电动汽车充电站的优化配置方法研究》第五章《单向/双向V2G环境下分布式电源与电动汽车充电站联合配置方法》…
阅读更多...
0206-1-网络层

0206-1-网络层

第 4 章 网络层 网络层提供的两种服务 虚电路服务 数据报服务 概要: 虚电路服务与数据报服务的对比 网际协议 IP 网际协议 IP 是 TCP/IP 体系中两个最主要的协议之一。与 IP 协议配套使用的还有四个协议&#xff1a; 地址解析协议 ARP (Address Resolution Protocol)逆地…
阅读更多...
数据分析师SQL面试准备(part2)

数据分析师SQL面试准备(part2)

1. 2. 3. 4.
阅读更多...
推荐文章
最新文章

Copyright @ 2022~2023