Question

matlab符号运算与赋值

如上图，在生成符号矩阵G的时候有警告，有什么更好的办法生成符号矩阵么？

如上图，fun是一个符号函数，然后x1,x2是自变量，对x1,x2求梯度。这里的[x1,x2]换成X（ = [x1,x2]），就会报错，这就使得变量个数改变的时候还要手动修改这里的[x1.x2,...].有没有办法不用这么麻烦的修改呢？

如上图，在对符号函数F进行赋值的时候，需要一个一个的进行赋值，调用多次subs函数，繁琐不说，F（i）这里还有警告，弄的挺烦人的，有没有好方法方便的对F中的各种符号赋值呢？
　　你好！我的程序是这样的：

　　主要的问题是：1、循环生成符号矩阵和向量，下面会有波浪线存在（第7、12、13行）；2、循环将符号进行赋值，红色波浪线存在（34、40行）；3、第38行，求梯度的时候 jacobian(fun,[x1,x2]) 这里的[x1,x2]如果用已经生成的向量符号X代替就会报错，想不明白为啥。。这样的结果是还得定义两个符号（第3行），并且变量个数变的时候还得手动改程序，比较麻烦；4、34、40行用了好几重subs，有没有简便点的方法，同时一次赋几个值。
问题比较多，帮忙看一下，小弟对matlab用法不太熟，求指导！分数会追加，谢谢！如果用得着，下面是复制过来的代码：





Answer 1

你这其实是三个问题，请把上述问题所涉及到的测试代码贴出来（如果受字数限制，可打包传到网盘）。不然为了回答你这个问题，还需要做上好半天无价值的文字录入工作。

当然，我也可以不做测试，只把思路大致说明一下，但你未必能看明白。

追问

你好！我把问题详细的写在了问题补充那里 从“你好！我的程序是这样的”开始。如果需要，后面我再把代码给你，谢谢！

追答

晕。。。你补充的是哪些我能看到。
关键就是要你把代码发给我，不然回答这个本来并不困难的问题，却要录入一大堆代码，不是浪费生命吗？

追问

已私信~

追答

算了，没有你那几个需要load的文件，择要回答吧：

 

问题1和2其实是相同的问题，就是数组没有提前分配存储空间，在循环过程中MATLAB会自动调整数组大小，这通常会影响效率，所以用波浪线提醒你。预置（preallocation）是用来解决这种问题的，例如：

G = sym(zeros(n));
for i =1:1:n %生成符号矩阵G
    for j = 1:1:n
        G(i,j) = sym (['g',num2str(i),num2str(j)]);
    end
end
X = sym(zeros(n,1));
Y = sym(zeros(n,1));
for i = 1:1:n %生成符号向量X,Y
    X(i) = sym(['x',num2str(i)]);
    Y(i) = sym(['y',num2str(i)]);
end
Z = sym('z');

不过，如果n不大，其实没必要预置。

 

问题3的原因在于，你对X、Y转置时，其实是“共轭转置”，导致

X = 
 conj(x1)
 conj(x2)

不是独立变量，当然会出错。正常的转置应该是类似X.'这样的（注意多了个小数点）。

 

对问题4的说明：

不知道你的具体数据形式，所以给点一般性说明：对于多个变量的代入，可以用类似subs(f,{x,y},{x0,y0})这样的形式；如果f的所有变量在工作区中都有具体值，也可以更简单地用subs(f)就可以了。需要注意的是，代入的所有变量都应该维数相同（或者都是标量，如果是数组，维数必需相同）。