MATLAB作业拉格朗日三阶样条插值函数牛顿插值

Lagrange插值M函数
syms x
x0=[0,1,2];
y0=[1,2,3];
n=length(x0);
for i=1:n
a=1;
for j=1:n
if j~=i
a=expand(a*(x-x0(j)));
end
end
vdisk
b=1;
for k=1:n
if k~=i
b=b*(x0(i)-x0(k));
end
end
A(i)=expand(a/b);
end
L=0;
for p=1:n
L=L+y0(p)*A(p);
end
L
>> Language
L =x + 1
三阶样条插值M函数
function m=naspline(x,y,dy0,dyn,xx)
n=length(x)-1;
h=diff(x);lemda=h(2/n)./(h(1:n-1)+h(2:n));mu=1-lemda;
g=3*(lemda.*diff(y(1:n))./h(1:n-1)+mu.*diff(y(2:n+1))./h(2:n)); g(1)=g(1)-lemda(1)*dy0;g(n-1)=g(n-1)-mu(n-1)*dyn;
dy=nachase(lemda,2*ones(1:n-1),mu,g);
m=[dy0;dy;dyn];
if nargin>=5
s=zeros(size(xx));
for i=1:n
if i==1,
kk=find(xx<=x(2));
elseif i==n
kk=find(xx>x(n));
else
kk=find(xx>x(i)&xx<=x(i+1));
end
xbar=(xx(kk)-x(i))/h(i);
s(kk)=alpha0(xbar)*y(i)+alpha1(xbar)*y(i+1)+...
h(i)*beta0(xbar)*m(i)+h(i)*beta1(xbar)*m(i+1);
end
m=s;
end
function x=nachase(a,b,c,d)
n=length(a);
for k=2:n
b(k)=b(k)-a(k)/b(k-1)*c(k-1);
d(k)=d(k)-a(k)/b(k-1)*d(k-1);
end
x(n)=d(n)/b(n);
for k=n-1:-1:1
x(k)=(d(k)-c(k)*x(k+1))/b(k);
end
x=x(:);
function y=alpha0(x)
y=2*x.^3-3*x.^2+1;
function y=alpha1(x)
y=-2*x.^3+3*x.^2;
function y=beta0(x)
y=x.^3-2*x.^2+x;
function y=beta1(x)
y=x.^3-x.^2;
电力线网络摄像机
naspline([-1 0 1],[-1 0 1],0,-1)
ans =
1.7500
-1.0000
>> naspline([-1 0 1],[-1 0 1],0,-1,-1:0.25:1)
ans =
-1.0000 -0.9258 -0.7188 -0.4023 0 0.4492 0.8438    1.0664
1.0000
ans =
Columns 1 through 5
-1.0000 -0.9258 -0.7188 -0.4023 0 Columns 6 through 9
0.4492 0.8438    1.0664    1.0000
牛顿插值多项式
function yi=Newton(x,y,xi)陶瓷咖啡壶
n=length(x);
m=length(y);
if n~=m
error;
return;
end
%计算均差表Y
Y=zeros(n);
Y(:,1)=y';
for k=1:n-1
for i=1:n-k
if abs(x(i+k)-x(i))<eps
error;
return;
end
Y(i,k+1)=(Y(i+1,k)-Y(i,k))/(x(i+k)-x(i));
end
end
%计算牛顿插值公式
yi=0;
for i=1:n
z=1;
FOSYfor k=1:i-1
z=z*(xi-x(k));
end
yi=yi+Y(1,i)*z;
end
End
>>format compact
>> x=pi*[1/6 1/4 1/3];y=[0.5 0.7071 0.866];xi=2*pi/9;
0.550.60.650.70.750.450.50.550.60.650.7
0.75
>> yi=Newton(x,y ,xi)
yi =
0.6434
>>fplot(…sin ‟,[pi/6 pi/4 pi/3]);hold on;
>>plot(x,y ,‟o ‟xi,0.6434,‟rv ‟);hold off;
钢格板压焊机牛顿迭代M 函数
function x=nanewton(fname,dfname,x0,e,N)
if nargin<5,N=500;end
if nargin<4,e=1e-4;end
x=x0;x0=x+2*e;k=0;
while abs(x0-x)>e&k<N,
k=k+1;
x0=x;x=x0-feval(fname,x0)/feval(dfname,x0);
disp(x)
end
if k==N,warning('已达迭代次数上限');end
>> fun=inline('x^3-x-1');dfun=inline('3*x^2-1');
>> nanewton(fun,dfun,1.5,0.5e-3)
1.3478
1.3252
1.3247
ans =
1.3247
0.550.60.650.70.750.450.5
0.55
0.6
0.65
0.7
0.75
生物总黄酮是指黄酮类化合物,是一大类天然产物,广泛存在于植物界,是许多中草药的有效成分。在自然界中最常见的是黄酮和黄酮醇,其它包括双氢黄(醇)、异黄酮、双黄酮、黄烷醇、查尔酮、橙酮、花苷及新黄酮类等。
简介
近年来,由于自由基生命科学的进展,使具有很强的抗氧化和消除自由基作用类黄酮受到空前的重视。类黄酮参与了磷酸与花生四烯酸的代谢、蛋白质的磷酸化、钙离子的转移、自由基的清除、抗氧化活力的增强、氧化还原作用、螯合作用和基因的表达。它们对健康的好处有:(  1 )抗炎症(  2 )抗过敏(  3 )抑制细菌(  4 )抑制寄生虫(  5 )抑制病毒(  6 )防治肝病(7 )防治血管疾病(8 )防治血管栓塞(9 )防治心与脑血管疾病(10 )抗肿瘤(11 )抗化学毒物等。天然来源的生物黄酮分子量小,能被人体迅速吸收,能通过血脑屏障,能时入脂肪组织,进而体现出如下功能:消除疲劳、保护血管、防动脉硬化、扩张毛细血管、疏通微循环、活化大脑及其他脏器细胞的功能、抗脂肪氧化、抗衰老。近年来国内外对茶多酚、银杏类黄酮等的药理和营养性的广泛深入的研究和临床试验,证实类黄酮既是药理因子,又是重要的营养因子为一种新发现的营养素,对人体具有重要的生理
保健功效。目前,很多著名的抗氧化剂和自由基清除剂都是类黄酮。例如,茶叶提取物和银杏提取物。葛根总黄酮在国内外研究和应用也已有多年,其防治动脉硬化、治偏瘫、防止大脑萎缩、降血脂、降血压、防治糖尿病、突发性耳聋乃至醒酒等不乏数例较多的临床报告。从法国松树皮和葡萄籽中提取的总黄酮" 碧萝藏"-- (英文称PYCNOGENOL )在欧洲以不同的商品名实际行销应用25 年之久,并被美国FDA 认可为食用黄酮类营养保健品,所报告的保健作用相当广泛,内用称之为" 类维生素" 或抗自由基营养素,外用称之为" 皮肤维生素" 。进一步的研究发现碧萝藏的抗氧化作用比VE 强50 倍,比VC 强20 倍,而且能通过血脑屏障到达脑部,防治中枢神经系统的疾病,尤其对皮肤的保健、年轻化及血管的健康抗炎作用特别显著。在欧洲碧萝藏已作为保健药物,在美国作为膳食补充品(相当于我国的保健食品),风行一时。随着对生物总黄酮与人类营养关系研究的深入,不远的将来可能证明黄酮类化合物是人类必需的微营养素或者是必需的食物因子。性状:片剂。
功能主治与用法用量
downlinker
功能主治:本品具有增加脑血流量及冠脉血流量的作用,可用于缓解高血压症状(颈项强痛)、心绞痛及突发性耳聋,有一定疗效。用法及用量:口服:每片含总黄酮60mg,每次5片,1日3次。
不良反应与注意
不良反应和注意:目前,暂没有发现任何不良反应.

本文发布于:2024-09-24 11:26:01,感谢您对本站的认可!

本文链接:https://www.17tex.com/tex/1/208208.html

版权声明:本站内容均来自互联网,仅供演示用,请勿用于商业和其他非法用途。如果侵犯了您的权益请与我们联系,我们将在24小时内删除。

标签:防治   黄酮   类黄酮   作用
留言与评论(共有 0 条评论)
   
验证码:
Copyright ©2019-2024 Comsenz Inc.Powered by © 易纺专利技术学习网 豫ICP备2022007602号 豫公网安备41160202000603 站长QQ:729038198 关于我们 投诉建议