求matlab求定积分的问题~%积分函数function f=jifen(x)a=6378245;b=6356863;e=sqrt((a^2+b^2)/(a^2));f=a^2*(1-e^2)*(cos(x))/((1-((e^2)*(sin(x))^2))^2);%积分函数完在命令窗口输入[S,n]=quad('jifen',0,54*pi/180)后回车出现错误信

来源:学生作业帮助网 编辑:六六作业网 时间:2024/11/16 01:32:13
求matlab求定积分的问题~%积分函数functionf=jifen(x)a=6378245;b=6356863;e=sqrt((a^2+b^2)/(a^2));f=a^2*(1-e^2)*(cos

求matlab求定积分的问题~%积分函数function f=jifen(x)a=6378245;b=6356863;e=sqrt((a^2+b^2)/(a^2));f=a^2*(1-e^2)*(cos(x))/((1-((e^2)*(sin(x))^2))^2);%积分函数完在命令窗口输入[S,n]=quad('jifen',0,54*pi/180)后回车出现错误信
求matlab求定积分的问题~
%积分函数
function f=jifen(x)
a=6378245;
b=6356863;
e=sqrt((a^2+b^2)/(a^2));
f=a^2*(1-e^2)*(cos(x))/((1-((e^2)*(sin(x))^2))^2);
%积分函数完
在命令窗口输入
[S,n]=quad('jifen',0,54*pi/180)
后回车
出现错误信息:
Error using ==> mpower
Matrix must be square.
Error in ==> jifen at 5
f=a^2*(1-e^2)*(cos(x))/((1-((e^2)*(sin(x))^2))^2);
Error in ==> quad at 62
y = f(x,varargin{:});
%错误信息完毕
本文是matlab初学者,只想通过matlab运算得出一个数值.请各位达人赐教!
另外,如果积分的上限是 通过有规律的 采样得到的数值,
希望得出一系列的 积分值S,那么如何引用S的值呢?

求matlab求定积分的问题~%积分函数function f=jifen(x)a=6378245;b=6356863;e=sqrt((a^2+b^2)/(a^2));f=a^2*(1-e^2)*(cos(x))/((1-((e^2)*(sin(x))^2))^2);%积分函数完在命令窗口输入[S,n]=quad('jifen',0,54*pi/180)后回车出现错误信

f=a^2*(1-e^2)*(cos(x))/((1-((e^2)*(sin(x))^2))^2);
改为
f=a^2*(1-e^2)*(cos(x))./((1-((e^2)*(sin(x)).^2)).^2);
即可.
结果:
S =
-2.4497e+015
n =
10017