如何在parfor循环中使用之前定义过的符号变量而不会出错

shuddkk

for i=1:2*n+2                                                              %生成2N+2阶的矩阵
/ P, k% o6 f. Y3 {9 s! S8 u! k    if i==1
. ^2 b4 t, ~1 g. p! I        f(i)=sigmar1(i)+P0;                                                 %第1层公式
$ S# |' D4 P* V0 C" p# Y        paRFor j=1:2*n+2
3 c, Q% Y: Q& M- F, ?) F3 f
5 {& I6 b5 ?0 @2 X8 K            if j<=n
0 F6 B8 a% D  P                KK(i,j)=equationsToMatrix(f(i),['D_',num2str(j)]);          %将线性方程组写成矩阵的形式
            elseif j<=2*n
                KK(i,j)=equationsToMatrix(f(i),['E_',num2str(j-n)]);
            elseif j==2*n+1
                KK(i,j)=equationsToMatrix(f(i),epsilon0);
& z* Z- q/ h# U3 [# @            else
                KK(i,j)=equationsToMatrix(f(i),gama0);
, }  |0 @6 e/ b/ B  M) |/ d$ g            end
5 P1 m; I2 H$ }$ I' A        end
    elseif i<=n
/ Y4 N( o3 ]2 y2 R7 B( r        f(i)=Ur1(i-1)-Ur2(i);                                             %第2层到n层公式
8 ], P8 N0 y* R' u; o3 L3 W0 y        for j=1:2*n+2
            if j<=n
' v7 [% @# i  E( ]6 ^0 \+ A- w                KK(i,j)=equationsToMatrix(f(i),['D_',num2str(j)]);
            elseif j<=2*n
                KK(i,j)=equationsToMatrix(f(i),['E_',num2str(j-n)]);
: F4 o$ W% E( U/ x/ t1 G            elseif j==2*n+1
                KK(i,j)=equationsToMatrix(f(i),epsilon0);
            else
/ X: [7 d' k7 }; }$ l$ L. U8 [                KK(i,j)=equationsToMatrix(f(i),gama0);
5 U! ]! s6 L* U' ^: o8 m            end
        end
- R; s. z% b/ u  H* l1 C    elseif i<=2*n-1
        f(i)=sigmar1(i-n)-sigmar2(i-n+1);                                 %第N+1层到2N-1层公式
        for j=1:2*n+2
            if j<=n
0 |" @  Z0 f" ]0 s5 _                KK(i,j)=equationsToMatrix(f(i),['D_',num2str(j)]);          %生成2N+2阶的矩阵
% V$ O: M: L  N            elseif j<=2*n
; {! j# `! o( s- T9 b3 f& h. r4 `                KK(i,j)=equationsToMatrix(f(i),['E_',num2str(j-n)]);
            elseif j==2*n+1
8 @) P9 u& I7 R. d: G                KK(i,j)=equationsToMatrix(f(i),epsilon0);
            else
                KK(i,j)=equationsToMatrix(f(i),gama0);
            end
: l5 F! e4 ]  F        end
    elseif i==2*n
        f(i)=sigmar1(n);                                                   %第2*N层公式
; T7 R+ k" x7 K& W        for j=1:2*n+2
0 R2 ?. O8 H( j% S$ s4 Z- U2 t* ?            if j<=n
                KK(i,j)=equationsToMatrix(f(i),['D_',num2str(j)]);
            elseif j<=2*n
2 K' L" ]' i2 L% Z0 {                KK(i,j)=equationsToMatrix(f(i),['E_',num2str(j-n)]);
, C! U0 `* ~8 `& m$ H" z9 D            elseif j==2*n+1
                KK(i,j)=equationsToMatrix(f(i),epsilon0);
7 {  g" }5 v- C+ R3 l/ {$ b' A            else
                KK(i,j)=equationsToMatrix(f(i),gama0);
7 y& Q7 _$ L$ q( ~( l            end
        end
4 ?. p/ k& x- `: a    elseif i==2*n+1
7 B8 L* P3 J6 z( L: j: I( V) n0 H$ f        for j=1:2*n+2
$ q! [( Z7 p# S            if j<=n
                KK(i,j)=xxD(j);
            elseif j<=2*n
                KK(i,j)=xxE(j-n);
3 |& j5 Y, W, v. @            elseif j==2*n+1
: i1 c& S: n5 ]3 @% d                KK(i,j)=tempepsilonx;
            else
                KK(i,j)=tempgamax;
            end
        end
# ~; D! f- X6 t/ X( X( ~    elseif i==2*n+2
        for j=1:2*n+2
. f9 a, k" F; i" l: i  ]& I            if j<=n
                KK(i,j)=yyD(j);
            elseif j<=2*n
5 f3 [4 p& B! s8 @, m' C                KK(i,j)=yyE(j-n);
% f# n4 e% o  q3 f: }            elseif j==2*n+1
                KK(i,j)=tempepsilony;
- j, b' k8 h% m* n5 T            else
                KK(i,j)=tempgamay;
: H! u) `5 k! @5 v            end
: X3 Y  N" k7 o& N, Q, O        end
    end
end
3 o. y6 B2 g' m& L6 \想把其中的一部分for循环改成parfor循环，但是在运行时会出现错误
Analyzing and transferring files to the workers ...done.
# U  E  u. r, x& o错误使用 Matrix_Of_Coefficient_K (line 195)
在 'gama0' 的工作进程上引发了 UndefinedFunction 错误。这可能是因为在这
: q( g6 C5 T# _7 |些工作进程上无法访问包含 'gama0' 的文件。使用 addAttachedFiles(pool,
* y6 v3 a) V6 J9 n, t. ~files) 指定要附加的必要文件。有关详细信息，请参阅
$ L& {% B# [; Z# z7 h9 k# z9 t: s"parallel.Pool/addAttachedFiles" 的文档。
7 V) ]9 f  M! X1 p
原因:
    未定义函数或变量 'gama0'。

kekek

帮你顶一下

llbnmo

我的解决办法是运算部分 弄成函数   函数里面需要用的是 数据用 函数输入方式传递，这样就避免很多的 par 运行过程中的 参数隔离检查报错，而且弄成函数后 循环体看着也简单明了。

小小鲁班

来学习一下