Python和MATLAB谐波生成导图

🎯要点

  1. 绘制三次谐波生成透射功率谱、对数对数图表示
  2. 半导体曲面二次谐波生成,分析判断材料特性
  3. 谐波均值估计计算边际似然(贝叶斯统计)
  4. 二次谐波散射分析胶体染料分子结构
  5. 交流电谐波波形傅立叶分析
  6. 分析旋转各向异性谐波
  7. 高次谐波非线性相互作用过程分析

谐波生成导图

在这里插入图片描述

通用谐波方程

Python 脚本建模通用谐波方程。请注意,方波仅具有奇次谐波。
y ( t ) = ∑ k = 1 N + 1 k − 1 sin ⁡ ( 2 ω k f t ) y ( t )=\sum_{k=1}^{ N +1} k^{-1} \sin (2 \omega kft ) y(t)=k=1N+1k1sin(2ωkft)

import matplotlib.pylab as plt
import numpy as np
import argparseVERSION = '0.2'def factor(shape, k, i):if shape == "triangle":return (1/(k*k) * (-1)**i)else:return (1/k)if __name__ == '__main__':parser = argparse.ArgumentParser(description='plot harmonics',epilog='Version: ' + VERSION)parser.add_argument('-f','--frequency',default=4,type=int, help='specify the frequency',action='store')parser.add_argument('-n','--harmonics',default=0,type=int, help='specify the number of harmonics',action='store')parser.add_argument('-t','--type',default='odd',type=str, help='specify the type of harmonics (even, odd, or all)',action='store')parser.add_argument('-s','--shape',default='square',type=str, help='specify the shaping factor (square, triangle)',action='store')args = parser.parse_args()f = args.frequencyif args.type == 'even':odd = 0mult = 2elif args.type == 'odd':odd = 1mult = 2else:   odd = 0mult = 1t = np.linspace(0, 1, num=8000)y = np.zeros(8000)for i in range(int(args.harmonics)+1):k = i * mult + oddyh = factor(args.shape,k,i) * np.sin(2 * np.pi * k *  f * t)y = y + yhplt.plot(t, y)plt.xlabel('time')plt.ylabel('harmonics {0}'.format(args.harmonics))plt.axis('tight')plt.show()print('finished')

脚本命令行参数:

$ python3 harmonic.py -h
usage: harmonics.py [-h] [-f FREQUENCY] [-n HARMONICS] [-t TYPE] [-s SHAPE]plot harmonicsoptional arguments:-h, --help            show this help message and exit-f FREQUENCY, --frequency FREQUENCYspecify the frequency-n HARMONICS, --harmonics HARMONICSspecify the number of harmonics-t TYPE, --type TYPE  specify the type of harmonics (even, odd, or all)-s SHAPE, --shape SHAPEspecify the shaping factor (square, triangle)

生成基波(正弦波)

$ python3 harmonic.py

生成基波(正弦波)+ 2 次谐波。我们可以看到方波正在形成,并带有大的波纹。

$ python3 harmonic.py -n 2

让我们调高谐波,看看我们开始近似理想的方波。

$ python3 harmonic.py -n 1000

产生三角波。请注意,由于谐波的缩放比例不同,方程现在有所不同。

$ python3 harmonic.py -n 1000 -s triangle

👉更新:亚图跨际

本文来自互联网用户投稿,该文观点仅代表作者本人,不代表本站立场。本站仅提供信息存储空间服务,不拥有所有权,不承担相关法律责任。如若转载,请注明出处:http://www.mzph.cn/diannao/52167.shtml

如若内容造成侵权/违法违规/事实不符,请联系多彩编程网进行投诉反馈email:809451989@qq.com,一经查实,立即删除!

相关文章

梧桐数据库(WuTongDB):数据库技术中LL算法详解

LL 算法是一种自顶向下的语法分析算法,广泛用于构建解析器。LL 分析器逐个读取输入符号,从左到右分析(Left-to-Right),并使用最左推导(Leftmost Derivation)来生成语法树。因此,LL 分…

appium下载及安装

下载地址:https://github.com/appium/appium-desktop/releases 双击安装就可以

深入学习SQL优化的第三天

目录 聚合函数 排序和分组 聚合函数 1251. 平均售价 表:Prices------------------------ | Column Name | Type | ------------------------ | product_id | int | | start_date | date | | end_date | date | | price | int …

桌球厅助教陪练系统源码开发和行业市场分析

台球助教陪练系统:引领智能化运动体验 作为一款专为台球爱好者设计的智能陪练系统,我们的目标是通过技术创新,让每位用户都能享受到个性化、高效的学习体验。无论是初学者还是寻求突破的高手,都能在我们的平台上找到适合自己的陪…

MybatisPlus的LambdaQueryWrapper用法

LambdaQueryWrapper<Tag> 是 MyBatis-Plus 框架中的一个功能强大的查询构造器&#xff0c;它用于构建 SQL 查询条件&#xff0c;特别是针对 Lambda 表达式的支持&#xff0c;使得代码更加简洁、类型安全。在这个例子中&#xff0c;LambdaQueryWrapper<Tag> 被用来构…

基于YOLO V8的PCB缺陷检测识别系统(python源码+Pyqt5界面+数据集+训练代码)

数据集准备&#xff1a;收集并标注PCB缺陷的图像。模型训练&#xff1a;使用YOLO v8框架训练一个模型来识别这些缺陷。GUI开发&#xff1a;利用PyQt5创建一个用户友好的图形界面。模型部署&#xff1a;在GUI中集成训练好的模型&#xff0c;使用户能够上传PCB图像并得到缺陷检测…

IOS 11 通用Base控制器封装

整体规划 BaseController&#xff1a;把viewDidLoad逻辑拆分为三个方法&#xff0c;方便管理。 BaseCommonController&#xff1a;不同项目可以复用的逻辑&#xff0c;例如&#xff1a;设置背景颜色方法等 BaseLogicController&#xff1a;本项目的通用逻辑&#xff0c;主要…

实现 FastCGI

CGI的由来&#xff1a; 最早的 Web 服务器只能简单地响应浏览器发来的 HTTP 请求&#xff0c;并将存储在服务器上的 HTML 文件返回给浏 览器&#xff0c;也就是静态 html 文件&#xff0c;但是后期随着网站功能增多网站开发也越来越复杂&#xff0c;以至于出现动态技 术&…

【Pyhthon读取 PDF文件表格 ,转为 CSV/TSV/JSON文件】

tabula-py tabula-py 是一个将 PDF 表格转换为 pandas DataFrame 的工具。 tabula-py 是 tabula-java 的包装器&#xff0c;需要您的机器上有 java。 tabula-py 还允许您将 PDF 中的表格转换为 CSV/TSV 文件。 tabula-py 的 PDF 提取准确度与 tabula-java 或 tabula app 相…

JavaScript基础知识(六)

相关api介绍 数字类型 - parseInt - parseFloat 字符串类型 .length 返回字符串长度 字符串可以通过下标的方式来获取值 字符串的相关api都不会改变原来字符串,并且可以返回一个新的字符串,若要获取其值则需要声明新的变量获取 str.indexOf(char[,startIndex]) 查询st…

8/21作业

一、 非阻塞型IO 让我们的read函数不再阻塞&#xff0c;无论是否读取到消息&#xff0c;立刻返回 1.1 fcntl函数 原型&#xff1a;int fcntl(int fd, int cmd, ... /* arg */ ); 调用&#xff1a;int flag fcntl(描述符,F_GETFL) fcntl(描述符&#xff0c;F_SETFL&…

洛谷p10892题解

题目背景 AzureHair 在 NOIP 2022 中被喵了个喵创死&#xff0c;于是患上了不治之症——T2 恐惧症&#xff0c;于是他在 NOIP 2023 中果断跳过了 T2 并杠 T3 两小时无果&#xff0c;遗憾离场&#xff0c;他的同学决定帮他治疗这种不治之症。 在他的同学给他治愈了 T2 恐惧症后…

机器学习-KNN 算法

一.K-近邻(KNN) K-近邻&#xff08;K-Nearest Neighbors, 简称 KNN&#xff09;是一种基于实例的学习算法&#xff0c;主要用于分类和回归问题。KNN 的工作原理直观且简单&#xff0c;它基于相似性进行预测&#xff0c;也就是说给定一个新的数据点&#xff0c;KNN 算法会查找距…

JUC知识点总结

JUC应用场景 1. 网页服务器处理并发请求 当一个网页服务器需要处理大量并发请求时&#xff0c;可以使用多线程来提高处理效率。 import java.io.IOException; import java.net.ServerSocket; import java.net.Socket; ​ public class WebServer {public static void main(S…

8月21日笔记

Frp Frp(Fast e Reverse ) Proxy) 是一款简单&#xff0c;好用&#xff0c;稳定的隧道工具。Frp 使用 Go语言开发&#xff0c;支持跨平台&#xff0c;仅需下载对应平台的二进制文件即可执行&#xff0c;没有额外依赖。它是一款高性能的反向代理应用&#xff0c;可以轻松地进行…

Spring DI 数据类型—— set 方法注入

首先新建项目&#xff0c;可参考 初识IDEA、模拟三层--控制层、业务层和数据访问层 一、spring 环境搭建 &#xff08;一&#xff09;pom.xml 导相关坐标 <?xml version"1.0" encoding"UTF-8"?> <project xmlns"http://maven.apache.or…

http连接未释放导致生产故障

凌晨4点运维老大收到报警&#xff08;公司官网页面超时&#xff0c;上次故障因为运维修改nginx导致官网域名下某些接口不可用后&#xff0c;运维在2台nginx服务器上放了检测程序&#xff0c;检测官网页面&#xff09;&#xff0c;运维自己先看了看服务器相关配置&#xff0c;后…

Java实现STL中的全排列函数next_permutation()

目录 一、引言 二、全排列函数next_permutation() 三、next_permutation()的使用 四、Java实现next_permutation() 五、使用next_permutation()实现全排列 一、引言 相信很多小伙伴们都做过全排列的算法题&#xff0c;输入一个n&#xff0c;输出1~n的全排列。对于这个问题…

k8s相关命令

一、Kubectl是什么 控制K8S通信的命令工具。 格式&#xff1a; kubectl [command] [TYPE] [NAME] [FLAGS] 二、Kubectl命令 1.查看容器命令 kubectl get pod kubectl get pod -n default #查看当前的命名空间的pod kubectl get pod -A #查看所有命名空间的pod kubectl de…

.Net插件开发开源框架

在.NET开发中&#xff0c;有许多开源框架可以用于插件开发&#xff0c;以下是一些最常见的框架&#xff1a; MEF&#xff08;Managed Extensibility Framework&#xff09; MEF是一个用于创建可插拔软件应用程序的库&#xff0c;它可以在不修改原始应用程序的情况下扩展应用程…