Python学习手册(第四部分函数)第20章 迭代和解析 第二部分

本文详细介绍了Python中的生成器概念,包括生成器函数和生成器表达式的使用方式,以及它们如何节省内存并允许计算时间分散。文章还展示了如何使用生成器重写内置函数如zip和map,并对比了生成器与列表在迭代行为上的区别。

重访迭代器:生成器

生成器提供了工具在需要的时候才产生结果,而不是立即产生结果。

  1. 生成器函数:编写为常规的def语句,但是使用yield语句一次返回一个结果,而不是立即产生结果。
  2. 生成器表达式:类似于上一小节的列表解析,但是,它们返回按需产生结果的一个对象,而不是构建一个结果列表。

由于二者都不会一次性构建一个列表,它们节省了内存空间,并且允许计算时间分散到各个结果请求。

生成器函数应用 

>>> def gensquares(N):
...     for i in range(N):
...             yield i**2
...
>>> for i in gensquares(5):
...     print(i,end=":")
...
0:1:4:9:16:>>> x = gensquares(2)
>>> x
<generator object gensquares at 0x00000156F3A47D00>
>>> next(x)
0
>>> next(x)
1
>>> next(x)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
StopIteration

生成器表达式:迭代器遇上列表解析

>>> [x**2 for x in range(4)]
[0, 1, 4, 9]
>>> (x**2 for x in range(4))
<generator object <genexpr> at 0x00000156F3A47D00>
>>> G = (x**2 for x in range(4))
>>> next(G)
0
>>> next(G)
1
>>> next(G)
4
>>> next(G)
9
>>> next(G)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
StopIteration
>>> for num in (x**2 for x in range(4)):
...     print("%s,%s"%(num,num/2.0))
...
0,0.0
1,0.5
4,2.0
9,4.5
>>> sum(x**2 for x in range(4))
14
>>> sorted(x**2 for x in range(4))
[0, 1, 4, 9]
>>> sorted((x**2 for x in range(4)),reverse=True)
[9, 4, 1, 0]
>>> import math
>>> list(map(math.sqrt,(x**2 for x in range(4))))
[0.0, 1.0, 2.0, 3.0]

生成器函数vs生成器表达式

生成器是单迭代器

>>> G = (c*2 for c in "spam")
>>> I = iter(G)
>>> next(I)
'ss'
>>> next(I)
'pp'
>>> next(I)
'aa'
>>> I1= iter(G)
>>> next(I1)
'mm'
>>> next(I1)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
StopIteration


# 列表等内置类型支持多个迭代器并且在一个活动迭代器中传递并反映他们的原处修改
>>> L = [1,2,3,4]
>>> I1 = iter(L)
>>> next(I1)
1
>>>
>>> next(I1)
2
>>> I2 = iter(L)
>>> next(I2)
1
>>> next(I2)
2
>>> del L[2:]
>>> next(I1)
Traceback (most recent call last):
  File "<stdin>", line 1, in <module>
StopIteration

编写自己的zip和map

def myzip(*seqs):
    seqs = [list(S) for S in seqs]
    while all(seqs):
        yield tuple(S.pop(0) for S in seqs)

def mymapPad(*seqs,pad=None):
    seqs = [list(S) for S in seqs]
    while any(seqs):
        yield tuple((S.pop(0) if S else pad) for S in seqs)

S1,S2 = 'abc','xyz123'
print(list(myzip(S1,S2)))
print(list(mymapPad(S1,S2)))
print(list(mymapPad(S1,S2,pad=99)))

#结果如下
[('a', 'x'), ('b', 'y'), ('c', 'z')]
[('a', 'x'), ('b', 'y'), ('c', 'z'), (None, '1'), (None, '2'), (None, '3')]
[('a', 'x'), ('b', 'y'), ('c', 'z'), (99, '1'), (99, '2'), (99, '3')]

#或者函数也可以这样写,结果一样
def myzip(*seqs):
    minlen = min(len(S) for S in seqs)
    return [tuple(S[i] for S in seqs) for i in range(minlen)]
def mymapPad(*seqs,pad=None):
    maxlen  = max(len(S) for S in seqs)
    index = range(maxlen)
    return [tuple((S[i] if len(S) > i else pad) for S in seqs) for i in index]

对迭代的各种方式进行计时

# File mytimer.py

import time,sys
if sys.platform[:3] == 'win':
    timefunc = time.clock
else:
    timefunc = time.time

def trace(*args):pass

def timer(func,*pargs,**kargs):
    _reps = kargs.pop('_reps',1000)
    trace(func,pargs,kargs,_reps)
    repslist = range(_reps)
    start = timefunc()
    for i in repslist:
        ret = func(*pargs,**kargs)
    elapsed = timefunc() - start
    return (elapsed,ret)

def best(func,*pargs,**kargs):
    _reps = kargs.pop('_reps',50)
    best = 2**32
    for i in range(_reps):
        (time,ret)=timer(func,*pargs,_reps=1,**kargs)
        if time < best:best = time
    return (best,ret)
# File timeseqs.py

import sys,mytimer
reps = 10000
repslist = range(reps)

def forLoop():
    res = []
    for x in repslist:
        res.append(x+10)
    return res

def listComp():
    return [x+10 for x in repslist]

def mapCall():
    return  list(map((lambda x:x+10),repslist))

def genExpr():
    return list(x+10 for x in repslist)

def genFunc():
    def gen():
        for x in repslist:
            yield x+10
    return list(gen())

print(sys.version)
for tester in (mytimer.timer,mytimer.best):
    print('<%s>'%tester.__name__)
    for test in (forLoop,listComp,mapCall,genExpr,genFunc):
        elapsed,result = mytimer.timer(test)
        print('-'*33)
        print('%-9s:%.5f => [%s...%s]'%(test.__name__,elapsed,result[0],result[-1]))
# file timetests.py

from mytimer import timer,best
def power(X,Y):return X**Y

print(timer(power,2,32))
print(timer(power,2,32,_reps=1000000))
print(timer(power,2,100000)[0])

print(best(power,2,32))
print(best(power,2,100000)[0])
print(best(power,2,100000,_reps=500)[0])

函数陷阱

本地变量是静态检测的

 

 

内容概要:本文围绕“基于分布式模型预测控制的多个固定翼无人机一致性控制”展开,利用Matlab代码实现相关算法的仿真,旨在通过分布式控制策略实现多架固定翼无人机在复杂动态环境中的协同飞行与一致性控制。研究结合模型预测控制(MPC)方法,构建适用于多无人机系统的分布式优化框架,重点解决了通信受限、信息延迟及无中心化指挥条件下的协同稳定性问题。内容涵盖固定翼无人机的动力学建模、分布式MPC优化求解机制、一致性协议设计、通信拓扑结构分析以及仿真验证全过程,确保多机系统在保持队形一致的同时完成协同任务。; 适合人群:具备自动控制理论、无人机系统建模或多智能体协同控制基础,从事智能无人系统、集群控制、自动化与机器人等领域研究的研究生、科研人员及工程技术人员。; 使用场景及目标:①应用于多无人机协同编队飞行、集群侦察、分布式任务执行等实际工程场景;②为分布式MPC算法在多智能体系统中的一致性控制提供可复现的Matlab仿真案例,推动先进控制理论向工程实践转化;③服务于科研论文复现、算法验证、控制系统课程设计与毕业课题参考。; 阅读建议:建议读者结合文中提供的Matlab代码逐模块运行与调试,重点关注分布式MPC在不同通信拓扑下对一致性收敛性能的影响,并可通过调整预测时域、权重矩阵与噪声参数等方式深化对算法鲁棒性与适应性的理解。
内容概要:本文提出了一种基于变分模态分解(VMD)、麻雀搜索算法(SSA)优化与长短期记忆网络(LSTM)相结合的光伏功率预测模型(VMD-SSA-LSTM),旨在提升光伏发电预测的精度与鲁棒性。该方法首先利用VMD对原始非平稳光伏功率序列进行自适应分解,获得一系列具有更稳定特征的本征模态分量(IMFs),有效降低数据复杂性与噪声干扰;随后引入麻雀搜索算法(SSA)对LSTM网络的关键超参数(如学习率、隐层节点数等)进行全局寻优,克服传统试凑法效率低、易陷入局部最优的问题,显著提升模型收敛速度与泛化能力;最后,构建多个LSTM子模型分别预测各模态分量,并将结果重构得到最终的光伏功率预测值。该混合模型充分融合了VMD在信号预处理中的优异分解性能、SSA在参数优化中的高效搜索能力以及LSTM在捕捉时间序列长期依赖关系上的强大建模优势,实现了对复杂气象因素影响下光伏出力波动的高精度拟合与预测。; 适合人群:具备一定电力系统、新能源发电或时间序列预测基础知识,熟悉MATLAB编程环境,从事光伏功率预测、智能电网调度、可再生能源集成、负荷预测等领域研究的科研人员、工程技术人员及高校研究生。; 使用场景及目标:①应用于光伏电站的短期与超短期功率预测,为电网安全调度、电力市场交易、储能系统配置及需求侧响应提供精准数据支撑;②解决传统单一预测模型(如ARIMA、BPNN、单一LSTM)在处理非平稳、强波动性光伏数据时存在的精度不足、稳定性差等问题;③为风电、负荷等其他非平稳时序预测问题提供一种有效的“分解-优化-预测”混合建模范式与技术实现路径。; 阅读建议:建议读者结合文中提供的完整MATLAB代码,深入理解VMD信号分解、SSA优化算法流程及LSTM网络构建的每一个技术环节,通过实际历史数据进行模型复现与对比实验(如与VMD-LSTM、SSA-LSTM等模型比较),掌握参数调优技巧与模型性能评估方法,从而真正掌握该先进混合预测模型的核心思想与应用精髓。
评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值