早教吧作业答案频道 -->数学-->
matlab 求解微分方程组已知u1如何求u2i1=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;
题目详情
matlab 求解微分方程组
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;
已知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;
▼优质解答
答案和解析
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
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 求解微分方程组...的网友还看了以下:
随机过程题目:设X是一连续随机变量,具有分布F,证明:(a)F(x)服从(0,1)上的均匀分布随机 2020-04-13 …
定义在R上的函数F(X)满足关系式F(二分之一+X)+F(二分之一减区X)=2,则F(八分之一)+ 2020-04-26 …
已知函数f(x)满足对任意的x属于R都有f(二分之一+x)+f(二分之一减x)=2成立,则f(八分 2020-04-26 …
一道简单的大一高数问题已知f(π)=1,f(x)二阶导数连续,且∫上x下0[f(x)+f""(x) 2020-05-13 …
把线段P1P2在P3点平分,线段P2P3在P4点平分,P3P4在P5点平分,这样不断地平分法,当n 2020-05-14 …
f(x)=Asin(ωx+φ)问题.简单,有分!其中A,ω,φ在图形分别代表什么.f(x)=sin 2020-05-14 …
(高一函数) f(x)-f(-x) f(-x)-f(x) f(x)+f(-x) f(x)f(-x) 2020-05-16 …
请问氢键是怎么判断的?是不是一定要H-F,H-N和H-O才算是氢键.如果C-F这个不是真的键,是2 2020-05-21 …
f(-|x|)与f(|x|)分别由f(x)怎样得来?这是我的理解,你们仔细看看对不对.函数y=f( 2020-06-03 …
22、已知f{x}是定义在上的奇函数,且f(1)=1,若a,b,a+b,有.⑴证明:函数f(x)在 2020-06-08 …