数学实验:通过Matlab软件操作,对函数性质有进一步了解 f(x)=sin^2(1.006)x 100
数学实验:通过Matlab软件操作,对函数性质有进一步了解f(x)=sin^2(1.006)x,x∈[0,π]1,在自己设定的区间内画出图形;2.分割区间成n=100,n...
数学实验:通过Matlab软件操作,对函数性质有进一步了解
f(x)=sin^2(1.006)x,x∈[0,π]
1,在自己设定的区间内画出图形;
2.分割区间成n=100,n=200和n=1000个等长子区间并且求出在每个子区间中点的函数值;
3.计算由2产生的函数值的平均值;
4.用对于n=1000的分割,由3计算出的均值,对于x解方程f(x)=均值
5.用高数有关知识解释上面3与4的结果
——————
这是数学实验报告的综合练习,明天就要交了,跪求大神帮忙…急!!!!!T_T 展开
f(x)=sin^2(1.006)x,x∈[0,π]
1,在自己设定的区间内画出图形;
2.分割区间成n=100,n=200和n=1000个等长子区间并且求出在每个子区间中点的函数值;
3.计算由2产生的函数值的平均值;
4.用对于n=1000的分割,由3计算出的均值,对于x解方程f(x)=均值
5.用高数有关知识解释上面3与4的结果
——————
这是数学实验报告的综合练习,明天就要交了,跪求大神帮忙…急!!!!!T_T 展开
3个回答
展开全部
f = @(x)(sin(1.006*x).^2);
format long;
m = [100,200,1000];
for i=1:3
%区间分n段等长子区间
n = m(i);
xn = (0:n-1)/n*pi; %xn(i)为第i段区间左端点值
d = pi/n; %子区间宽度
xnm = xn+d/2; %区间中点
ynm = f(xnm); %区间中点的函数值
S(i) = sum(ynm); %求和
end
%求平均值
Smean = sum(S)/sum(m) %函数均值
%解方程
x = fsolve(@(x) sin(1.006*x)^2-Smean,0.1)
更多追问追答
追问
这个确定是用Matlab做的吗?
追答
你太逗了
展开全部
%f(x)=sin^2(m/100)x,x属于[0,pi],m=403
%(1)-----------------
x=0:0.01:pi;
m=403;
f = sin(m*x/100).^2;
plot(x,f);
%(2)---------------
n=100;
x1 = 0:pi/n:pi;
x1 = x1 + pi/2*n;
x1 = x1(1:n-1);
f1 = sin(m*x/100).^2;
n=200;
x2 = 0:pi/n:pi;
x2 = x2 + pi/2*n;
x2 = x2(1:n-1);
f2 = sin(m*x/100).^2;
n=1000;
x3 = 0:pi/n:pi;
x3 = x3 + pi/2*n;
x3 = x3(1:n-1);
f3 = sin(m*x/100).^2;
%(3)--------------------------
m1 = (99*mean(f1)+199*mean(f2)+999*mean(f3))/(99+199+999);
disp('由2产生的函数值的平均值')
disp(m1)
%(4)-----------------------
xx = fsolve(@(x) sin(m*x/100)^2-m1,0.1);
disp('对于x解方程f(x)=均值')
disp(xx)
%(1)-----------------
x=0:0.01:pi;
m=403;
f = sin(m*x/100).^2;
plot(x,f);
%(2)---------------
n=100;
x1 = 0:pi/n:pi;
x1 = x1 + pi/2*n;
x1 = x1(1:n-1);
f1 = sin(m*x/100).^2;
n=200;
x2 = 0:pi/n:pi;
x2 = x2 + pi/2*n;
x2 = x2(1:n-1);
f2 = sin(m*x/100).^2;
n=1000;
x3 = 0:pi/n:pi;
x3 = x3 + pi/2*n;
x3 = x3(1:n-1);
f3 = sin(m*x/100).^2;
%(3)--------------------------
m1 = (99*mean(f1)+199*mean(f2)+999*mean(f3))/(99+199+999);
disp('由2产生的函数值的平均值')
disp(m1)
%(4)-----------------------
xx = fsolve(@(x) sin(m*x/100)^2-m1,0.1);
disp('对于x解方程f(x)=均值')
disp(xx)
追问
你这是复制来的,不一样,我看过了
已赞过
已踩过<
评论
收起
你对这个回答的评价是?
推荐律师服务:
若未解决您的问题,请您详细描述您的问题,通过百度律临进行免费专业咨询