高手帮忙，2个关于xy，并含修正贝塞尔函数的方程，怎样用matlab求解xy，谢谢~~~

第二类修正贝塞尔函数的零阶和一阶，分别怎样用matlab计算~

1、首先双击matlab软件图标，打开matlab软件，可以看到matlab软件的操作窗口。

2、使用syms命令，创建七个符号变量a、b、c、d、x、y、z。

3、使用符号变量c，创建函数A，其中A=sin(c)。

4、使用函数diff(A)，求解函数A的一阶微分。

5、使用符号变量d，创建函数B，其中 B=13*d^6。

6、使用函数diff(B)，求解函数B的一阶微分。

7、使用符号变量z，创建函数C，其中C=z^2*sin(z)。

请问你的题目抄全了吗，方程应该含有等号啊，可以运用隐函数求导的方法把y1、y2表示出来。

其实，这个问题本身并不困难，要点只在于怎样计算贝塞尔函数。
MATLAB提供了计算贝塞尔函数的函数，详情请参见另一个问题的回答。

 

首先，试图使用符号数学工具箱求解析解，代码如下：

syms x y
eq1=(cos(x)-1/3*(cos(x))^3)+2*(1000+683-2*2480)/ ...
    (3*(1000-683))-y*(sin(x))^2/0.000115-2*sin(x)* ...
    sin(x-0.959931)/0.000115^2/(9.8*(1000-683)/0.071);
eq2=sin(0.959931-x)*besselk(0,(9.8*(1000-683)/0.071)^ ...
    0.5*0.000115*sin(x))+(9.8*(1000-683)/0.071)^0.5*y* ...
    besselk(1,(9.8*(1000-683)/0.071)^0.5*0.000115*sin(x));
[x,y]=solve(eq1,eq2)

得到的结果为：

x =
7.2408611253328779611417644360562
y =
-.85870621389208280165235516006804e-6

这是因为无法求得解析解，因而调用了数值方法求解得到的结果（注意，在不同MATLAB版本中的处理方式可能存在差别，我使用Maple内核的2008a求解得到上述结果，但使用MuPad内核的2012b则得到复数解）。

 

由于这里求得的x不符合0-pi区间的要求，所以，考虑直接用数值方法求解方程（使用优化工具箱fsolve函数）。代码如下：

function zd

x0 = [pi/2; -1e-5];
options=optimset('Display','iter');
x = fsolve(@eqs, x0, options);
fprintf('x = %.6g, y = %.6g
', x);
fprintf('Eq1(x, y) = %.6g, Eq1(x, y) = %.6g
', eqs(x));

function f = eqs(X)
x = X(1);
y = X(2);
f(1)=(cos(x)-1/3*(cos(x))^3)+2*(1000+683-2*2480)/ ...
    (3*(1000-683))-y*(sin(x))^2/0.000115-2*sin(x)* ...
    sin(x-0.959931)/0.000115^2/(9.8*(1000-683)/0.071);
f(2)=sin(0.959931-x)*besselk(0,(9.8*(1000-683)/0.071)^ ...
    0.5*0.000115*sin(x))+(9.8*(1000-683)/0.071)^0.5*y* ...
    besselk(1,(9.8*(1000-683)/0.071)^0.5*0.000115*sin(x));

得到的结果：

x = 0.957676, y = -8.58706e-007
Eq1(x, y) = -5.15143e-014, Eq1(x, y) = -1.38778e-017

我觉得没有太多需要的说明的，所以就不多写了，有问题请追问。

 

另外，我注意到，前面用符号数学工具箱求出的解与使用fsolve得到的结果相差2*pi，有兴趣可自行验证。

 

==================================

由于百度知道系统抽风，另一个问题的答案“正在提交中”，我把主要内容贴到这里，供参考。

 

MATLAB提供了计算贝塞尔函数的函数，具体包括：

• besselj - 第一类贝塞尔函数，或简称贝塞尔函数；

• bessely - 第二类贝塞尔函数，又称诺伊曼函数（Neumann function）；

• besseli - 第一类修正贝塞尔函数；

• besselk - 第二类修正贝塞尔函数；

• besselh - 第三类贝塞尔函数，又称汉克尔函数（Hankel function）。

 

这几个函数的调用语法基本相同，例如

J = besselj(nu,Z)
J = besselj(nu,Z,1)
[J,ierr] = besselj(nu,Z)

其中，nu为贝塞尔函数的阶数，Z为函数自变量。阶数必须为实数，但Z可以是复数。

 

就你的问题而言，非常简单，K0(x)、K1(x)在MATLAB中的表达式分别为besselk(0,x)、besselk(1,x)。

另外值得一提的是，上述函数是MATLAB基本模块提供的特殊函数（也就是说不需要任何附加工具箱），采用数值方法计算；而符号数学工具箱则提供了第一和第二类的4个贝塞尔函数，名称和调用方式都与基本模块的4个函数完全一致，但支持微分、积分等符号运算。

购房人给别人发给如风

---

长武县15244022668： 请写出一个含有字母xy且包含加、减、乘、除、乘方的代数式 - ？
宿炒格列： 3xy(2)+7xy(2)3xy(2)-5xy÷xy.(2)表示平方或乘方,乘号因有x省略掉,希望能帮到你

长武县15244022668： 请高手帮我设计一条方程式,里头只有两个未知数XY,当X取任何数时,Y都是等于2,要复杂一点点的! - ？
宿炒格列：[答案] Y= ( X^6 + 6*X^5 + 15*X^4 + 20*X^3 + 15*X^2 + 6*X + 1 - (X+1)^6 )^2011 + 2 其中A^n形式代表A的n次幂

长武县15244022668： 写出两个含X y的二次多项式,使它们的和为xy,则这两个多项式可以是( ) ( )哪两个? - ？
宿炒格列： -xy+1 与2xy-1

长武县15244022668： 关于xy的多项式7x^2 - 3xy+6x^2+mxy+6+2xy+4x^2+5m中不含xy项,求这个多项式的常数项 - ？
宿炒格列： 依照题意,就是说-3xy+mxy+2xy=(-3+m+2)xy是不存在的,那么就是说-3+m+2=0,即m=1.再看常数项是6+5m=6+5=11.所以常数项为11.

长武县15244022668： 试至少写两个只含有字母x、y的多项式,且满足下列条件: - ？
宿炒格列： yx^5-xy+xy^2 yx^5-xy^2+xy yx^5-x^2y+xy^2

长武县15244022668： 口袋妖怪xy两个队伍大家帮忙选修改一下 - ？
宿炒格列： 整体思路差不多,做一些怪物的修改就行了.1. 放墙手选气球雷皇或者命玉雷皇,速度、耐久、输出要远远好于超能猫,而且气球雷皇不怕众多速度围巾的地震物攻手,确保第一回合能放出墙来,有了觉冰之后还能反杀,还能搞定火鸡的天敌家门鸟.2. mega手建议换一个,既然是火鸡接力接力的终点最好是mega袋龙,保证打击面.3. 地龙建议带速度围巾或者抗冰果.

长武县15244022668： 已知关于x、y的多项式 - ？
宿炒格列： x的平方-1/3kxy-3y的平方+3xy-9中不含xy的项.试求k的值.急啊,求各位大神帮忙!!! ∴-k/3+3=0; k=9;您好,很高兴为您解答,skyhunter002为您答疑解惑 如果本题有什么不明白可以追问,如果满意记得采纳 如果有其他问题请采纳本题后另发点击向我求助,答题不易,请谅解,谢谢. 祝学习进步

长武县15244022668： 请教初一数学题 关于x y的二次多项式最少有 项 最多有 项 - ？
宿炒格列： 根据题意可知此多项式中不含同类项,所以它不会有无数个项. 因为此题是关于XY的二次多项式,所以,在这个多项式中只会有二次项、一次项和常数项. 其中二次项有X^2、Y^2和XY三项,一次项有X和Y两项,常数项有一项,共计有6项. 因此,此多项式最多有6项.

长武县15244022668： 关于44+XXY染色体异常发生在???不懂啊,老师讲了还是不懂...生物高手帮帮忙~~~ - ？
宿炒格列： 44+XXY异常发生的原因如下:(1)该个体的母亲产生的异常卵子22+XX(正常卵子为22+X):这种情况又分为:A.卵原细胞在减一分裂的后期,两个X同时进入一个次级卵母细胞,另一个次级卵母细胞不含X染色体,含有两个X染色体的次级...

长武县15244022668： 《口袋妖怪XY》如何修改出闪光和幻之神奇宝贝? - ？
宿炒格列： 你好,修改需要去购买一个powersave修改器,大概200元左右,操作似乎也挺复杂的……