Yes. Let $(q_k,p_k)$ be a symplectic basis to $\mathbb{R}^{2n}$ so that $J(q_k)=p_k$ and $J(p_k)=-q_k$ for all $k$. Consider matrices $A_k$ and $B_{k,l}$ such that

- $A_k (q_k)=p_k \quad$and$\quad A_k (p_k) = -q_k$
- $A_k (q_l)=q_l \quad$and$\quad A_k (p_l) = p_l$ for $l \neq k$
- $B_{k,l}(q_k)=q_l, \quad B_{k,l}(q_l)=q_k, \quad B_{k,l}(p_k)=p_l,\quad B_{k,l}(p_l)=p_k\quad$ and $B_{k,l}$ is the identity away from the subspace spanned by the vectors $q_k,p_k,q_l,p_l$.

The matrices $A_k$ and $B_{k,l}$ are symplectic, as can easily be checked. If $E_k$ is the subspace spanned by the vectors $q_k,p_k$, we have $E_k=\text{Ker}(A_k^2+I)$. Also, note that $A_k^{T}=A_k^{-1}$ and $B_{k,l}^T=B_{k,l}$.

Now, since these matrices are symplectic, we can subsitute them into the formula for $N$ to get $A_k^{-1} N A_k = N$ and $B_{k,l} N B_{k,l}= N$. Taking the first equation and applying it to the vector $q_k$, we have $A_k^{-1}NA_k(q_k)=N(q_k)$, so $N(p_k)=A_kN(q_k)$. Similarly, $N(q_k)=-A_kN(p_k)$. In particular, $A_k^2 N(q_k)=-N(q_k)$ and $A_k^2N(p_k)=-N(p_k)$. Hence, $N(q_k), N(p_k) \in E_k$, and so there are numbers $a_k,b_k,c_k,d_k$ such that
$$ N(q_k)=a_k q_k+b_k p_k \\ N(p_k)=c_k q_k+d_k p_k .$$
Applying $A_k$ to either of these equations and using $N(q_k)=-AN(p_k)$ or $N(p_k)=AN(q_k)$, respectively, we have that $a_k=d_k$ and $c_k=-b_k$. Applying $B_{k,l}$ to these formulas for each $k,l$ shows that $a_k=a_l=:a$ and $b_k=b_l =: b$. So putting this all together, we get $N=aI+bJ$. Since $N$ is skew symmetric, it follows that $a=0$, and hence the result.