用matlab解下列方程组(含积分)

2021-08-14 05:14:29 字数 5091 阅读 3446

1楼:匿名用户

从原理上说,tx、ty应该可以是随时间变化的函数,但那样就成了一个泛函问题。如果tx、ty取常数值方程有解,而且没有其它要求,当然应该使用常数解。

现在假设tx、ty为常数,可以使用fsolve求解:

m=2400;

q=7500/2490;

a=1.63;

%tx=x1,ty=x2,t=x3

f=@(x)[quadl(@(t)x(1)./(m-q*t),0,x(3))-1692;

quadl(@(t)x(2)./(m-q*t)-a,0,x(3))-57;

sqrt(x(1)^2+x(2)^2)-7500];

x=fsolve(f,[5304,5304,700]);

fprintf('tx=%.2f\nty=%.2f\nt=%.2f\n',x)

得到的结果:

tx=6876.27

ty=2994.48

t=417.07

2楼:匿名用户

你给的问题还不是很清楚

初始条件,结束条件 也没有给清楚

你的两个微分方程应该是关于两个速度 vx vy的,他们的导数就是加速度,其实就是两个方向的牛顿方程dvx f vx---- = -------- ------ (1)

dt m-q*t v

dvy f vy---- = -------- ------ - a (2)

dt m-q*t v

其中 v = sqrt(vx^2+vy^2) 是速度总值然后就是要给出 边界上的条件和约束

例如vx(0) vy(0) 的值(初始两个方向的速度)最后时刻应该满足的条件

这些条件,参数你都要给齐 ,而且注意参数的单位,最后都用标准国际单位这样才能够正确的得到解

matlab求解含积分的方程组

3楼:匿名用户

用matlab求解含积分的方程组,可以这样实现:

1、先用integral函数,求解定积分,即

integral(@(t)sqrt(a^2.*sin(t).^2+b^2.*cos(t).^2),0,pi/2)

2、再用vpasolve函数,求解β、θ值,即

[beta1,theta1]=vpasolve(eq1,eq2, [beta1,theta1])

这里:eq1是指方程一,eq2是指方程二

3、详细的运行**为

a=15;b=14;r=14.5;

syms theta1 beta1

eq1=integral(@(t)sqrt(a^2.*sin(t).^2+b^2.*cos(t).^2),0,pi/2)+beta1*pi*r/180;

eq2=(a*sin(theta1))/(b*cos(theta1));

[beta1,theta1]=vpasolve(eq1==100,eq2==tan(beta1), [beta1,theta1])

这里:beta1是指β,theta1是指θ

运行结果

用matlab求一个含有积分的方程组 5

4楼:匿名用户

用fsolve求一个含有积分的方程组,可以这样来写:

第一步:建立自定义函数 func()

第二步:用fsolve求解其近似值

[k,fval] =fsolve(@(k)func(k),k0) %k为系数a、b,即a=k(1),b=k(2),fval是每个等式的值

第三步:编程运行

运行结果:a=9.01704582797954;b =25.897553850857

执行**及运行过程。

5楼:云南美食吧

(1830658703514901060,替换网页地址question后面的那串数字即可)。

这个问题应该很难求出解析解,应该使用fsolve求数值解。

参考**(题中所给条件u=1/6疑为u=g/6之误):

m = 2400;q = 7500/2490;g = 9.8;a = g/6;% tx = x1, ty = x2, t = x3f = @(x) [quadl(@(t)x(1)./(m-q*t),0,x(3))-1692; quadl(@(t)x(2).

/(m-q*t)-a,0,x(3))-57; sqrt(x(1)^2+x(2)^2)-7500];x = fsolve(f,[5304,5304,700]);fprintf('tx = %.2f\nty = %.2f\nt = %.

2f\n',x)

求出的结果:

tx = 6873.97ty = 2999.75t = 417.17

6楼:狮子中国博物馆

7(2x-1)-3(4x-1)=4(3x 2)-1; (5y 1) (1-y)= (9y 1) (1-3y); 20% (1-20%)(320-x)=320×40% 2(x-2) 2=x 1 2(x-2)-3(4x-1)=9(1-x) x/3 -5 = (5-x)/2 2(x 1) /3=5(x 1) /6 -1 (1/5)x 1 =(2x 1)/4 (5-2)/2 - (4 x)/3 =1 x/3 -1 = (1-x)/2 (x-2)/2 - (3x-2)/4 =-1 11x 64-2x=100-9x 15-(8-5x)=7x (4-3x) 3(x-7)-2[9-4(2-x)]=22 3/2[2/3(1/4x-1)-2]-x=2 2(x-2)-3(4x-1)=9(1-x) 11x 64-2x=100-9x 15-(8-5x)=7x (4-3x) 3(x-7)-2[9-4(2-x)]=22 3/2[2/3(1/4x-1)-2]-x=2 2(x-2) 2=x 1 7(2x-1)-3(4x-1)=4(3x 2)-1(5y 1) (1-y)= (9y 1) (1-3y)[ (- 2)-4 ]=x 220% (1-20%)(320-x)=320×40%2(x-2) 2=x 1 6。

2(x-2)-3(4x-1)=9(1-x) 7。11x 64-2x=100-9x 15-(8-5x)=7x (4-3x) 3(x-7)-2[9-4(2-x)]=22 3/2[2/3(1/4x-1)-2]-x=25x 1-2x=3x-23y-4=2y 187x*13=57z/93=41 15x 863-65x=54 58y*55=274892(x 2) 4=92(x 4)=103(x-5)=184x 8=2(x-1)3(x 3)=9 x6(x/2 1)=129(x 6)=632 x=2(x-1/2)8x 3(1-x)=-27 x-2(x-1)=1x/3 -5 = (5-x)/2 2(x 1) /3=5(x 1) /6 -1 (1/5)x 1 =(2x 1)/4 (5-2)/2 - (4 x)/3 =1 15x-8(5x 1。

怎么用matlab解含有字母系数的方程组的解,举个简单例子

7楼:大野瘦子

用法以这个为例:

x+a*y=10

x-b*y=1

其中x,y为变量,a,b为字母系数.

只要在matlab中输入

syms x,y,a,b

[x y]=solve('x+a*y=10','x-b*y=1','x','y')

即可求出解

x =(a + 10*b)/(a + b)

y =9/(a + b)

对于函数solve的具体用法,可以通过输入help solve来学习。

matlab中方程求解的基本命令

1.roots(p) %求多项式的根,其中p是多项式向量。

例求x3-x2+x-1=0的根

解:>>roots([1,-1,1,-1])

注: [1,-1,1,-1]在matlab中表示多项式 x3-x2+x-1

2.solve(fun) %求方程fun=0的符号解,如果不能求得精确的符号解,可以计算可变精度的数值解

例:用solve求方程x9+x8+1=0的根

解:>>solve(‘x^9+x^8+1’)

给出了方程的数值解(32位有效数字的符号量)

3.solve(fun,var) %对指定变量var求代数方程fun=0的符号解。

例:解方程 ax2+bx2+c=0

解:>>syms a b c x;

>>f=a*x^2+b*x+c;

>>solve(f)

如果不指明变量,系统默认为x,也可指定自变量,比如指定b为自变量

>>symsa b c x;

>> f=a*x^2+b*x+c;

>>solve(f,b)

4.fsolve(fun,x0) %求非线性方程fun=0在估计值x0附近的近似解。

例:用fsolve求方程x=e-x在0附近的根

解:>>fsolve(‘x-exp(-x)’,0)

5.fzero(fun,x0) %求函数fun在x0附近的零点

例:求方程x-10x+2=0在x0=0.5附近的根

解:>>fzero(‘x-10^x+2’,0.5)

8楼:匿名用户

举个简单例子,解方程组

x+a*y=10

x-b*y=1

其中x,y为变量,a,b为字母系数。

只要在matlab中输入

syms x,y,a,b

[x y]=solve('x+a*y=10','x-b*y=1','x','y')

即可求出解

x =(a + 10*b)/(a + b)y =9/(a + b)

对于函数solve的具体用法,可以通过输入help solve来学习。

希望我的回答能够解决您的疑问,谢谢。

9楼:匿名用户

例:解方程组a*x+2*y=4,4*x+b*y=7matlab**:syms a b x y;

z=solve('a*x+2*y=4','4*x+b*y=7')即可。

利用定积分的定义计算下列定积分,利用定积分定义计算下列积分

1楼 匿名用户 写成a 1,b 2也没错,但是此时函数f x 根号 x ,而不是根号 1 x 。你再好好看看。 利用定积分定义计算下列积分 2楼 社工制编组 这个题目很基础的,多看几遍书一定能做出来,要相信自己 3楼 匿名用户 s of pilot knobs th 利用定积分定义计算下列极限 4楼...

用定积分几何意义怎么求?求过程,利用定积分的几何意义,计算下列定积分

1楼 一切都是缘 定积分的几何意义是 在x轴上方的叫正面积 下方的叫负面积 利用定积分的几何意义,计算下列定积分 2楼 我不是他舅 y 9 x x y 9 且y 9 x 0 所以是圆在x轴上方的部分 所以是半圆 且积分限 3到3,所以是整个半圆 半径是3 所以原式 9 2 3楼 桥梁abc也懂生活 ...