数值分析实验报告 matlab,数值分析matlab实验报告..doc
数值分析matlab实验报告.
实验2.1 多项式差值的振荡现象
一、实验内容
,考虑区间[-1,1]的一个等距划分,分点为,i=0,1,2,...,n,.其中,li(x),i=0,1,2,...,n是Lagrange插值基函数.
选择不断增大的分点数目n=2,3,...,画出原函数f(x)及插值多项式函数Ln(x)在[-1,1]上的图像,比较并分析实验结果.
选择其他的函数,例如定义在区间[-5,5]上的函数,,重复上述的实验看其结果如何.
二、实验程序
function chapter2
promps={'请选择试验函数,若选f(x),请输入f,若选好h(x),请输入h,若选g(x),请输入g:'};
result=inputdlg(promps,'charpt 2',1,{'f'});
Nb_f=char(result);
if(Nb_f~='f'&&Nb_f~='h'&&Nb_f~='g')
errordlg('试验函数选择错误!');
return;
end
result=inputdlg({'请输入插值多项式的次数N:'},'charpt_2',1,{'10'});
Nd=str2num(char(result));
if(Nd<1)
errordlg('插值多项式的次数输入错误!');
return;
end
switch Nb_f
case'f'
f=inline('1./(1+25*x.^2)');a=-1;b=1;
case'h'
f=inline('x./(1+x.^4)');a=-5;b=5;
case'g'
f=inline('atan(x)');a=-5;b=5;
end
x0=linspace(a,b,Nd+1);y0=feval(f,x0);
x=a:0.1:b;y=Lagrange(x0,y0,x);
clf;
fplot(f,[a b],'rx');
hold on;
plot(x,y,'b--');
xlabel('x');ylabel('y=f(x) x and y=Ln(x) --');
2.Lagrange函数
function y=Lagrange(x0,y0,x)
n=length(x0);
m=length(x);
for i=1:m
z=x(i);
s=0;
for k=1:n
p=1.0;
for j=1:n
if (j~=k)
p=p.*(z-x0(j))/(x0(k)-x0(j));
end
end
s=s+p*y0(k);
end
y(i)=s;
end
三、实验结果及分析
1) 选择不断增大的分点数目n,原函数f(x)及插值多项式函数Ln(x)在[-1,1]上的图像。
随着提高插值多项式次数,可以提高逼近的精度,但是次数的增加,在区间两端点附近与原函数偏离很远,即出现了Runge现象。
2) 选择不断增大的分点数目n,原函数h(x)及插值多项式函数Ln(x)在[-1,1]上的图像。
选择不断增大的分点数目n,原函数g(x)及插值多项式函数Ln(x)在[-1,1]上的图像。
同样,随着提高插值多项式次数,可以提高逼近的精度,但是次数的增加,在区间两端点附近与原函数偏离很远,即出现了Runge现象。
实验3.1最小二乘拟合
一、实验内容
为基的多项式最小二乘拟合程序,并用于对下表中数据作3次多项式最小二乘拟合.
xi-1.0-0.50.00.51.01.52.0yi-4.447-0.4520.5510.048-0.4470.5494.552取权数,求拟合曲线中的参数、平方误差,并作离散数据的拟合函数的图形.
二、实验程序
function chapter3
x0=-1:0.5:2;
y0=[-4.447 -0.452 0.551 0.048 -0.447 0.549 4.552];
n=3;
alph=polyfit(x0,y0,n);
y=polyval(alph,x0);
r=(y0-y)*(y0-y)';
x=-1:0.01:2;
y=polyval(alph,x);
plot(x,y,'k--');
xlabel('x');ylabel('y0 * and polyfit. y-.');
hold on;
plot(x0,y0,'*');
title('离散数据的多项式拟合');
grid
数值分析实验报告 matlab,数值分析matlab实验报告..doc相关推荐
- 数值分析复化求积matlab,MATLAB数值分析实验二(复合梯形、辛普森和龙贝格求积,以及二重积分计算等)...
1.理解如何在计算机上使用数值方法计算定积近似值; 2.学会复合梯形.复合Simpson和龙贝格求积分公式的编程与应用. 3.探索二重积分在矩形区域的数值积分方法. 佛山科学技术学院 实 验 报 告 ...
- 控制工程matlab实验报告小结,控制工程MATLAB实验报告.doc
控制工程MATLAB实验报告 浙江科技学院 机电系统开放性实验 设计题目: MATLAB控制机电系统实分析 专 业: 机 械 设 计 制 造及 其 自 动 化 班 级: 机 制 122 学生姓名: 王 ...
- MATLAB环境认识实验报告,matlab实验报告ljg
matlab实验报告ljg 第 1 页 共 3 页电子信息工程学系实验报告课程名称: MATLAB 程序设计 实验项目名称: MATLAB 数值计算 实验时间:2011.4.26 班级:电信 081 ...
- matlab叶子分割实验,基于MATLAB进行树叶面积测量实验报告
科学计算与数据处理实验报告学 号 姓名 实验 名称 基于Matlab的植物叶面积数字摄影图像处理 实验目的1.学习用matlab解决实际问题; 2.学习用matlab对图像进行背景去除; 3.学习用m ...
- MATLAB语言算法实验报告,机械工程实验——matlab实验报告.doc
机械工程实验教学中心 - PAGE 20 - 机械工程实验教学中心 实验指导书 实验名称 基于Matlab的信号处理实验 课程名称 自选综合实验 一.实验目的及要求 实验目的 通过基于Matlab的信 ...
- 一般单纯形法的matlab程序,实验报告(单纯形法的matlab程序)
<实验报告(单纯形法的matlab程序)>由会员分享,可在线阅读,更多相关<实验报告(单纯形法的matlab程序)(5页珍藏版)>请在人人文库网上搜索. 1.实验一:线性规划单 ...
- matlab操作入门实验报告,matlab操作实验报告
实验一matlab基本操作 一.实验目的 熟悉matlab的安装与启动:熟悉matlab用户界面:熟悉matlab功能.建模元素:熟悉matlab优化建模过程. 二.实验设备与工具 1.计算机 2.m ...
- 实验一熟悉matlab环境,数字信号处理报告实验一:熟悉MATLAB环境.doc
数字信号处理报告实验一:熟悉MATLAB环境.doc 实验一熟悉MATLAB环境一 实验目的1. 熟悉MATLAB的主要操作命令.2. 学会简单的矩阵输入和数据读写.3. 掌握简单的绘图命令.4. 用 ...
- matlab m文件的编写,Matlab实验报告(四)M文件的编写.doc
Matlab实验报告(四)M文件的编写.doc 实验目的 学习MATLAB屮的关系运算和逻辑运算,掌握它们的表达形式和川法. 掌握MATLAB中的选择结构和循环结构. 学会用MATLAB进行M文件的编 ...
- matlab矩阵处理实验报告,matlab实验报告一二三
<matlab实验报告一二三>由会员分享,可在线阅读,更多相关<matlab实验报告一二三(37页珍藏版)>请在金锄头文库上搜索. 1.2015 秋秋 2013 级级MATLA ...
最新文章
- 某程序员为方便老婆工作,写了一些小工具给老婆用!如今老婆要离职,公司老板却要求必须把工具留下!网友:跟他收费!...
- java 分配cpu使用率_Java-JVM-CPU使用率达到100%
- (009) java后台开发之堆和栈的区别
- 递归概念分类注意事项
- 不支持对系统目录进行即席更新_「目录」让你的文档结构一目了然
- 基于上一篇AS项目依赖库问题的优化解决方案
- mysql身份验证失败_SMTP身份验证失败PAM-MySQL无法进行身份验证
- Python进阶|聊聊异常处理
- 荣耀20系列国行版正式发布 拍照全球第三 屏幕开孔全球最小
- [渝粤教育] 中国地质大学 中外美术史 复习题
- [转]安装完成后机器重新启动运行 JBuilder 时却每次只是看到 L
- html5开发播放器,larkplayer: 插件化的 HTML5 播放器
- python六大数据类型的定义_python六大类标准数据类型和数据类型转换
- 华为交换之间静态路由配置实例
- C盘太小,调整磁盘分区大小
- 机器学习入门实战——决策树算法实战Titanic数据集
- Python学习笔记(5),Battleship 游戏
- 从今天起,我决定去掉cnzz统计
- 计算机领域影响因子3.5什么水平,计算机学科SCI(包括SCI-E)影响因子排名
- 教你在WORD中如何制作三线表