在赋值 A(:) = B 中，A 和 B 中的元素数目必须相同。

purpose_857

用matlab求三元一次微分，出现了如下问题，请大神帮忙看下：
问题：在赋值 A(: ) = B 中，A 和 B 中的元素数目必须相同。

/ d; f3 ?  C/ @7 v: H7 R出错 longgekuta_lunwen2 (line 41)
    y(i+1) =y(i)+h/6*(K1+2*K2+2*K3+K4);
5 x3 v6 M  O+ G- i" h/ v+ `) C$ h- F
代码如下：clc;
( t( t+ C- W7 |9 D1 pclear all;
, ^( Y* N7 N# `; o" ?global r midu E Es u c q t w i n m1 m2 A B U1 S1 o j %o为弹塑性界面速度（C)与缓冲层弹性波速（Cd）的比值 f为应变值
3 O- k7 N0 u& D2 `; W8 S0 h* Qr=0.2; %落石半径
- p6 G' X. U4 Umidu=1640;%密度
E=35e6; %弹性模量
! D1 O0 m5 _! h2 y* rEs=8e6; %压缩模量
2 e3 D# l$ Y, D3 L" S* o1 ^u=0.32; %泊松比
; W0 X, s: D/ E! Q% fc=12e3; %粘聚力
q=30 ; %内摩擦角
t=sind(q);
( @3 r! ~( L" ^- y. g/ L2 @2 v+ `w=3*sind(q)/sqrt(3+sind(q)^2); %屈服强度
0 @4 i% N3 v! m" t0 a8 Ei=3*c*cosd(q)/sqrt(3+sind(q)^2);%剪切强度

n=3/(3+2*w);
4 e2 C3 s$ h; r7 A; p2 qm1=6*w/(3+2*w);
7 C+ g5 i6 H: B9 Xm2=6*i/(3+2*w*Es);
. V* t* J5 ~: J  C' O2 P! d: z# w
6 `8 S8 x( w# }6 t; Z% Q0 D" FA=6*(1+u)*(1-2*u)*o^2*i/((3+2*w)*(1-o)*E*((1-2*u)*(3-6*w/(3+2*w))*(1+o)-(6*w*o^2*(1+u))/(3+2*w)));
B=-3/2*A;
U1=3*(1-o^2)*A;
/ O$ I$ y, x) L  M4 Y, _5 m7 b  DS1=-E/(Es*(1+u)*(1-2*u))*((1+u)*A*o+2*(1-2*u)*B/(3*o^2)+2*u*B);
4 W/ z& k% o3 o2 v
h=j-1;
, v5 g* a: s5 e& w6 O4 wN =10;
1 ?( f. m6 D+ l8 I( B. K0 a) gx1=1;  %初始值1
y1=U1; %初始值1
! i* W$ ]! S4 E& d3 uz1=S1; %初始值1
/ L( F# x+ |* o4 ~4 M  yx2=j;  %初始值2
y2=j;  %初始值2
x = ones(1,N);
y = ones(1,N)*0;
z = ones(1,N)*0;
: Z% x3 f) B/ n3 _
* D0 ~1 L8 g$ R  C$ b1 f. tfor i=1:N
! t- i5 B' G  w    [K1,L1] = F(1,y(i),z(i))
    [K2,L2] = F((j+1)/2,y(i)+h*K1/2,z(i)+h*L1/2)
    [K3,L3] = F((j+1)/2,y(i)+h*K2/2,z(i)+h*L2/2)
    [K4,L4] = F(j,y(i)+h*K3,z(i)+h*L3)
2 V! w1 w$ z* Q; `6 p2 K' y, ~
# p" {* z) Y5 Z8 ?    y(i+1) =y(i)+h/6*(K1+2*K2+2*K3+K4);
  t8 B; ~8 f4 {2 {) q    z(i+1)= z(i)+h/6*(L1+2*L2+2*L3+L4);
; R. \4 _- a/ J5 t7 F% B4 x/ p. ?end
& w& [0 K. i5 ]' d: z, x
figure(1);
9 s  b7 J$ `- ~+ B4 U- T" rcomet3(x,y,z,0.1)
3 {( O  x4 |& u, ]! M# b+ cfunction [K,L] = F(X,U,S)
global r midu E Es u c q t w i n m1 m2 A B U1 S1 C   %C为弹塑性界面速度
K = (2*U/X+3*(m1*S+m2)*(X-U)/(X*(n-3*S+2*i/Es)))./((3*n*midu*C^2)*(X-U)^2/(Es*(n-3*S+2*i/Es)^2)-1);
8 |  _/ _. w) S, f2 }$ BL = ((2*U/X+3*(m1*S+m2)*(X-U)/(X*(n-3*S+2*i/Es)))./((3*n*midu*C^2)*(X-U)^2/(Es*(n-3*S+2*i/Es)^2)-1))*n*midu*C^2*(X-U)/(Es*(n-3*S+2*i/Es))-(m1*S+m2)/X;
end
3 W! ?( N& t3 s9 j请大神赐教，感谢!

peerless2021

只是单纯从代码角度出发，你这里的o, C需要赋值，否则就会导致空集，仅供参考

purpose_857

whatever_ 发表于 2023-3-1 11:27
仅供参考，主程序里的j，o，C三个变量未定义，且最后一个变量未在主程序里声明为全局变量

这里 o, C 是未知量，给 j 赋初始值的情况下，这里微分能求出来吗？
. s0 D; P7 v6 ]6 b  Y9 H" h

whatever_

仅供参考，主程序里的j，o，C三个变量未定义，且最后一个变量未在主程序里声明为全局变量