如何用Matlab程序编写等式两边同时定积分的问题

作者&投稿:寸的 (若有异议请与网页底部的电邮联系)
求MATLAB大神,算一个定积分的程序~

如果是写matlab语句的话,如下:
syms x
t = int(1/sqrt(25*(x^(1/3)-1)-23*(x-x^(1/3))),0,1)
不过,最终算不出结果,不是程序编写的问题,而是这个公式不是显示积分,因此是不能直接积分的~
你看看你的问题能不能有其他条件将计算公式转化一下,或许才能积分~

以f(x)=e^2x+sin(x+π/3) ,积分下限:a=0,积分上限:b=π/4 为例。
1.首先,建立被积函数M文件。
点击New,选择Function.

2.在编辑器中输入指令:
function f=f(x)
f=exp(2*x).*sin(x+pi/3);
(注意,用 .* )

3.保存函数文件到自己的工作路径,比如G:\matlab\work。
输入指令:
cd G:\matlab\work
是该路径成为当前路径。

4.调用积分函数squad进行计算。输入指令:
F=quad('f(x)',0,pi/4)
如图,回车后可得到计算结果1.8612。
其他函数,只要修改函数文件中的表达式即可。

function y=f(x1,y1,x2)
dx=(x2-x1)/100;
x=x1:dx:x2;
n=length(x);
for i=1:(n-1)
dy(i)=3/(x(i)+1)*dx;
end
y=y1+sum(dy);
%想法就是定积分的概念,将积分区间等分,取左端点值*dx为每份小区间的面积。

原式可变为:
dy/dx=1/3/(x+1)
y(x=0)=0;

由于表达式较简单,可以直接用dsolve
dsolve('Dy=1/3/(x+1)','y(0)=0','x')


matlab中B(k,k)是什么意思,下面是全程序。
楼主是要改错么?你的B定义成向量了,不是矩阵,所以B(i,k)会出错。是不是你定义错了,想这么定义:A=[1 2 1; 2 2 3; -1 -3 0];b=[0; 3; 2];若还有问题请追问

matlab求解方程组, x^2+y^2+z^2=100 a*x+b*y+c*z=10 其中
方程无解。可以很简单地画出两个方程的曲线来:ezplot('-0.050+0.015+0.070-2.*(d.*0.0215).^0.5-2.*(p.*0.0504).^0.5',[-0.02 0.02])hold onezplot('-0.015+d+p',[-0.02 0.02])由图可见,两条线无交点,也就意味着方程组无解。或者换另一种方式来分析:把方程1...

matlab编程……
\/\/int a = 1, b = 1, c = 1;\/\/for (;a<=100; a++)\/\/{ \/\/for (b=a;b<=100;b++)\/\/{ \/\/for (c=b;c<=100;c++)\/\/{ \/\/if (c*c == a*a+b*b)\/\/{ \/\/cout<<"a="<<a<<";\\tb="<<b<<";\\tc="<<c<<endl;\/\/} \/\/} \/\/} \/\/} \/\/return 0;\/\/} \/\/...

matlab求解矩阵方程组:A+B=C;A-B=D;其中,C、D已知,都是21*21的矩阵,请...
两式相加可得,A=(C+D)\/2 想减为:B=(C-D)\/2

matlab编程中b(【1 3】,2)=b([3 1],2)是不是错误的?
不是啊 比如说 b=[1,2,3;4,5,6;7,8,9]b([1 3],2)=b([3,1],2);disp(b)其实就是说把b的第二列的第三行元素和第一行元素对调,即2和8对调 b = 1 2 3 4 5 6 7 8 9 1 8 3 4 5 6 7 2 9 ...

同一个matlab程序 matlab 2013b运行出错 matlab2012a运行正确,谁能帮我...
超出了递归调用的次数。您可以试着增加允许的递归次数,例如:set(0,'RecursionLimit',1000)看能否运行。导致这个问题的可能原因:1、由于软件版本差异导致递归的出口条件得不到满足;2、程序中有随机性因素,不同次的运行结果可能不同。具体要看代码才能判断。

跪求大神解答matlab编程:例如a=1,b=2矩阵X[1]=[ 1 1 ;2 2] X[2]=...
>> Xa=[ 1 1 ;2 2];Xb=[2 2 ;3 3];Xab=Xa+Xb Xab = 3 3 5 5

matlab 编程问题
题1:基本运算 >> A=[22 33 4;1 3 31;54 2 32];B=[1 4 -3;3 2 3;-3 1 8];>> A+6*B ans =[28 57 -14; 19 15 49; 36 8 80]>> A+B*eye(3,3) %eye——单位矩阵 ans =[23 37 1; 4 5 34; 51 3 40]>> A*B...

关于乘除法数字谜的matlab程序ABCDEF*B=EFABCD,怎么写啊,我编写的实在...
076923*3=230769 可以用两种方法:1、笔算:把上式看成:ABCDEF*F=EFABCD,然后令F等于0、1、2、3...9,然后通过上算式的关系依次推出直至结果出来!2、编程序算:可以用turbo.c 2.0编写C语言程序,用FOR循环做,这需要一定的编程能力。结果如上,LZ可以验证一下。

求教matlab:已知x=1,y=2,先求解a=x+y,b=x*y,而后要求解t=a+b,这个程 ...
x = 1;y = 2;a = x+yb = x*yt = a+ba = 3 b = 2 t = 5

卢龙县19326007672: 求教如何用Matlab程序编写等式两边同时定积分的问题dx/(x+1)=3*dy x=(0,x) y=(0,y)完成等式两边积分后,得到一个关于x和y关系的等式. -
禾弯久强:[答案] 原式可变为: dy/dx=1/3/(x+1) y(x=0)=0; 由于表达式较简单,可以直接用dsolve dsolve('Dy=1/3/(x+1)','y(0)=0','x')

卢龙县19326007672: 如何用Matlab程序编写等式两边同时定积分的问题 -
禾弯久强: 原式可变为: dy/dx=1/3/(x+1) y(x=0)=0;由于表达式较简单,可以直接用dsolve dsolve('Dy=1/3/(x+1)','y(0)=0','x')

卢龙县19326007672: matlab如何拟合等式两边均有因变量y的等式中的参数. -
禾弯久强: 用matlab如何拟合等式两边均有因变量y的等式中的参数,可以这样来处理,将x、y看成x向量变量,即x=[x1 x2].(x1=x,x2=y),其自定义函数应这样来写,func=@(a,x)a(1)*x(:,1)*(1-a(2)*x(:,2))^a(2)/(1-a(2)*x(:,2)+x(:,2))^(a(2)-1)-x(:,2) 有了自定义函数,就可以可以按常规方法去拟合了.

卢龙县19326007672: matlab如何计算例如 x^2=x+5,等式两边都有相同的自变量,谢谢大家了 -
禾弯久强: 命令窗口输入:syms x;solve('x^2-x-5=0') ans = 1/2+1/2*21^(1/2) 1/2-1/2*21^(1/2) 具体的可以help solve.

卢龙县19326007672: 求MATLAB解方程编程 -
禾弯久强: 思路:对等式两边求导,变成二阶微分方程,然后求解.clc;clear v=dsolve('D2v=25/8*v')%v(0)=0f1=subs(v,'t',0)-0;f2=subs(v,'t','t0')-80;[C1,C2]=solve(f...

卢龙县19326007672: 多元等式方程怎么用matlab编程 -
禾弯久强: 你没给出m、k的具体值,所以我假设了,程序如下:m1=1;m2=2;m3=3;k1=0.1;k2=0.2;k3=0.3;f=@(x)(m1*k1*x(1)+m2*k2*x(2)+m3*k3*x(3))./(k1*x(1)+k2*x(2)+k3*x(3));x=[102132];y=f(x)

卢龙县19326007672: matlab求含有未知矩阵的等式 -
禾弯久强: 原式化为 A'PA-EPE+Q=0 => (A'-E)P(A-E)=-Q 然后两边左乘inv(A'-E),右乘inv(A-E) P就出来了, 至于Matlab代码很简单 自己试试

卢龙县19326007672: 现有两个符号等式,怎样用matlab求出p1和pn啊...求大神给出程序!(等式在追问里面)
禾弯久强: ep1 = '1-2*p1+c+(pn-c)*(y*k-y*q)/k*q/(k-q)=0'; ep2 = 't*(k-q-pn+(t*pn-y*k+y*q+y*p1*k-y*p1*q)/t/k*q)/(k-q)+(pn-c)*t*(-1+1/k*q)/(k-q)=0'; [p1,p2] = solve(ep1,ep2) p1 = -k*(-2*p1^2-1+3*p1-c+p1*c)/(-pn*k+c*k+2*pn^2-3*pn*c+c^2) p2 = k*(1-2*p1+c)/q/(-pn+c)

卢龙县19326007672: 如何用Matlab编写下列方程,请教高手指点! -
禾弯久强: 功能X1 X2] = myfunction3(A,B,C) δ= B * B-4 * A * C 如果三角洲> 0X1 =(B + SQRT(三角洲))/(2 *); χ2=(-β-的SQRT(δ))/(2 *);elseif的增量== 0X1 =-B /( 2); (χ2= X1); 其他 X1 =-b的/(2 *)+((开方(δ))/(2 *)); /> X2 = B /(2 *) - ((SQRT(三角洲))/(2 *)); 年底 您的程序有几个错误,所以这就像一个

卢龙县19326007672: 如何用MATLAB编写数学公式 -
禾弯久强: W=[1,2,3,4,5,6,7,8,9,10]; Y=[1,2,3,4,5,6,7,8,9,10] ;W2=W.^2 WY=W.*Y

本站内容来自于网友发表,不代表本站立场,仅表示其个人看法,不对其真实性、正确性、有效性作任何的担保
相关事宜请发邮件给我们
© 星空见康网