Navigation Menu

Skip to content

Instantly share code, notes, and snippets.

View ina111's full-sized avatar

ina111 ina111

View GitHub Profile
# -*- coding: utf-8 -*-
import sys
reload(sys)
import platform
# デフォルトの文字コードを変更する.
sys.setdefaultencoding('utf-8')
import numpy as np
import matplotlib.pyplot as plt
import matplotlib.font_manager
Sorry, something went wrong. Reload?
Sorry, we cannot display this file.
Sorry, this file is invalid so it cannot be displayed.
@ina111
ina111 / matplotlib_PDF_ja.py
Last active March 21, 2016 12:04
MacでMatplotlibで日本語の画像とPDFを作る
# -*- coding: utf-8 -*-
"""
Mac OSXのmatplotlibで日本語フォント(Osaka)を使うサンプル
python2.7(anaconda-2.1.0)
cf. http://matplotlib.org/users/customizing.html
"""
# デフォルトの文字コードを変更する.
import sys
reload(sys)
sys.setdefaultencoding('utf-8')
@ina111
ina111 / visible_range.py
Last active December 3, 2016 14:35
衛星などの上空にある物体の地上からの可視範囲のグラフ化
# -*- coding: utf-8 -*-
# 衛星高度と可視範囲のグラフ化
#
# cf. 半揚 稔雄(2014) 「ミッション解析と軌道設計の基礎」
# Copyright (c) 2016 Takahiro Inagawa
# Released under the MIT license
import sys
reload(sys)
import platform
sys.setdefaultencoding('utf-8')
@ina111
ina111 / Bokeh_sample.py
Last active September 11, 2017 06:43
pythonのBokehというインタラクティブな可視化ライブラリの基本的な機能のサンプルです。ここではbokeh.plottingインターフェイスを使用しています。
# -*- coding: utf-8 -*-
# Author : Takahiro Inagawa(@ina111)
import io
import numpy as np
from jinja2 import Template
from bokeh.embed import components
from bokeh.models import Range1d
% -----------
% Optimum Design of Nonplanar wings-Minimum Induced Drag
% for A Given Lift and Wing Root Bending Moment (NAL TR-797)
%
% Created by Takahiro Inagawa on 2013-2-28.
% Copyright (c) 2013 Takahiro Inagawa. All rights reserved.
% -----------
clear all
tic; %Computation time measurement
# -*- coding: utf-8 -*-
# -----------
# Optimum Design of Nonplanar wings-Minimum Induced Drag
# for A Given Lift and Wing Root Bending Moment (NAL TR-797)
#
# Created by Takahiro Inagawa on 2018-03-24.
# Copyright (c) 2018 Takahiro Inagawa. All rights reserved.
# -----------
import numpy as np

##概要 JAXA(旧NAL)の1984年の航空宇宙技術研究所報告(NAL TR-797)

「非平面翼の最適設計-揚力と翼根曲げモーメントを与えた時の最小誘導抵抗-」

の数値計算部分をMatlab(Octaveで動作確認)で実装しました. 人力飛行機を念頭にした数値を入れてあります. 論文は汎用的なものなのでConstantの部分の数字を変えれば汎用的に使えます.

論文中では非平面翼の計算を行なっていますが,ここに置いているスクリプトは平面翼です.

This file has been truncated, but you can view the full file.
{
"cells": [
{
"cell_type": "markdown",
"metadata": {},
"source": [
"## REFPROPのPythonラッパーの使い方\n",
"\n",
"本家NISTから公開されているラッパーの使い方のjupyter notebookはこちら\n",
"* https://nbviewer.jupyter.org/github/usnistgov/REFPROP-wrappers/blob/master/wrappers/python/notebooks/Tutorial.ipynb\n",
@ina111
ina111 / high_altitude_wind_from_service.py
Last active September 8, 2021 07:09
気象庁の高層風のページからデータをスクレイピングして、プロット
# coding: UTF-8
# 気象庁の高層風のページからデータ取得してプロット
# cf. http://www.data.jma.go.jp/obd/stats/etrn/upper/index.php
import os
import urllib2
from BeautifulSoup import BeautifulSoup
import datetime
import pandas as pd
import numpy as np