i tried to find k which depends on w , but i cannot preallocate p3,p4.....,please help?

Question

Rosalinda 2023 年 10 月 18 日

0
リンク

この質問への直接リンク

https://jp.mathworks.com/matlabcentral/answers/2035419-i-tried-to-find-k-which-depends-on-w-but-i-cannot-preallocate-p3-p4-please-help

コメント済み: Rosalinda 2023 年 10 月 30 日

syms k 
%declaring variables
c44=43;
e15=-0.48;
e11=7.61;
u11=100;
d11=2.6;
p=5700;
pe=3512;
c44e=533.34;
e11e=5.02;
csh=0.38969499199;
q=1.602*(1/10^(19));
uo=10^(20);
i=sqrt(-1);
%for loop for w to k
w=0:0.2:3;
 n=length(w);
 z1=zeros(1,n);
s=sym(zeros(1,n));
s3=sym(zeros(1,n));
b=sym(zeros(1,n));
c=sym(zeros(1,n));
y=sym(zeros(1,n));
z=sym(zeros(1,n));
s4=sym(zeros(1,n));
s5=sym(zeros(1,n));
s6=sym(zeros(1,n));
p3=sym(zeros(1,n));
p4=sym(zeros(1,n));
p5=sym(zeros(1,n));
p6=sym(zeros(1,n));
q3=sym(zeros(1,n));
q4=sym(zeros(1,n));
q5=sym(zeros(1,n));
q6=sym(zeros(1,n));
 for j=1:n
 s(j)=sqrt(1-(w(j)/k*csh)^2);
a=d11*(c44*e11+e15^2);
b(j)=e11*d11*p*w(j)^2/k^2+w(j)*(c44*e11+e15^2)*i/k^2-c44*q*uo*u11/k^2;
c(j)=(e11*w(j)*i/k^2-q*uo*u11/k^2)*p*w(j)^2/k^2;
s3(j)=sqrt(1+(-b(j)-sqrt(b(j)^2-4*a*c(j)))/2*a);
s4(j)=-s3(j);
s5(j)=sqrt(1+(-b(j)+sqrt(b(j)^2-4*a*c(j)))/2*a);
s6(j)=-s5(j);
p3(j)=-e15*(s3(j)^2-1)/(c44*(s3(j)^2-1)+p*w(j)^2/k^2);
p4(j)=-e15*(s4(j)^2-1)/(c44*(s4(j)^2-1)+p*w(j)^2/k^2);
p5(j)=-e15*(s5(j)^2-1)/(c44*(s5(j)^2-1)+p*w(j)^2/k^2);
p6(j)=-e15*(s6(j)^2-1)/(c44*(s6(j)^2-1)+p*w(j)^2/k^2);
q3(j)=-uo*u11*(s3(j)^2-1)/(d11*(s3(j)^2-1)+w(j)*i/k^2);
q4(j)=-uo*u11*(s4(j)^2-1)/(d11*(s4(j)^2-1)+w(j)*i/k^2);
q5(j)=-uo*u11*(s5(j)^2-1)/(d11*(s5(j)^2-1)+w(j)*i/k^2);
q6(j)=-uo*u11*(s6(j)^2-1)/(d11*(s6(j)^2-1)+w(j)*i/k^2);
A=[e15*exp(k),-e15*exp(-k),(c44*p3(j)+e15)*s3*exp(s3(j)*k),(c44*p4(j)+e15)*s4*exp(s4(j)*k),(c44*p5(j)+e15)*s5(j)*exp(s5*k),(c44*p6(j)+e15)*s6*exp(s6(j)*k),0,0;
    e11*exp(k),e11*exp(-k),(e15*p3(j)-e11)*s3(j)*exp(s3(j)*k),(e15*p4(j)-e11)*s4(j)*exp(s4(j)*k),(e15*p5(j)-e11)*s5(j)*exp(s5(j)*k),(e15*p6(j)-e11)*s6(j)*exp(s6(j)*k),0,0;
    -q*uo*u11*exp(k),q*uo*u11*exp(-k),-q*exp(s3(j)*k)*(uo*u11*s3(j)-d11*q3(j)*s3(j)),-q*exp(s4(j)*k)*(uo*u11*s4(j)-d11*q4(j)*s4(j)),-q*exp(s5(j)*k)*(uo*u11*s5(j)-d11*q5(j)*s5(j)),-q*exp(s6(j)*k)*(uo*u11*s6(j)-d11*q6(j)*s6(j)),0,0;
    -q*uo*u11,q*uo*u11,-q*(uo*u11*s3(j)-d11*q3(j)*s3(j)),-q*(uo*u11*s4(j)-d11*q4(j)*s4(j)),-q*(uo*u11*s5(j)-d11*q5(j)*s5(j)),-q*(uo*u11*s6(j)-d11*q6(j)*s6(j)),0,0;
    e15,-e15,(c44*p3(j)+e15)*s3(j),(c44*p4(j)+e15)*s4(j),(c44*p5(j)+e15)*s5(j),(c44*p6(j)+e15)*s6(j),-c44e*s,0;
    0,0,p3(j),p4(j),p5(j),p6(j),-1,0;
    -e11,e11,(e15*p3(j)-e11)*s3(j),(e15*p4(j)-e11)*s4(j),(e15*p5(j)-e11)*s5(j),(e15*p6(j)-e11)*s6(j),0,e11e;
    1,1,1,1,1,1,0,-1];
y(j)=det(A);
z(j)=vpasolve(y(j),k);
z1(j)=real(z);
 end
Error using sym/cat>checkDimensions
CAT arguments dimensions not consistent.

Error in sym/cat>catMany (line 33)
[resz, ranges] = checkDimensions(sz,dim);

Error in sym/cat (line 25)
    ySym = catMany(dim, args);

Error in sym/vertcat (line 19)
    ySym = cat(1,args{:});
 

0 件のコメント
-2 件の古いコメントを表示-2 件の古いコメントを非表示

サインインしてコメントする。

サインインしてこの質問に回答する。

Answer 1

Torsten 2023 年 10 月 18 日

0
リンク

この回答への直接リンク

https://jp.mathworks.com/matlabcentral/answers/2035419-i-tried-to-find-k-which-depends-on-w-but-i-cannot-preallocate-p3-p4-please-help#answer_1336114

編集済み: Torsten 2023 年 10 月 18 日

MATLAB Online で開く

Use

A=[e15*exp(k),-e15*exp(-k),(c44*p3(j)+e15)*s3(j)*exp(s3(j)*k),(c44*p4(j)+e15)*s4(j)*exp(s4(j)*k),(c44*p5(j)+e15)*s5(j)*exp(s5(j)*k),(c44*p6(j)+e15)*s6(j)*exp(s6(j)*k),0,0;
    e11*exp(k),e11*exp(-k),(e15*p3(j)-e11)*s3(j)*exp(s3(j)*k),(e15*p4(j)-e11)*s4(j)*exp(s4(j)*k),(e15*p5(j)-e11)*s5(j)*exp(s5(j)*k),(e15*p6(j)-e11)*s6(j)*exp(s6(j)*k),0,0;
    -q*uo*u11*exp(k),q*uo*u11*exp(-k),-q*exp(s3(j)*k)*(uo*u11*s3(j)-d11*q3(j)*s3(j)),-q*exp(s4(j)*k)*(uo*u11*s4(j)-d11*q4(j)*s4(j)),-q*exp(s5(j)*k)*(uo*u11*s5(j)-d11*q5(j)*s5(j)),-q*exp(s6(j)*k)*(uo*u11*s6(j)-d11*q6(j)*s6(j)),0,0;
    -q*uo*u11,q*uo*u11,-q*(uo*u11*s3(j)-d11*q3(j)*s3(j)),-q*(uo*u11*s4(j)-d11*q4(j)*s4(j)),-q*(uo*u11*s5(j)-d11*q5(j)*s5(j)),-q*(uo*u11*s6(j)-d11*q6(j)*s6(j)),0,0;
    e15,-e15,(c44*p3(j)+e15)*s3(j),(c44*p4(j)+e15)*s4(j),(c44*p5(j)+e15)*s5(j),(c44*p6(j)+e15)*s6(j),-c44e*s(j),0;
    0,0,p3(j),p4(j),p5(j),p6(j),-1,0;
    -e11,e11,(e15*p3(j)-e11)*s3(j),(e15*p4(j)-e11)*s4(j),(e15*p5(j)-e11)*s5(j),(e15*p6(j)-e11)*s6(j),0,e11e;
    1,1,1,1,1,1,0,-1];
y(j)=det(A)
z(j)=vpasolve(y(j),k);
z1(j)=real(z(j));

instead of

A=[e15*exp(k),-e15*exp(-k),(c44*p3(j)+e15)*s3*exp(s3(j)*k),(c44*p4(j)+e15)*s4*exp(s4(j)*k),(c44*p5(j)+e15)*s5(j)*exp(s5*k),(c44*p6(j)+e15)*s6*exp(s6(j)*k),0,0;
    e11*exp(k),e11*exp(-k),(e15*p3(j)-e11)*s3(j)*exp(s3(j)*k),(e15*p4(j)-e11)*s4(j)*exp(s4(j)*k),(e15*p5(j)-e11)*s5(j)*exp(s5(j)*k),(e15*p6(j)-e11)*s6(j)*exp(s6(j)*k),0,0;
    -q*uo*u11*exp(k),q*uo*u11*exp(-k),-q*exp(s3(j)*k)*(uo*u11*s3(j)-d11*q3(j)*s3(j)),-q*exp(s4(j)*k)*(uo*u11*s4(j)-d11*q4(j)*s4(j)),-q*exp(s5(j)*k)*(uo*u11*s5(j)-d11*q5(j)*s5(j)),-q*exp(s6(j)*k)*(uo*u11*s6(j)-d11*q6(j)*s6(j)),0,0;
    -q*uo*u11,q*uo*u11,-q*(uo*u11*s3(j)-d11*q3(j)*s3(j)),-q*(uo*u11*s4(j)-d11*q4(j)*s4(j)),-q*(uo*u11*s5(j)-d11*q5(j)*s5(j)),-q*(uo*u11*s6(j)-d11*q6(j)*s6(j)),0,0;
    e15,-e15,(c44*p3(j)+e15)*s3(j),(c44*p4(j)+e15)*s4(j),(c44*p5(j)+e15)*s5(j),(c44*p6(j)+e15)*s6(j),-c44e*s,0;
    0,0,p3(j),p4(j),p5(j),p6(j),-1,0;
    -e11,e11,(e15*p3(j)-e11)*s3(j),(e15*p4(j)-e11)*s4(j),(e15*p5(j)-e11)*s5(j),(e15*p6(j)-e11)*s6(j),0,e11e;
    1,1,1,1,1,1,0,-1];
y(j)=det(A);
z(j)=vpasolve(y(j),k);
z1(j)=real(z);

1 件のコメント
-1 件の古いコメントを表示-1 件の古いコメントを非表示

Rosalinda 2023 年 10 月 30 日

MATLAB Online で開く

syms k 
%declaring variables
c44=43;
e15=-0.48;
e11=7.61;
u11=100;
d11=2.6;
p=5700;
pe=3512;
c44e=533.34;
e11e=5.02;
csh=0.38969499199;
q=1.602*(1/10^(19));
uo=10^(20);
i=sqrt(-1);
%for loop for w to k
w=0:0.2:3;
 n=length(w);
 z1=zeros(1,n);
s=sym(zeros(1,n));
s3=sym(zeros(1,n));
b=sym(zeros(1,n));
c=sym(zeros(1,n));
y=sym(zeros(1,n));
z=sym(zeros(1,n));
s4=sym(zeros(1,n));
s5=sym(zeros(1,n));
s6=sym(zeros(1,n));
p3=sym(zeros(1,n));
p4=sym(zeros(1,n));
p5=sym(zeros(1,n));
p6=sym(zeros(1,n));
q3=sym(zeros(1,n));
q4=sym(zeros(1,n));
q5=sym(zeros(1,n));
q6=sym(zeros(1,n));
 for j=1:n
 s(j)=sqrt(1-(w(j)/k*csh)^2);
a=d11*(c44*e11+e15^2);
b(j)=e11*d11*p*w(j)^2/k^2+w(j)*(c44*e11+e15^2)*i/k^2-c44*q*uo*u11/k^2;
c(j)=(e11*w(j)*i/k^2-q*uo*u11/k^2)*p*w(j)^2/k^2;
s3(j)=sqrt(1+(-b(j)-sqrt(b(j)^2-4*a*c(j)))/2*a);
s4(j)=-s3(j);
s5(j)=sqrt(1+(-b(j)+sqrt(b(j)^2-4*a*c(j)))/2*a);
s6(j)=-s5(j);
p3(j)=-e15*(s3(j)^2-1)/(c44*(s3(j)^2-1)+p*w(j)^2/k^2);
p4(j)=-e15*(s4(j)^2-1)/(c44*(s4(j)^2-1)+p*w(j)^2/k^2);
p5(j)=-e15*(s5(j)^2-1)/(c44*(s5(j)^2-1)+p*w(j)^2/k^2);
p6(j)=-e15*(s6(j)^2-1)/(c44*(s6(j)^2-1)+p*w(j)^2/k^2);
q3(j)=-uo*u11*(s3(j)^2-1)/(d11*(s3(j)^2-1)+w(j)*i/k^2);
q4(j)=-uo*u11*(s4(j)^2-1)/(d11*(s4(j)^2-1)+w(j)*i/k^2);
q5(j)=-uo*u11*(s5(j)^2-1)/(d11*(s5(j)^2-1)+w(j)*i/k^2);
q6(j)=-uo*u11*(s6(j)^2-1)/(d11*(s6(j)^2-1)+w(j)*i/k^2);
A=[e15*exp(k),-e15*exp(-k),(c44*p3(j)+e15)*s3*exp(s3(j)*k),(c44*p4(j)+e15)*s4*exp(s4(j)*k),(c44*p5(j)+e15)*s5(j)*exp(s5*k),(c44*p6(j)+e15)*s6*exp(s6(j)*k),0,0;
    e11*exp(k),e11*exp(-k),(e15*p3(j)-e11)*s3(j)*exp(s3(j)*k),(e15*p4(j)-e11)*s4(j)*exp(s4(j)*k),(e15*p5(j)-e11)*s5(j)*exp(s5(j)*k),(e15*p6(j)-e11)*s6(j)*exp(s6(j)*k),0,0;
    -q*uo*u11*exp(k),q*uo*u11*exp(-k),-q*exp(s3(j)*k)*(uo*u11*s3(j)-d11*q3(j)*s3(j)),-q*exp(s4(j)*k)*(uo*u11*s4(j)-d11*q4(j)*s4(j)),-q*exp(s5(j)*k)*(uo*u11*s5(j)-d11*q5(j)*s5(j)),-q*exp(s6(j)*k)*(uo*u11*s6(j)-d11*q6(j)*s6(j)),0,0;
    -q*uo*u11,q*uo*u11,-q*(uo*u11*s3(j)-d11*q3(j)*s3(j)),-q*(uo*u11*s4(j)-d11*q4(j)*s4(j)),-q*(uo*u11*s5(j)-d11*q5(j)*s5(j)),-q*(uo*u11*s6(j)-d11*q6(j)*s6(j)),0,0;
    e15,-e15,(c44*p3(j)+e15)*s3(j),(c44*p4(j)+e15)*s4(j),(c44*p5(j)+e15)*s5(j),(c44*p6(j)+e15)*s6(j),-c44e*s,0;
    0,0,p3(j),p4(j),p5(j),p6(j),-1,0;
    -e11,e11,(e15*p3(j)-e11)*s3(j),(e15*p4(j)-e11)*s4(j),(e15*p5(j)-e11)*s5(j),(e15*p6(j)-e11)*s6(j),0,e11e;
    1,1,1,1,1,1,0,-1];
y(j)=det(A);
z(j)=vpasolve(y(j),k);
z1(j)=real(z(j));
 end
Error using sym/cat>checkDimensions
CAT arguments dimensions not consistent.

Error in sym/cat>catMany (line 33)
[resz, ranges] = checkDimensions(sz,dim);

Error in sym/cat (line 25)
    ySym = catMany(dim, args);

Error in sym/vertcat (line 19)
    ySym = cat(1,args{:});