| 查看: 1365 | 回复: 7 | ||||
| 本帖产生 1 个 博学EPI ,点击这里进行查看 | ||||
21056651木虫 (小有名气)
|
[求助]
请教matlap
|
|||
|
2.(25分)给定数据 x 1.30 1.32 1.34 1.36 1.38 f(x) 3.60210 3.90330 4.25560 4.67344 5.17744 请用MATLAB求最小二乘拟合二次多项式,并作出拟合曲线。再用组合Simpson公式计算 的近似值。 3.(15分)用Gauss-Seidel迭代通过MATLAB求解下列线性方程组的解。 -2x1+2x2+3x3=12 -4x1+2x2+x3=12 x1+2x2+3x3=16 4.(25分)给定数据 x 0 2 3 5 f(x) 1 -3 -4 2 通过MATLAB求出f(x)的3次Lagrange插值多项式和3次Newton插值多项式。 5.(25分)用MATLAB计算填料吸收塔的总传质单元数NOG 逆流操作的填料吸收塔,在温度20℃、压力101.325kPa的条件下,用水洗涤含有5.5%SO2的空气,使SO2下降到0.5%。当液气比为40时,求取以气相浓度为基准的总传质单元数NOG。若以气相浓度为基准的总传质单元高度HOG=0.69m,进而求填料塔高度。进塔水中不含SO2,20℃时水对SO2的溶解度数据如表中。 20℃时水对SO2的溶解度数据 x×103 1.96 1.40 0.846 0.562 0.422 0.281 0.141 0.056 y×103 51.3 34.2 18.6 11.2 7.63 4.21 1.58 0.658 要求写出程序和运算结果,选修课,实验没顾上。愿高手相助。 |
» 猜你喜欢
投稿Elsevier的杂志(返修),总是在选择OA和subscription界面被踢皮球
已经有8人回复
自荐读博
已经有7人回复
自然科学基金委宣布启动申请书“瘦身提质”行动
已经有4人回复
求个博导看看
已经有18人回复
lg2880
铁杆木虫 (正式写手)
【鬼迷日眼】
- 博学EPI: 19
- 应助: 13 (小学生)
- 金币: 6535.5
- 散金: 471
- 红花: 5
- 帖子: 724
- 在线: 287.4小时
- 虫号: 765154
- 注册: 2009-05-07
- 性别: GG
- 专业: 无机非金属类光电信息与功

2楼2010-07-03 19:55:59
shp011222
铁杆木虫 (著名写手)
- 博学EPI: 21
- 应助: 0 (幼儿园)
- 金币: 4895.8
- 散金: 20
- 红花: 3
- 帖子: 1230
- 在线: 139小时
- 虫号: 865435
- 注册: 2009-10-08
- 专业: 环境工程
3楼2010-07-03 22:08:16
21056651
木虫 (小有名气)
- 应助: 0 (幼儿园)
- 金币: 2493.8
- 散金: 30
- 红花: 1
- 帖子: 59
- 在线: 148.7小时
- 虫号: 919517
- 注册: 2009-12-03
- 性别: GG
- 专业: 环境化学
4楼2010-07-03 23:03:15
mapingke
木虫 (著名写手)
- 应助: 1 (幼儿园)
- 金币: 2215.5
- 帖子: 1289
- 在线: 77.5小时
- 虫号: 972682
- 注册: 2010-03-16
- 专业: 计算数学与科学工程计算
5楼2010-07-04 00:36:39
6楼2010-07-04 01:05:59
lcg2009
荣誉版主 (著名写手)
-

专家经验: +1 - 博学EPI: 7
- 应助: 2 (幼儿园)
- 贵宾: 1.53
- 金币: 14228.9
- 散金: 2100
- 红花: 20
- 沙发: 3
- 帖子: 2415
- 在线: 643.2小时
- 虫号: 899699
- 注册: 2009-11-10
- 性别: GG
- 专业: 知识管理
- 管辖: 教师之家

7楼2010-07-04 09:41:59
lcg2009
荣誉版主 (著名写手)
-

专家经验: +1 - 博学EPI: 7
- 应助: 2 (幼儿园)
- 贵宾: 1.53
- 金币: 14228.9
- 散金: 2100
- 红花: 20
- 沙发: 3
- 帖子: 2415
- 在线: 643.2小时
- 虫号: 899699
- 注册: 2009-11-10
- 性别: GG
- 专业: 知识管理
- 管辖: 教师之家
【答案】应助回帖
21056651(金币+5): 2010-07-04 23:14:30
|
拉格朗日函数的定义 function yi=agrange(x,y,xi) %Lagrange插值多项式,其中, %x为向量,全部的插值节点; %y为向量,插值节点处的函数值; %xi为标量,被估计函数的自变量: %yi为xi处的函数估计值. n=length(x); m=length(y); %输入的插值点与它的函数值应有相同的个数 if n~=m error('The length of X must be equal!'); return; end p=zeros(1,n); for k=1:n t=ones(1,n); for j=1:n if j~=k %输入的插值节点必须互异 if abs(x(k)-x(j))< eps error('the DATA is error!'); break; return; end end t(1,j)=(xi-x(j))/(x(k)-x(j)); t(1,k)=1; end p(k)=prod(t); end yi=sum(y.*p); 时间太忙就做这么做吧!其他题目建议楼主先百度再自己慢慢调试。 |

8楼2010-07-04 09:44:05







回复此楼