monhime

N = 300; %r方向の格子数
M = 300; %z方向の格子数
RMAX = 30; %r座標の最大値
ZMAX = 30; %z座標の最大値
Conv = 1e-6; %収束判定誤差
dr = RMAX / N; %r方向の刻み幅
dz = ZMAX / M; %z方向の刻み幅
phi = zeros(N,M); %0埋め
flag = -10*ones(N,M); %非境界部分の-10埋め

monhime

N = 800; %r方向の格子数
M = 800; %z方向の格子数
RMAX = 30; %r座標の最大値
ZMAX = 30; %z座標の最大値
Conv = 1e-6; %収束判定誤差
dr = RMAX / N; %r方向の刻み幅
dz = ZMAX / M; %z方向の刻み幅
phi = zeros(N,M); %0埋め 
flag = -10*ones(N,M); %非境界部分の-10埋め

monhime

%  tdata: 0から20まで100等分した値
tdata = linspace(0, 20, 100);

% ydata: 計測して得られた物体の位置（を想定）
% ここでは理想のモデル式に乱数のノイズを足しています
ydata = 12.5*exp(-0.2*tdata) .* sin(tdata) + 0.3*randn(size(tdata));


% 係数ベクトルtと、tdataを入力とする関数を作成
% 考えているモデル式の係数の部分を係数で置き換えます

monhime

% 10から10^6までログスケールで等間隔な100個の点ω [rad/s]
omega = logspace(1,6,100)

% Gdata: 計測して得られた利得（を想定）
% ここでは理想のモデル式に乱数のノイズを足しています
Gdata = 20*log10(abs(1000./(i*omega+1000))) + randn(100,1)


G = @(k,omega) 20*log10(abs(k(1)./(i*omega + k(2))))

monhime

% V: 印加した電圧のベクトル
V = [74.121; 88.9452; 93.8866; 96.3573; 97.59265; 98.828; 101.2987; 103.7694; 106.2401; 108.7108; 111.1815; 113.6522; 123.535] 
% P: その電圧での火花放電率  
P = [0; 0; 0.2; 0.2; 0.2; 0.1; 0.4; 0.7; 0.7; 0.8; 1; 1; 1]

% 係数ベクトルkと、入力vの関数を作成
% normcdf(x,m,v)は平均値m, 標準偏差vの正規累積分布のxのときの値を出力します
% k(1)が平均値, k(2)が標準偏差を表します
p = @(k,vdata) normcdf(vdata,k(1),k(2)); 

monhime

\documentclass[10pt, a4paper, titlepage, uplatex]{jsarticle} 
\usepackage{amsmath,ams symb}
\usepackage[dvipdfmx]{graphicx} 
\usepackage{circuitikz}

\begin{document}
	\begin{figure}[htbp]
		\begin{center}
		\begin{circuitikz}[american currents]
			\draw(0,0)

monhime

\documentclass[10pt, a4paper, titlepage, uplatex]{jsarticle} 
\usepackage{amsmath,ams symb}
\usepackage[dvipdfmx]{graphicx} 
\usepackage{circuitikz}

\begin{document}
	\begin{figure}[htbp]
		\begin{center}
		\begin{circuitikz}[american currents]
			\draw(0,0)

monhime

\documentclass[10pt, a4paper, titlepage, uplatex]{jsarticle} 
\usepackage{amsmath,ams symb}
\usepackage[dvipdfmx]{graphicx} 
\usepackage{circuitikz}

\begin{document}
	\begin{figure}[htbp]
%三相交流発電機
		\begin{center}
		\begin{circuitikz}[american currents]

monhime

\documentclass[10pt, a4paper, titlepage, uplatex]{jsarticle} 
\usepackage{amsmath,ams symb}
\usepackage[dvipdfmx]{graphicx} 
\usepackage{circuitikz}

\begin{document}
	\begin{figure}[htbp]
		\begin{center}
			\begin{circuitikz}[american currents]
				\draw

monhime

\documentclass[10pt, a3paper, titlepage, uplatex]{jsarticle} 
\usepackage{amsmath,ams symb}
\usepackage[dvipdfmx]{graphicx} 
\usepackage{circuitikz}

\begin{document}
	\begin{figure}[htbp]
		\begin{center}
		\begin{circuitikz}[american currents]
			\draw(0,0)