Matlab解决前循环方程式的系统更快

Question

我正在研究Matlab中的代码，它涉及到一个方程组的系统图形解决方案。现在，我知道这样做的唯一方法是这样的：

E = 200000; v = .3;  
S = [1/E, -v/E, -v/E; 
-v/E, 1/E, -v/E; 
-v/E, -v/E, 1/E]; 
e1 = linspace(0,eu,52); 
sig10 = zeros(size(e1)); 
syms e2; syms e3; syms sig1; 

a = 0; 
for k = 1:52 
    sig = (S^-1)*[e1(k); e2; e3]; 
    eq1 = sig(3,1) == 0; 
    eq2 = sig(2,1) == a*sig1; 
    eq3 = sig(1,1) == sig1; 
    [solsig, sol2, sol3] = solve([eq1, eq2, eq3], [sig1, e2, e3]); 
    sig10(1,k) = solsig; 
end 

还有其他几个这样的循环。它可以工作，但是需要永久执行。我试图象征性地解决循环外的系统问题，但我无法弄清楚如何在需要的地方放回k变量。

syms tous; syms sigs1; syms sigs2; 
sigxhill = zeros(1,361); 

for t=0:360 
    alpha = t; 
    eq1 = sigsx == sig1*cosd(alpha)^2+sigs2*sind(alpha)^2; 
    eq2 = 0 == sig1*sind(alpha)^2+sigs2*cosd(alpha)^2; 
    eq3 = tous == (sig1 - sigs2)*sind(alpha)*cosd(alpha); 
    eq4 = 2*sig0^2 == F*sigys^2+G*sigsx^2+H*(sigsx-sigsy)^2+2*N*tous^2; 
    sol1 = solve([eq1, eq2, eq3, eq4], [sigsx, tous, sigs1, sigs2]); 
    sigxhill(1,t+1) = sol1.sigsx(1); touhill(1,t+1) = sol1.tous(1); 
    sig1hill(1,t+1) = sol1.sigs1(1); sig2hill(1,t+1) = sol1.sigs2(1); 
end 

没有任何人有任何想法：此外，由于这是一个矩阵的问题，这里的另一个代码，竟然放弃了我类似的问题，因为我知道有可能是一个简单的方法来进行此事？

Answer 1

不幸的是我不能执行你的代码，因为我没有安装符号数学工具箱。但是，在放入系统情况下的好办法，开始是使用MATLAB探查：

profile on; 
E = 200000; v = .3;  
S = [1/E, -v/E, -v/E; 
-v/E, 1/E, -v/E; 
-v/E, -v/E, 1/E]; 
e1 = linspace(0,eu,52); 
sig10 = zeros(size(e1)); 
syms e2; syms e3; syms sig1; 

a = 0; 
for k = 1:52 
    sig = (S^-1)*[e1(k); e2; e3]; 
    eq1 = sig(3,1) == 0; 
    eq2 = sig(2,1) == a*sig1; 
    eq3 = sig(1,1) == sig1; 
    [solsig, sol2, sol3] = solve([eq1, eq2, eq3], [sig1, e2, e3]); 
    sig10(1,k) = solsig; 
end 
profile viewer; 

这应该给你的代码的部分，其中最需要的执行时间一些额外的信息。

也见一些详细信息，文档：

http://de.mathworks.com/help/matlab/ref/profile.html