diff --git a/run.m b/run.m
index 823edb0..089c2ba 100644
--- a/run.m
+++ b/run.m
@@ -2,7 +2,7 @@ clear
 clc
 % yalmip('clear')
 addpath('.\Powerflow')
-[~, ~, ~, ~,Volt,Vangle,Y,Yangle,r,c,newwordParameter,PG,QG,PD,QD,Balance]=pf('s1047.dat', '0');
+[~, ~, ~, ~,Volt,Vangle,Y,Yangle,r,c,newwordParameter,PG,QG,PD,QD,Balance]=pf('ieee301.dat', '0');
 %% 开始生成量测量
 sigma=0.03;% 标准差
 %% 电压
diff --git a/鍏紡/鍐呯偣娉曞叕寮�.tex b/鍏紡/鍐呯偣娉曞叕寮�.tex
new file mode 100644
index 0000000..dd2c1ee
--- /dev/null
+++ b/鍏紡/鍐呯偣娉曞叕寮�.tex
@@ -0,0 +1,104 @@
+\documentclass[10pt,a4paper,final]{report}
+\usepackage[utf8]{inputenc}
+\usepackage{amsmath}
+\usepackage{amsfonts}
+\usepackage{amssymb}
+\usepackage{fontspec}%浣跨敤xetex
+\setmainfont[BoldFont=榛戜綋]{瀹嬩綋}               % 浣跨敤绯荤粺榛樿瀛椾綋
+\XeTeXlinebreaklocale "zh"                      % 閽堝涓枃杩涜鏂
+\XeTeXlinebreakskip = 0pt plus 1pt minus 0.1pt  % 缁欎簣TeX鏂涓�瀹氳嚜鐢卞害
+\linespread{1.5}                                % 1.5鍊嶈璺�
+\begin{document}
+鐢变簬涔嬪墠鐨凣auss-Newton瀵�1047鑺傜偣鏀舵暃鎬т笉濂斤紙dX涓�0鐨勬椂鍊欙紝鏈�浼樻潯浠朵笉涓�0锛夛紝鍑嗗鏀圭敤鍐呯偣娉曟眰瑙ｃ��
+\begin{equation}
+f(x)=[z-h(x)]^T W [z-h(x)]
+\end{equation}
+鏈�浼樻潯浠剁瓑浠蜂负
+\begin{equation}
+\bigtriangledown f(x)=J^T W [z-h(x)]=0
+\end{equation}
+鍏朵腑$J$鏄�$h(x)$鐨�$Jacobi$鐭╅樀
+\newline
+瀵逛簬绾胯矾鍔熺巼
+绾胯矾鏈夊姛鍔熺巼Jacobi
+\begin{equation}
+\begin{aligned}
+\frac{\partial P_{ij}}{\partial V_1}=
+2V_1G_{ij}-V_2[cos(\theta_1 - \theta_2)G_{ij}+sin (\theta_1 - \theta_2)B_{ij}]
+\end{aligned}
+\end{equation}
+
+\begin{equation}
+\begin{aligned}
+\frac{\partial P_{12}}{\partial V_2}=
+-V_1[cos(\theta_1 - \theta_2)G_{ij}+sin (\theta_1 - \theta_2)B_{ij}]
+\end{aligned}
+\end{equation}
+
+
+鍒╃敤Newton娉曞涓婂紡杩涜灞曞紑锛屽啓鎴愬厓绱犲舰寮忋�備互涓や釜鍙橀噺涓轰緥
+\begin{equation}
+\bigtriangledown f= [ \bigtriangledown f_1\,, ... \,,\bigtriangledown f_n]^T
+\end{equation}
+\begin{equation}
+\begin{aligned}
+\bigtriangledown^2 f_1=(\omega_1 \frac{ \partial f_1 }{\partial x_1} \frac{ \partial f_1 }{\partial x_1}
++\omega_2 \frac{ \partial f_2 }{\partial x_1} \frac{ \partial f_2 }{\partial x_1} ) \Delta x_1 \\
++ (\omega_1 \frac{ \partial f_1 }{\partial x_1} \frac{ \partial f_1 }{\partial x_2}
++\omega_2 \frac{ \partial f_2 }{\partial x_1} \frac{ \partial f_2 }{\partial x_2} ) \Delta x_2 \\
++( \omega_1 f_1 \frac{ \partial ^2 f_1 }{ \partial x_1 \partial x_1 } + \omega_2 f_2 \frac{ \partial ^2 f_2 }{ \partial x_1 \partial x_1 } ) \Delta x_1 \\
++ ( \omega_1 f_1 \frac{ \partial ^2 f_1 }{ \partial x_1 \partial x_2 } + \omega_2 f_2 \frac{ \partial ^2 f_2 }{ \partial x_1 \partial x_2 } ) \Delta x_2
+\end{aligned}
+\end{equation}
+鍒╃敤matlab鐨勮鍙ュ彲浠ヨ〃绀轰负
+\begin{equation}
+J ^T W J + \tilde{Q}
+\end{equation}
+鍏朵腑
+\begin{equation}
+J ^T W J=
+(\omega_1 \frac{ \partial f_1 }{\partial x_1} \frac{ \partial f_1 }{\partial x_1}
++\omega_2 \frac{ \partial f_2 }{\partial x_1} \frac{ \partial f_2 }{\partial x_1} ) \Delta x_1 \\
++ (\omega_1 \frac{ \partial f_1 }{\partial x_1} \frac{ \partial f_1 }{\partial x_2}
++\omega_2 \frac{ \partial f_2 }{\partial x_1} \frac{ \partial f_2 }{\partial x_2} ) \Delta x_2
+\end{equation}
+
+\begin{equation}
+\label{闈炵煝閲忓寲浜岄樁瀵兼暟}
+\tilde{Q}=
+( \omega_1 f_1 \frac{ \partial ^2 f_1 }{ \partial x_1 \partial x_1 } + \omega_2 f_2 \frac{ \partial ^2 f_2 }{ \partial x_1 \partial x_1 } ) \Delta x_1 \\
++ ( \omega_1 f_1 \frac{ \partial ^2 f_1 }{ \partial x_1 \partial x_2 } + \omega_2 f_2 \frac{ \partial ^2 f_2 }{ \partial x_1 \partial x_2 } ) \Delta x_2
+\end{equation}
+寮�(\ref{闈炵煝閲忓寲浜岄樁瀵兼暟})鍙互鐢∕atlab璇█鍦ㄤ笉闇�瑕佸惊鐜殑鎯呭喌涓嬪鐞嗗ソ銆傚洜涓�
+% 绾胯矾鍔熺巼浜岄樁瀵兼暟
+\begin{equation}
+\begin{aligned}
+\frac{\partial^2 P_{12}}{\partial V_1^2}&=
+\frac{-2}{k^2}B_{12}\\
+&=\frac{-2B_{12}}{k^2}
+\end{aligned}
+\end{equation}
+
+\begin{equation}
+\begin{aligned}
+\frac{\partial^2 P_{12}}{\partial V_1 \partial V_2 }&=
+\frac{-1}{k}[sin(\theta_1 - \theta_2)G_{ij}-cos (\theta_1 - \theta_2)B_{ij}] \\
+&=
+\frac{-[sin(\theta_1 - \theta_2)G_{ij}-cos (\theta_1 - \theta_2)B_{ij}]}{k}
+\end{aligned}
+\end{equation}
+
+\begin{equation}
+\begin{aligned}
+\frac{\partial^2 P_{12}}{\partial V_2^2}&=0
+\end{aligned}
+\end{equation}
+鍒╃敤
+sparse([1 1],[1 1],[$ \omega_1 f_1 \frac{ \partial ^2 f_1 }{ \partial x_1 \partial x_1 } \quad \omega_2 f_2 \frac{ \partial ^2 f_2 }{ \partial x_1 \partial x_1 }  $] ) 鍙互琛ㄧず
+$
+\omega_1 f_1 \frac{ \partial ^2 f_1 }{ \partial x_1 \partial x_1 } + \omega_2 f_2 \frac{ \partial ^2 f_2 }{ \partial x_1 \partial x_1 }
+$
+鍦ㄦ捣妫煩闃典腑涓� $\Delta x_1$ 瀵瑰簲鐨勫厓绱狅紝涔熷氨鏄涓�琛岀涓�鍒楃殑鍏冪礌锛屽叾浠栫殑绫讳技銆傜敤杩欑鍔炴硶鍙互涓嶅埄鐢ㄧ煝閲忓寲鍏紡浠荤劧渚濋潬Matlab璇彞瀹炵幇鐭㈤噺鍖栬繍绠椼��
+
+\end{document}
+