声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3034|回复: 1

[其他相关] 用三次样条插值法解非线性方程组,求高手matlab编程

[复制链接]
发表于 2006-7-1 22:56 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
<P ><FONT face="Times New Roman">2</FONT>.已知数据点<v:shapetype><FONT face="Times New Roman"> <v:stroke joinstyle="miter"></v:stroke><v:formulas><v:f eqn="if lineDrawn pixelLineWidth 0"></v:f><v:f eqn="sum @0 1 0"></v:f><v:f eqn="sum 0 0 @1"></v:f><v:f eqn="prod @2 1 2"></v:f><v:f eqn="prod @3 21600 pixelWidth"></v:f><v:f eqn="prod @3 21600 pixelHeight"></v:f><v:f eqn="sum @0 0 1"></v:f><v:f eqn="prod @6 1 2"></v:f><v:f eqn="prod @7 21600 pixelWidth"></v:f><v:f eqn="sum @8 21600 0"></v:f><v:f eqn="prod @7 21600 pixelHeight"></v:f><v:f eqn="sum @10 21600 0"></v:f></v:formulas><v:path connecttype="rect" gradientshapeok="t" extrusionok="f"></v:path><lock aspectratio="t" v:ext="edit"></lock></FONT></v:shapetype>,作三次插值样条函数,然后求x~<FONT face="Times New Roman">处y~=s(x~)的值<v:shape> <v:imagedata></v:imagedata></v:shape></FONT></P>
<P >边界条件:</P>
<P ><FONT face="Times New Roman">1)  </FONT>自然三次样条:</P>
<P ><FONT face="Times New Roman">2)  </FONT>给出两端点的一阶导数<v:shape><FONT face="Times New Roman"> <v:imagedata></v:imagedata></FONT></v:shape></P>
<P >说明:基本方程组用<I ><FONT face="Times New Roman">m</FONT></I>表示法;方程组求解用追赶法。</P>
<P >Example:<p></p></P>
<P >给定数据表</P>
<TABLE  cellSpacing=0 cellPadding=0 border=0>

<TR >
<TD  vAlign=top width=68>
<P ><I ><FONT face="Times New Roman">x<SUB>i</SUB></FONT></I></P></TD>
<TD  vAlign=top width=68>
<P ><FONT face="Times New Roman">0</FONT></P></TD>
<TD  vAlign=top width=68>
<P ><FONT face="Times New Roman">2</FONT></P></TD>
<TD  vAlign=top width=68>
<P ><FONT face="Times New Roman">4</FONT></P></TD>
<TD  vAlign=top width=68>
<P ><FONT face="Times New Roman">6</FONT></P></TD>
<TD  vAlign=top width=68>
<P ><FONT face="Times New Roman">10</FONT></P></TD></TR>
<TR >
<TD  vAlign=top width=68>
<P ><I ><FONT face="Times New Roman">y<SUB>i</SUB></FONT></I></P></TD>
<TD  vAlign=top width=68>
<P ><FONT face="Times New Roman">0</FONT></P></TD>
<TD  vAlign=top width=68>
<P ><FONT face="Times New Roman">16</FONT></P></TD>
<TD  vAlign=top width=68>
<P ><FONT face="Times New Roman">36</FONT></P></TD>
<TD  vAlign=top width=68>
<P ><FONT face="Times New Roman">54</FONT></P></TD>
<TD  vAlign=top width=68>
<P ><FONT face="Times New Roman">82</FONT></P></TD></TR></TABLE>
<P ><FONT face="Times New Roman">a.      </FONT>确定自然三次样条函数,并求<FONT face="Times New Roman"><I >S</I>(4.65)</FONT>;</P>
<P ><FONT face="Times New Roman">b.     </FONT>确定满足边界条件s'(0)=8  和 s'(0)的三次样条函数,并求<FONT face="Times New Roman"><I >S</I>(3)</FONT>和<FONT face="Times New Roman"><I >S</I>(8).</FONT></P>
回复
分享到:

使用道具 举报

发表于 2006-7-4 15:02 | 显示全部楼层
这样的问题,可以先将系数矩阵A生成,再用
x=A\b得到结果

matlab的方程求解自动会兼顾到效率


Sample TextSample Text
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-1 08:32 , Processed in 0.073176 second(s), 18 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表