2021-10-01から1ヶ月間の記事一覧

汎用旋盤の解説

www.youtube.com

レーベンバーグ=マーカート法

www.youtube.com

Matlabでベクトル場Aの全微分dAに関する公式の計算

以下を計算で確かめる。 ただし、Sは以下で定義する。 等式の右辺、1/2(rot A) x dx + S • dx を計算し、dAに一致することを示す。 clear; close all; clc; syms Ax Ay Az ddx ddy ddz dx dy dz real; A = [Ax, Ay, Az]'; nabla=[ddx, ddy, ddz]'; gradA=nab…

AlfredでチートシートPDFを開く

以下のPDF File Searcherを少し改造した。cheat コマンド(あるいはOption+C)で起動するようにし、Dropbox/Cheatsheet以下を検索するようにした。 www.alfredapp.com

テンソル積とクロネッカー積の違い

テンソル積はa⊗b、クロネッカー積はA⊗bとなっており、クロネッカー積では左に行列が来る。 math.1book.info

ダイクストラ法で最短経路を求める

qiita.com

RASPフォーマット(.engファイル)の解析

showEng.py #!/usr/bin/env python import pandas as pd import matplotlib.pyplot as plt from scipy import integrate import sys import warnings warnings.filterwarnings("ignore") filename = sys.argv[1] metadata = pd.read_csv(filename, delimiter…

モデルロケットのエンジンの型番

https://www.apogeerockets.com/education/downloads/Newsletter486.pdf アルファベットはトータルインパルス(N・s) 数字は平均推力(N)

Pythonで定積分を求める

qiita.com

Q学習のアルゴリズム

https://www.tcom242242.net/entry/ai-2/%E5%BC%B7%E5%8C%96%E5%AD%A6%E7%BF%92/%E3%80%90%E5%BC%B7%E5%8C%96%E5%AD%A6%E7%BF%92%E3%80%81%E5%85%A5%E9%96%80%E3%80%91q%E5%AD%A6%E7%BF%92_%E8%BF%B7%E8%B7%AF%E3%82%92%E4%BE%8B%E3%81%AB/

SARSAとQ学習の違い

今さら聞けない強化学習(10): SarsaとQ学習の違い - Qiita

重心法によるクリスプ値の計算

import numpy as np import skfuzzy as fuzz from skfuzzy import control as ctrl import matplotlib.pyplot as plt import pdb credit = ctrl.Antecedent(np.arange(0, 11, 0.1), 'credit') credit['low'] = fuzz.trimf(credit.universe, [0,2,4]) credit[…

トラスの最適化 その2

github.com

トラスの最適化

He, Linwei, Matthew Gilbert, and Xingyi Song. "A Python script for adaptive layout optimization of trusses." Structural and Multidisciplinary Optimization 60.2 (2019): 835-847. Pythonスクリプト figshare.com from math import gcd, ceil impor…

焼き鈍し法(Simulated Annealing:SA)の実装

OneMax問題(10ビット)を解いてみる。 温度を線形に減少する。 import numpy as np import matplotlib.pyplot as plt R = 10 T = 10 N_bit = 10 x = [np.random.randint(100)%2 for _ in range(N_bit)] y = np.zeros_like(x) def evalOneMax(x): return sum…

pipを使う際のtips

結論:python -m pip install を使う 注意点: 1. sudoは使わない(権限の関連で問題が生じる) 2. pip installは使わない(PATHに読み込んでいるはずのsite-packages以下にインストールされたモジュールが読めないというエラーが発生する恐れがある)

山登り法の実装

OneMax問題(10ビット)を解く。 import numpy as np import matplotlib.pyplot as plt R = 500 N_bit = 10 x = [np.random.randint(100)%2 for _ in range(N_bit)] y = np.zeros_like(x) def evalOneMax(x): return sum(x) def mutate(x): idx = np.rando…

OpenRocketをPythonから呼び出す

ここを参考に。orhelperを使う。 wiki.openrocket.info 自分の環境では以下が必要だった。 import matplotlib matplotlib.use('Qt5Agg') 以下のサンプルコードsimple_plot.pyでは、OpenRocket-15.03.jarとsimple.orkをカレントディレクトリにおくことを想定…

粒子群最適化(PSO)の実装

OneMax問題(10ビット)を解いてみる。 import numpy as np import random import pdb import matplotlib.pyplot as plt T = 100 N_bit = 10 N = 20 x = np.random.rand(N,N_bit) v = np.random.rand(N,N_bit) g = np.zeros(N_bit) p = np.zeros((N,N_bit)) …

wxPythonを使う(PyenvでFramework Pythonのインストール)

qiita.com以下のコマンドでPyenvを使ってPythonをインストールする。 $ env PYTHON_CONFIGURE_OPTS="--enable-framework" pyenv install 3.9.2

Sfepyのサンプルコード

https://sfepy.org/doc-devel/tutorial.html #!/usr/bin/env python from __future__ import print_function from __future__ import absolute_import from argparse import ArgumentParser import numpy as nm import sys sys.path.append('.') from sfepy.…

Sfepyのインストール

追記:よくわからんけど色々やっていたら動いた。 最終的には export ETS_TOOLKIT="qt" が鍵だった。と思ったら今度はこの環境変数があると動かなくなった(unsetして消した)。謎が多い。 とりあえず色々エラーの原因になるようなのでAnaconda環境をすべて削…

Conda-forgeからのインストール

waregawa-log.hatenablog.com $ conda config --append channels conda-forge $ conda install hogehoge

DEAPで遺伝的プログラミング

import numpy from deap import algorithms, base, creator, tools, gp creator.create("FitnessMin", base.Fitness, weights=(-1.0,)) creator.create("Tree", gp.PrimitiveTree, fitness=creator.FitnessMin) pset = gp.PrimitiveSet(name="MAIN", arity=1…