m0_51078708 2020-09-24 15:33 采纳率: 0%
浏览 205

pyhton 用欧拉方法解状态方程x’=Ax+B,A和B都是矩阵

x’(t)=[0,-1;1,-2]x(t)+[0;2],x(0)=[1;2],t属于[0,1),在python中用欧拉法求解,T取0.1。

我写的代码
import numpy as np
from matplotlib.pyplot import *

Euler's method with integration step size T=0.1

T = 0.1 # step size
t = np.arange(0.,1,0.1) # time instants that the solution will be evaluated
u = np.array([[0];[1]])
x0[0]=np.array([[1];[2]])
def f(x0,u):
x=np.zeros(shape=(2,1))
x[0]=-1*x0[1]
x[1]=1*x0[0]-2*x0[1]+u
return np.squeeze(x)
If k in range(1,9):
x0[k]=x0[k-1]+f(x0[k-1])*T

我无论怎么改都会说维度有问题或者未命名函数,我实在是不知道怎么改了,求解答

  • 写回答

1条回答 默认 最新

  • threenewbee 2020-09-24 18:03
    关注
    评论

报告相同问题?

悬赏问题

  • ¥65 永磁型步进电机PID算法
  • ¥15 sqlite 附加(attach database)加密数据库时,返回26是什么原因呢?
  • ¥88 找成都本地经验丰富懂小程序开发的技术大咖
  • ¥15 如何处理复杂数据表格的除法运算
  • ¥15 如何用stc8h1k08的片子做485数据透传的功能?(关键词-串口)
  • ¥15 有兄弟姐妹会用word插图功能制作类似citespace的图片吗?
  • ¥15 latex怎么处理论文引理引用参考文献
  • ¥15 请教:如何用postman调用本地虚拟机区块链接上的合约?
  • ¥15 为什么使用javacv转封装rtsp为rtmp时出现如下问题:[h264 @ 000000004faf7500]no frame?
  • ¥15 乘性高斯噪声在深度学习网络中的应用