- 1、原创力文档(book118)网站文档一经付费(服务费),不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。。
- 2、本站所有内容均由合作方或网友上传,本站不对文档的完整性、权威性及其观点立场正确性做任何保证或承诺!文档内容仅供研究参考,付费前请自行鉴别。如您付费,意味着您自己接受本站规则且自行承担风险,本站不退款、不进行额外附加服务;查看《如何避免下载的几个坑》。如果您已付费下载过本站文档,您可以点击 这里二次下载。
- 3、如文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“版权申诉”(推荐),也可以打举报电话:400-050-0827(电话支持时间:9:00-18:30)。
- 4、该文档为VIP文档,如果想要下载,成为VIP会员后,下载免费。
- 5、成为VIP后,下载本文档将扣除1次下载权益。下载后,不支持退款、换文档。如有疑问请联系我们。
- 6、成为VIP后,您将拥有八大权益,权益包括:VIP文档下载权益、阅读免打扰、文档格式转换、高级专利检索、专属身份标志、高级客服、多端互通、版权登记。
- 7、VIP文档为合作方或网友上传,每下载1次, 网站将根据用户上传文档的质量评分、类型等,对文档贡献者给予高额补贴、流量扶持。如果你也想贡献VIP文档。上传文档
查看更多
实用数值计算基础
数值计算基础作业
学院:生物工程学院
专业:轻工技术与工程
学号:XXXXXXXXXXX
姓名: XXXX问题一:利用最小二乘法进行数据拟合
蛋白质是一类生物大生子,是生物体生命活动的直接承担者。现有一种微孔滤纸,可以使溶液中的蛋白质吸附在滤纸上,从而使样品中蛋白质的含量下降,现为了测定微孔滤纸吸附蛋白的能力采用考马斯亮蓝法进行测定。考马斯亮蓝G-250能够与蛋白质发生反应,溶液颜色由蓝色变为红色。现设置浓度梯度为0.02mg/ml、0.04mg/ml、0.06mg/ml、0.08mg/ml、0.10mg/ml的标准蛋白溶液,加入考马斯亮蓝后在波长为595nm处测定吸光度,每个样品测定三次,求平均值绘制标准曲线,求解曲线方程。数据记录如下:
标准蛋白浓度(mg/ml) 0 0.02 0.04 0.06 0.08 0.10
吸光值 0.000 0.168 0.270 0.406 0.532 0.646 0.000 0.166 0.272 0.407 0.529 0.639 0.000 0.163 0.261 0.412 0.532 0.634 平均值 0.000 0.166 0.268 0.408 0.531 0.640
现根据已有数据求出蛋白的标准曲线?
答:
根据描点的结果,将曲线拟合成
y = ax+b
列表如下:
i Xi yi xiyi Xi2 0 0.00 0.000 0.000 0.0000 1 0.02 0.166 0.00332 0.0004 2 0.04 0.268 0.01072 0.0016 3 0.06 0.408 0.02448 0.0036 4 0.08 0.531 0.04248 0.0064 5 0.10 0.640 0.06400 0.0100 Xi 0.3 yi 2.013 xiyi 0.145 xi2 0.022
法方程组为
解得 a=6 b= 0(2)最小二乘法线性拟合的MATLAB程序polyfit.m
function p = polyfit(x,y,n)
G = zeros(n+1,n+1);
for k = 0:n
for j = 0:n
G(k+1,j+1) = sum (x.^(k+j));
end
b(k+1) = sum (x.^k.*y);
end
a = G\b;
p = fliplr(a);
(3)程序运行结果
x=[0.00 0.02 0.04 0.06 0.08 0.10];
y=[0.000 0.166 0.268 0.408 0.531 0.640];
p =
6.335714 0.018714
所拟合曲线方程为:
y = 6.335714x +0.018714.
分析
在多项式拟合中,构得的法方程组往往是病态的,影响拟合的精确度,这时可以根据观察点构造一组正交多项式作为基函数,进行正交多项式拟合。
问题二:利用Lagrange插值多项式进行求值
在评价水产品的鲜度时挥发性盐基氮的含量是一个非常重要的指标,对虾在贮藏运输的过程中因为细菌及微生物的活动以及自身的生物化学反应使得肌体中的蛋白质和非蛋白氮分解产生氨以及胺类等碱性含氮化合物,这类化合物具有挥发性,其含量越高表明蛋白质破坏的程度越高,使得对虾的营养价值大大降低甚至产生食品安全问题。用多糖溶液处理对虾,检测其挥发性盐基氮的含量,验证多糖溶液是否能延长对虾的保质期。同时设置对照。检测第0、3、6、9、12天对虾中的挥发性盐基氮。数据记录如下表。
时间/天 0 3 6 9 12 挥发性盐基氮的含量
(mg/100g) 对照组 8.91 20.93 35.47 56.82 69.32 实验组 8.91 15.64 26.38 35.9 54.32
现根据已有数据试估计对照组第五天和实验组第七天样品中的挥发性盐基氮?
答:
取节点信息(3,20.93),(6,35.47)由Lagrange插值函数有
将t= 5代入有L1(5)=35.....
(2)取节点信息(6,26.38),(9,35.90)由Lagrange插值函数有
将t=7代入有L1(7)=29.5533
文档评论(0)