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

purpose_857

用matlab求三元一次微分，出现了如下问题，请大神帮忙看下：
问题：在赋值 A(: ) = B 中，A 和 B 中的元素数目必须相同。
- c2 G5 B0 l9 [, i( A; P% k7 R2 b/ |
出错 longgekuta_lunwen2 (line 41)
    y(i+1) =y(i)+h/6*(K1+2*K2+2*K3+K4);
( C$ {7 |' T! s* A) y
- Y+ w" q+ S! R1 u$ k代码如下：clc;
clear all;
( ?: T. G; [! ]7 qglobal r midu E Es u c q t w i n m1 m2 A B U1 S1 o j %o为弹塑性界面速度（C)与缓冲层弹性波速（Cd）的比值 f为应变值
& A5 s9 a8 E7 }" D) rr=0.2; %落石半径
" e, ~. Q! }7 W- w2 O3 n- B- T! smidu=1640;%密度
: F9 Q$ a( s% w6 U; R& q- `E=35e6; %弹性模量
: U. S" A& `- \# Z: r5 A) iEs=8e6; %压缩模量
u=0.32; %泊松比
c=12e3; %粘聚力
q=30 ; %内摩擦角
# R* U- \7 r/ `" U! l6 k3 s, W9 @' G' bt=sind(q);
w=3*sind(q)/sqrt(3+sind(q)^2); %屈服强度
i=3*c*cosd(q)/sqrt(3+sind(q)^2);%剪切强度
* W. Y/ v; l1 B8 d1 T) v
n=3/(3+2*w);
m1=6*w/(3+2*w);
5 S5 E6 H6 y8 Z& E6 h8 N, rm2=6*i/(3+2*w*Es);
6 r- \, ^) v8 c4 E' h9 l4 o9 s1 _
/ N& o9 ]- D% G. _* jA=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)));
% ~$ f1 i" J/ @6 b5 D3 i1 lB=-3/2*A;
- [! ?' |! t' fU1=3*(1-o^2)*A;
- t) Y0 }- v6 ]  F5 {S1=-E/(Es*(1+u)*(1-2*u))*((1+u)*A*o+2*(1-2*u)*B/(3*o^2)+2*u*B);

h=j-1;
N =10;
x1=1;  %初始值1
) k) B" d, Z) h$ |2 Y# J0 m/ L; A1 Dy1=U1; %初始值1
7 ?- U0 g  f5 z5 zz1=S1; %初始值1
x2=j;  %初始值2
3 j9 d# t4 b8 x0 k* q9 P- c9 Uy2=j;  %初始值2
$ |2 ?' c  p8 Q* u" G9 [x = ones(1,N);
1 |) s9 p2 V# Q. M( O; xy = ones(1,N)*0;
z = ones(1,N)*0;

for i=1:N
0 K  R* ?& K5 H& t  c1 b% c    [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)
- u; S& p% _( _8 A
& E" \* s: @& I) t4 E    y(i+1) =y(i)+h/6*(K1+2*K2+2*K3+K4);
& R  L" N3 L- ~5 z6 S5 E    z(i+1)= z(i)+h/6*(L1+2*L2+2*L3+L4);
% Y" ?5 J1 S7 ~. m+ x& }% ], Yend

+ \' U) @/ n6 F0 n2 Pfigure(1);
comet3(x,y,z,0.1)
1 x6 ?7 ?- [  p% R1 c0 bfunction [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);
L = ((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
; ~# u( e* W7 c请大神赐教，感谢!

whatever_

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

purpose_857

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

这里 o, C 是未知量，给 j 赋初始值的情况下，这里微分能求出来吗？

peerless2021

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