189 8069 5689
可以调用 diff 函数求导。
目前累计服务客户上千多家，积累了丰富的产品开发及服务经验。以网站设计水平和技术实力，树立企业形象，为客户提供成都做网站、成都网站设计、网站策划、网页设计、网络营销、VI设计、网站改版、漏洞修补等服务。创新互联始终以务实、诚信为根本，不断创新和提高建站品质，通过对领先技术的掌握、对创意设计的研究、对客户形象的视觉传递、对应用系统的结合，为客户提供更好的一站式互联网解决方案，携手广大客户，共同发展进步。
举例说明：
先定义符号 x、y 以及符号二元函数表达式 z，然后调用 diff 函数求偏导，代码如下：
clc;clear;
syms x y
z=x^2+y^2+exp(x*y);
z_x=diff(z,x,1)
z_y=diff(z,y,1)
z_x2=diff(z,x,2)
z_y2=diff(z,y,2)
z_xy=diff(z,x,y)
结果如下：
z_x =
2*x + y*exp(x*y)
z_y =
2*y + x*exp(x*y)
z_x2 =
y^2*exp(x*y) + 2
z_y2 =
x^2*exp(x*y) + 2
z_xy =
exp(x*y) + x*y*exp(x*y)
图中调用了5次diff函数，分别计算了 z 对 x 的一阶偏导，z 对 y 的一阶偏导，z 对 x 的二阶偏导，z 对 y 的二阶偏导，以及 z 对 x,y 都求一次偏导。
例子中所用的调用格式为： diff（f，var，n）
f 为符号表达式，也可以是符号函数（这个数据类型低版本的matlab没有）。n 为求导次数，缺省为1。var 为求导的符号变量，可以缺省（matlab会根据表达式自己选择一个），但不建议缺省，除非表达式只含有一个符号变量。
此外，matlab还允许 diff（f，n，var）与 diff（S，v1，v2，...，vn）的调用形式。
diff（f，v1，v2，...，vn） 会把表达式 f 对变量 v1，v2 等 n 个变量都求一次偏导，得到 f 的 n 阶偏导。
matlab求偏导步骤如下。
1、第一步，双击打开matlab。
2、第二步，定义两个变量a、b,在窗口中输入代码：syms a b。
3、第三步，定义一个多元函数用于求偏导数，输入代码：f=exp(a)*exp(b)，即f=a^2*b^2。
4、第四步，分别对a、b求解偏导数，输入两行代码：dfa=diff(f,a)和dfb=diff(f,b)。
5、第五步，查看两个偏导数的图像，
输入代码：subplot(2,2,1:2),ezsurf(f),subplot(223),ezsurf(dfa),subplot(224),ezsurf(dfb)。
6、第六步，点击运行，选择本地磁盘的位置保存程序（注意文件名需要用因为开头且不能有中文）。
7、第七步，再次点击运行查看该函数的偏导求解图像如下：
题名:科学与工程数值算法-Java版
作者：丁军, 杨丽丽
出版社:清华大学出版社,2003
这本书挺好的，是去年我学《数值分析》的法宝。上面有微分方程的数值算法代码，推荐参考一下
如何用matlab进行多元函数偏导数计算
可以调用 diff 函数求导。
举例说明：
先定义符号 x、y 以及符号二元函数表达式 z，然后调用 diff 函数求偏导，代码如下：
clc;clear;
syms x y
z=x^2+y^2+exp(x*y);
z_x=diff(z,x,1)
z_y=diff(z,y,1)
z_x2=diff(z,x,2)
z_y2=diff(z,y,2)
z_xy=diff(z,x,y)
结果如下：
z_x =
2*x + y*exp(x*y)
z_y =
2*y + x*exp(x*y)
z_x2 =
y^2*exp(x*y) + 2
z_y2 =
x^2*exp(x*y) + 2
z_xy =
exp(x*y) + x*y*exp(x*y)
把
for i = 2:100
a(i)=1; %a2=a3=a4=...=a_num=1
end
改成
for i = 2:100
eval(['a',num2str(i) '=1;'])
end
你原来的写法是把数组a的第2到第100个元素（即a(2)-a(100)）赋值为1，而符号变量a2-a100仍然存在。