急!如何用matlab矩形法求解定积分?
来源:学生作业帮 编辑:大师作文网作业帮 分类:综合作业 时间:2024/11/18 13:11:55
急!如何用matlab矩形法求解定积分?
我现需要用MATLAB矩形法,求解函数sinx在区间0到2pi的定积分数值,题目要求是将sinx图形在区间0到2pi内按矩形法分为1等份,2等份,3等份,4等份,5等份.直至1000等份,总共1000种等份分法,采用MATLAB计算出各种等份分法 函数sinx在区间内积分的数值解,将各个数值解用线连接起来,作出曲线图.将得到的曲线与标准值(及用高等数学知识算出的结果,好像是0...)所在的直线进行比较. 并通过图形分析随着等份数的增加,各数值解与标准解的关系.
请各位高手帮我编出计算的程序 和 作曲线图的程序!
本人感激不尽!
越详细越好!麻烦能够附带编程旁白,帮助我理解你的程序!
希望尽快帮忙!明天就要了!
我现需要用MATLAB矩形法,求解函数sinx在区间0到2pi的定积分数值,题目要求是将sinx图形在区间0到2pi内按矩形法分为1等份,2等份,3等份,4等份,5等份.直至1000等份,总共1000种等份分法,采用MATLAB计算出各种等份分法 函数sinx在区间内积分的数值解,将各个数值解用线连接起来,作出曲线图.将得到的曲线与标准值(及用高等数学知识算出的结果,好像是0...)所在的直线进行比较. 并通过图形分析随着等份数的增加,各数值解与标准解的关系.
请各位高手帮我编出计算的程序 和 作曲线图的程序!
本人感激不尽!
越详细越好!麻烦能够附带编程旁白,帮助我理解你的程序!
希望尽快帮忙!明天就要了!
昨天大牛:tanarri
不是已经回答了吗?
你到“我的提问”里面看一下吧.
不过好像等分数变了,我的方法是比较笨的方法
n=1:1000; %等分数
x0=0; %积分下界
x1=2*pi; %积分上界
Ans=0; %解析值
rect=[]; %用来存放各种等分下的矩形法结果
for nn=n %第一次是1等分,第二次是2等分,……
dx=(x1-x0)/nn; %nn等分时,x间隔
x=linspace(x0,x1,nn+1); %nn等分的点
y=sin(x); %计算出各等分点的函数值
temp=sum(y)*dx; %nn等分,积分值
rect=[rect,temp]; %存到数组里
end
plot(n,rect); %从图中可以看出,任何一个等分与解析值都比较接近,主要存在一些涨落(波动)
%注意,纵轴有一个10E-15因子,说明误差都比较小
res=rect-Ans; %误差
max(res) %最大误差
不是已经回答了吗?
你到“我的提问”里面看一下吧.
不过好像等分数变了,我的方法是比较笨的方法
n=1:1000; %等分数
x0=0; %积分下界
x1=2*pi; %积分上界
Ans=0; %解析值
rect=[]; %用来存放各种等分下的矩形法结果
for nn=n %第一次是1等分,第二次是2等分,……
dx=(x1-x0)/nn; %nn等分时,x间隔
x=linspace(x0,x1,nn+1); %nn等分的点
y=sin(x); %计算出各等分点的函数值
temp=sum(y)*dx; %nn等分,积分值
rect=[rect,temp]; %存到数组里
end
plot(n,rect); %从图中可以看出,任何一个等分与解析值都比较接近,主要存在一些涨落(波动)
%注意,纵轴有一个10E-15因子,说明误差都比较小
res=rect-Ans; %误差
max(res) %最大误差