03 2022 档案
摘要：x=np.arange(1,100) y = np.sin(0.1*x) def D_1(x,y): n = len(x) d = np.zeros(n) d[0] = (y[1]-y[0])/(x[1]-x[0]) d[n-1] = (y[n-1]-y[n-2])/(x[n-1]-x[n-2])
阅读全文
摘要：x = [1,2,3] y1 = np.array([4,5,6]) y2 = np.array([5,6,7]) y3 = np.array([7,6,5]) width = 0.2 Sum = y1+y2+y3 per = 100*y2/Sum #添加在中间，旋转90度 def add_mid(
阅读全文
摘要：```python def firstorder(T,K,t,In): # T 时间常数 # K 增益 Out = np.zeros(len(In)) for i in range(1,len(t)): dt = t[i]-t[i-1] s = T/dt Out[i] = K*In[i]/(s+1)
阅读全文
摘要：x = np.arange(0,101) wn = np.random.rand(101) y1 = -1+0.02*x+0.1*wn y2 = np.sin(0.5*x) y3 = 0.0001*x**2-0.03*x+1 y4 = 0.8*np.cos(x) fig,ax1 = plt.subp
阅读全文
摘要：import numpy as np import statsmodels.api as sm from gekko import GEKKO # Data x0 = np.array([4,5,2,3,-1,1,6,7]) x1 = np.array([3,2,3,4, 3,5,2,6]) y =
阅读全文
摘要：import numpy as np from scipy.stats import linregress import statsmodels.api as sm import matplotlib.pyplot as plt from gekko import GEKKO # Data x =
阅读全文
