matlab 解微分方程出错clear all;clc;load 15.mat;t=celiang(:,1);u1=celiang(:,2);syms u2 i5 ii x u5%定义u1,L1,L2,R1,C1,C3,ii 等vU1=celiang(:,2);k=length(u1);C1=1;h=5e-7;C3=1;R1=1;L1=1;R2=1;L2=1;C2=1;i1(1)=0;i1(2)=0;%在开始的两个点用的

来源：学生作业帮助网 编辑：作业帮 时间：2024/04/29 21:16:26

matlab 解微分方程出错clear all;clc;load 15.mat;t=celiang(:,1);u1=celiang(:,2);syms u2 i5 ii x u5%定义u1,L1,L2,R1,C1,C3,ii 等vU1=celiang(:,2);k=length(u1);C1=1;h=5e-7;C3=1;R1=1;L1=1;R2=1;L2=1;C2=1;i1(1)=0;i1(2)=0;%在开始的两个点用的
matlab 解微分方程出错
clear all;clc;
load 15.mat;
t=celiang(:,1);
u1=celiang(:,2);
syms u2 i5 ii x u5%定义u1,L1,L2,R1,C1,C3,ii 等
vU1=celiang(:,2);
k=length(u1);
C1=1;
h=5e-7;
C3=1;
R1=1;
L1=1;
R2=1;
L2=1;
C2=1;


i1(1)=0;i1(2)=0;%在开始的两个点用的是最简单的使之为零的方法.
for j=3:k-1
    i1(j)=(vU1(j+1)-vU1(j))*12*C1/(23*h)+16*i1(j-1)/23-5*i1(j-2)/23;
end
i1(k)=(vU1(k)-vU1(k-1))*C1/h;
i1=i1';


i3=i5+ii-i1;
fai=x;
p1 =  1.367e-014;       
p2 =-2.585e-027;       
p3 =  -8.37e-011;       
p4 =  3.585e-023;       
p5 =    3.2e-007;       
p6 =-1.298e-019;       
p7 =  -0.0002556;       
p8 =  1.143e-016;
ih= p1*x.^7 + p2*x.^6 + p3*x.^5 + p4*x.^4 + p5*x.^3 +p6*x.^2 + p7*x + p8;
i7=i3-ih;
i4=-100*i7;
u4=u5/100;
i2=-i4-i5;
%%%%%%%%%%%%%%%%%%%%%
s=dsolve(diff((u2-u1),t)==i5/C3,diff(i3,t)==(u1-R1*i3-u5)/L1,diff(fai,t)==u5,diff(i4,t)==(u2-u4- R2*i4)/L2);
u2=s.u2

已知u1如何求u2
i1=C1*du1/dt;%%含有微分
i5=C3*d(u2-u1)/dt;%%含有微分
i3=i5+ii-i1;
u5=u1-R1*i3-L1*di3/dt;%%含有微分
u5=dfai/dt;%%含有微分
fai=x;
p1 = 1.367e-014;
p2 = -2.585e-027;
p3 = -8.37e-011;
p4 = 3.585e-023;
p5 = 3.2e-007;
p6 = -1.298e-019;
p7 = -0.0002556;
p8 = 1.143e-016;
ih= p1*x.^7 + p2*x.^6 + p3*x.^5 + p4*x.^4 + p5*x.^3 + p6*x.^2 + p7*x + p8;
i7=i3-ih;
i4=-100*i7;
u4=u5/100;
u2=u4+R2*i4+L2*di4/dt;%%含有微分
i2+i4+i5=0;

matlab 解微分方程出错clear all;clc;load 15.mat;t=celiang(:,1);u1=celiang(:,2);syms u2 i5 ii x u5%定义u1,L1,L2,R1,C1,C3,ii 等vU1=celiang(:,2);k=length(u1);C1=1;h=5e-7;C3=1;R1=1;L1=1;R2=1;L2=1;C2=1;i1(1)=0;i1(2)=0;%在开始的两个点用的
clear all;clc;
load 15.mat;
t=celiang(:,1);
u1=celiang(:,2);
syms u2 i5 ii x u5%定义u1,L1,L2,R1,C1,C3,ii 等
vU1=celiang(:,2);
k=length(u1);
C1=1;
h=5e-7;
C3=1;
R1=1;
L1=1;
R2=1;
L2=1;
C2=1;
i1(1)=0;i1(2)=0;%在开始的两个点用的是最简单的使之为零的方法.
for j=3:k-1
i1(j)=(vU1(j+1)-vU1(j))*12*C1/(23*h)+16*i1(j-1)/23-5*i1(j-2)/23;
end
i1(k)=(vU1(k)-vU1(k-1))*C1/h;
i1=i1';
i3=i5+ii-i1;
fai=x;
p1 = 1.367e-014;
p2 =-2.585e-027;
p3 = -8.37e-011;
p4 = 3.585e-023;
p5 = 3.2e-007;
p6 =-1.298e-019;
p7 = -0.0002556;
p8 = 1.143e-016;
ih= p1*x.^7 + p2*x.^6 + p3*x.^5 + p4*x.^4 + p5*x.^3 +p6*x.^2 + p7*x + p8;
i7=i3-ih;
i4=-100*i7;
u4=u5/100;
i2=-i4-i5;
u21=u2-u1;
%%%%%%%%%%%%%%%%%%%%%
s=dsolve(diff(u21,t)==i5/C3,diff(i3,t)==(u1-R1*i3-u5)/L1,diff(fai,t)==u5,diff(i4,t)==(u2-u4- R2*i4)/L2);
u2=s.u2

matlab怎么解微分方程 matlab解偏微分方程, matlab怎么解微分方程? Matlab 出错 matlab出错, matlab 解微分方程出错clear all;clc;load 15.mat;t=celiang(:,1);u1=celiang(:,2);syms u2 i5 ii x u5%定义u1,L1,L2,R1,C1,C3,ii 等vU1=celiang(:,2);k=length(u1);C1=1;h=5e-7;C3=1;R1=1;L1=1;R2=1;L2=1;C2=1;i1(1)=0;i1(2)=0;%在开始的两个点用的 怎么 用matlab解符号偏微分方程 matlab 用eig和diag解微分方程 matlab解偏微分方程请把matlab语句给出来～ matlab 求和的出错 matlab求微分方程的问题?clear all;clc;syms k t;result1=int(sym(sin(t)^2/(1+k*sin(t))),t,0,pi)这个没有解析解啊,怎么办兄弟们? 为什么用matlab解微分方程的时候总是出错我输入的命令式dslove('Du=1+u^2','x')结果总是出现Undefined command/function 'dslove'.请问错误出在哪里? matlab 会解偏微分方程吗? matlab求微分方程 matlab解决微分方程问题 Matlab无法求出解析解的微分方程,怎样作出函数图像?比如这个微分方程：dequ='D2x+Dx-x+x^3=cos(t)'求不出解析解,报告出错.但我仍然想作出x关于t的函数图像,我想这时可以通过数值计算,而不是求出 请用Matlab求解下面二阶常微分方程的数值解 求大神帮忙用matlab求解微分方程的数值解