2019-01-01から1年間の記事一覧

ANSYSで対象性を用いた計算の全体結果を表示する(未完了)

studentcommunity.ansys.comANSYS WorkbenchのOptionでBeta optionsを有効にする SymmetryからGraphical expansion(Beta)のNum Repeatを2に、MethodをHalfに、⊿Zを1e-3に設定してみたが、表示結果は変わらなかった。

Bézier曲線とB-spline曲線

clear; close all; figure();grid(); xlim([-10 10]);ylim([-10 10]); [x,y]=ginput(4); x0=x(1);x1=x(2);x2=x(3);x3=x(4); y0=y(1);y1=y(2);y2=y(3);y3=y(4); Px1=[];Py1=[]; Px2=[];Py2=[]; for t=0:0.01:1 %Bézier B0=(1-t)^3; B1=3*(1-t)^2*t; B2=3*(1-t…

NURBSのシミュレータ

NURBS demo - WebGL based online evaluator for NURBS Curveshttp://geometrie.foretnik.net/files/NURBS-en.swf

ROS2のインストール

index.ros.org $ brew install python3 tinyxml2 tinyxml eigen pcre opencv openssl qt freetype assimp log4cxx

デジタルフィルタによる補間後加減速

Simulinkブロック Linear acc/dec 線形加減速(台形型加減速) Exponential acc/dec 指数型加減速 S-curve type acc/dec ベル型加減速

たたみこみによるデジタルフィルタの計算

clear; close all; h=[1 4 -2]; x=[2 -1 3 1]; hlength=length(h); xlength=length(x); ylength=hlength+xlength-1; hzero=[h zeros(1,xlength-1)]; xzero=[x zeros(1,hlength-1)]; y=zeros(1,ylength); for n=1:ylength-1 for k=0:n y(n+1)=y(n+1)+hzero(k+…

macOS CatalinaでPapers3がクラッシュして起動しない

以下から修正版をダウンロードできるreadcubesupport.freshdesk.com

MOSFETのゲート抵抗について

blog.handen.net

MOSFETの寄生容量について

電気双極子の作る電界 その2

clear all; close all; q01=2; q02=-2; e0=8.85*1e-12; k=1/(4*pi*e0); omega=1; r1=[-2,0]'; r2=[2,0]'; x=-5:5; y=-5:5; [X,Y]=meshgrid(x,y); a1=sqrt((X-r1(1)).^2+(Y-r1(2)).^2); a2=sqrt((X-r2(1)).^2+(Y-r2(2)).^2); myVideo = VideoWriter('myVideoF…

電気双極子モーメント

clear all; close all; q1=2; q2=-2; e0=8.85*1e-12; k=1/(4*pi*e0); x=-5:5; y=-5:5; r1=[-2,0]'; r2=[2,0]'; [X,Y]=meshgrid(x,y); a1=sqrt((X-r1(1)).^2+(Y-r1(2)).^2); a2=sqrt((X-r2(1)).^2+(Y-r2(2)).^2); E1x=q1*k./a1^3.*(X-r1(1)); E1y=q1*k./a1^3.…

無限直線電流による磁界の発生(アンペール/ビオ=サバールの法則)

clear all; close all; I=1; x=-5:5; y=-5:5; z=-5:5; [X,Y,Z]=meshgrid(x,y,z); Hx=-I/(2*pi)*Y./(X.^2+Y.^2); Hy=I/(2*pi)*X./(X.^2+Y.^2); Hz=0*X; plot3([0,0],[0,0],[z(1),z(end)],'r-','linewidth',5); hold on; quiver3(X,Y,Z,Hx,Hy,Hz,'b');

電気双極子の作る電界

clear all; close all; % q1=2; % q2=-2; e0=8.85*1e-12; k=1/(4*pi*e0); q1s=[1,2,3]; q2s=[-3,-2,-1]; x=-5:5; y=-5:5; r1=[-2,0]'; r2=[2,0]' [X,Y]=meshgrid(x,y); a1=sqrt((X-r1(1)).^2+(Y-r1(2)).^2); a2=sqrt((X-r2(1)).^2+(Y-r2(2)).^2); for i1=1:l…

点電荷の作る電界

clear all; close all; q=2; x=-10:10; y=-10:10; [X,Y]=meshgrid(x,y); Z=sqrt(100-X.^2-Y.^2); Ex=100*X./(X.^2+Y.^2+Z.^2).^(3/2); Ey=100*Y./(X.^2+Y.^2+Z.^2).^(3/2); Ez=100*Z./(X.^2+Y.^2+Z.^2).^(3/2); Ex(isnan(Ex))=0; Ey(isnan(Ey))=0; Ez(isnan(…

電気力線

clear all; close all; q=2; r=[[1,0]',[0,1]',[-1,0]',[0,-1]',[1.5,1.5]',[-1.5,1.5]',[-1.5,-1.5]',[1.5,-1.5]']; r2=arrayfun(@(x,y)x.^2+y.^2,r(1,:),r(2,:)) r1=sqrt(r2); e0=8.85*1e-12; E=1/(4*pi*e0)*q./r2./r1.*r; plot(0,0,'bo'); hold on; array…

静電気力

clear all; close all; qA=-2; qB=1; rA=[2,1,0]'; rB=[0,1,0]'; r=[rA rB]; e0=8.85*1e-12; k = 1/4/pi/e0; F = k*qA*qB/sqrt(sum((rA-rB).^2)); Fdummy=1;%correction for view F=Fdummy;%The real F is too huge to be shown, so just leave it v=rA-rB; …

2つのmp4動画をコマンドラインで結合

stackoverflow.com cat mylist.txt file '/path/to/file1' file '/path/to/file2' file '/path/to/file3' $ ffmpeg -f concat -safe 0 -i mylist.txt -c copy output.mp4

Youtubeで字幕付き動画をダウンロードする

自動生成された英語のキャプションを動画にくっつけてダウンロードしたい 1. youtub-dlで字幕と動画をダウンロード(字幕はsrt形式: --convert-subtitlesオプション) youtube-dl -f mp4 --write-auto-sub --convert-subtitles srt --sub-lang en https://ww…

MojaveでC関連のヘッダが見つからないエラー

github.com 以下を実行してヘッダをインストール open /Library/Developer/CommandLineTools/Packages/macOS_SDK_headers_for_macOS_10.14.pkg

ANSYSでFEMのテスト

バスの車内で緊急脱出時には窓ガラスの端を割れという表示を見かけた 中央と端でどれだけ応力集中が異なるかのシミュレーション端にforceをかけたとき von Mises応力 ひずみ 中央にforceをかけたとき von Mises応力 ひずみ 静荷重だと端の方が割れやすいとい…

Pintle injectorの設計

ここを参考に。 pintleinjector.blogspot.com pintleinjector.blogspot.com 1. 推力と推進剤の選定 例題では5000N, RP1-LOX 。仮にこれでやってみる。 2. 燃焼室圧力の選定 例題では1.5Mpa。 3. 最適混合比の同定 RP-1/LOXの凍結流のデータ 燃焼室圧力が異な…

ノズル開口比の選定

なぜ海面上剥離限界で設計するのか? 1)ノズル背圧Pa(大気圧)が上がっていくと、不足膨張状態となり、ノズル後方で衝撃波が発生することで膨張が完了する。どんどんノズル背圧をあげていくと、ノズル出口でちょうど衝撃波が立つ状態となる。これを剥離限…

Mojave+VM FusionでUSBが認識されない

Full disk accessをONにしたら認識されるようになった。 https://communities.vmware.com/thread/598488

Apacheが再起動でコケる問題

Apacheの再起動で以下のエラーが出る $ sudo /etc/init.d/httpd restart(98)Address already in use: make_sock: could not bind to address 0.0.0.0:8080番PortをListenしているプロセスがいるので、そのプロセスを停止してから再起動する https://ubuntufo…

MacのPreview.app, Keynote.appでファイルがリネームできない問題

長年の問題が解決した forums.macrumors.com以下をターミナルで実行する defaults write NSGlobalDomain NSAutomaticWindowAnimationsEnabled -bool YES

deviseの使い方

いつも忘れるのでメモ qiita.com

CourseraのDeep Learning Specialization

途中やめになってしまったCNNとSequence modelを完了 取り掛かりから1年かかってしまったが終えることができた https://www.coursera.org/account/accomplishments/specialization/4HXUWQWGAB3P

Super resolutionを試す

github.com入力 出力

ターミナルでGoogle driveからファイルを直接ダウンロードする

wgetを使うが、Cookie許可とConfirm codeの取得が必要 unix.stackexchange.com $ wget --save-cookies cookies.txt --keep-session-cookies --no-check-certificate 'https://docs.google.com/uc?export=download&id=FILEID' -O- | sed -rn 's/.*confirm=([0…

darknetでYOLO v3を試す(GPU, OpenCVつき)

darknetはこっちのレポジトリをつかう github.comGPUとOpenCVを1にしてmake(cudnnを1にするとビルドできるが動かなかった) makefile GPU=1 CUDNN=0 CUDNN_HALF=0 OPENCV=1 AVX=0 OPENMP=0 LIBSO=0 ZED_CAMERA=0コマンド dont_showをつけないといちいち表…