回 帖 发 新 帖 刷新版面

主题:有关解非线性方程组得问题

syms x p k;
r=p/(1-k*cos(4*x))
r1=diff(r,x)%中间变量
r11=diff(r,x,2)%中间变量
ft=int(r+(2*7.5*r*(r1^2-r*r11)/sqrt((r1^2+r^2)^3)))/(r+2*7.5*r/sqrt(r1^2+r^2))
fs=int(sqrt(r^2+r1^2))
fs1=subs(fs,x,eps)
fs2=subs(fs,x,0.5*pi-0.000001)
equ1=fs2-fs1-pi*1.5*10/6
ft1=subs(ft,x,eps)
ft2=subs(ft,x,0.5*pi)
equ2=ft2-ft1-pi/3
solve(equ1,equ2,k,p)
最后方程组解不出来,怎么回事啊?

回复列表 (共7个回复)

沙发


你的f表达式有问题,我换成别的表达式,能够显示图形

板凳


不好意思,回错地方了

3 楼


没关系,能帮忙就好了,呵呵
我替那个们谢谢你

4 楼


把你的ft表达式后面 sqr改称sqrt,你写错了

5 楼

原来角文本里是sqrt,不过没用阿

6 楼

可以出结果,我的运行情况如下:
syms x p k;
r=p/(1-k*cos(4*x))
r1=diff(r,x)%中间变量
r11=diff(r,x,2)%中间变量
ft=int(r+(2*7.5*r*(r1^2-r*r11)/sqrt((r1^2+r^2)^3)))/(r+2*7.5*r/sqrt(r1^2+r^2))
fs=int(sqrt(r^2+r1^2))
fs1=subs(fs,x,eps)
fs2=subs(fs,x,0.5*pi-0.000001)
equ1=fs2-fs1-pi*1.5*10/6
ft1=subs(ft,x,eps)
ft2=subs(ft,x,0.5*pi)
equ2=ft2-ft1-pi/3
solve(equ1,equ2,k,p)
 
r =
 
p/(1-k*cos(4*x))
 
 
r1 =
 
-4*p/(1-k*cos(4*x))^2*k*sin(4*x)
 
 
r11 =
 
32*p/(1-k*cos(4*x))^3*k^2*sin(4*x)^2-16*p/(1-k*cos(4*x))^2*k*cos(4*x)
 

7 楼

LS 用的是MATLAB 吗 ?怎么可以是符号解
MATLAB 不是据说只能为数值解的吗

我来回复

您尚未登录,请登录后再回复。点此登录或注册