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

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

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; …