案例 黄河小浪底调水调沙问题
日期 | 6.29 | 6.3 | 7.1 | 7.2 | ||||
时间 | 0 | 12 | 24 | 36 | 48 | 60 | 72 | 84 |
水流 | 1800 | 1900 | 2100 | 2200 | 2300 | 2400 | 2500 | 2600 |
含沙量 | 32 | 60 | 75 | 85 | 90 | 98 | 100 | 102 |
日期 | 7.3 | 7.4 | 7.6 | 7.8 | ||||
时间 | 96 | 108 | 120 | 132 | 144 | 156 | 168 | 180 |
水流 | 2650 | 2700 | 2720 | 2650 | 2600 | 2500 | 2300 | 2200 |
含沙量 | 108 | 112 | 115 | 116 | 118 | 120 | 118 | 105 |
日期 | 7.9 | 7.11 | 7.12 | 7.13 | ||||
时间 | 192 | 204 | 216 | 228 | 240 | 252 | 264 | 276 |
水流 | 2000 | 1850 | 1820 | 1800 | 1750 | 1500 | 1000 | 900 |
含沙量 | 80 | 60 | 50 | 30 | 26 | 20 | 8 | 5 |
现在,根据实验数据,建立数学模型研究一下问题: (1)给出任意时刻排沙量和总排沙量的方法; (2)给出排沙量和水流量的关系。
问题二代码
A=xlsread('C:\Users\86180\Desktop\含沙量.xls');
y=A(3,:);
r=A(4,:);
x=y.*r;
figure(1);
plot(y,x,'*'); //离散点
a=polyfit(y/1000,x/100,2);
y1=y/1000;
x1=a(1)*y1.^2+a(2)*y1+a(3); //二次多项式函数
x2=100*x1; //可以去掉
figure(2);
plot(y,x,'*',y,x2,'+-');
legend('实测值','拟合值');