How do I transform this Digital filter to analog

How do I transform this digital filter to analog one, using bilinear transformation from z-plane to s-plane? I also have the information that s=(2*(z-1))/(ts*(z+1)) where ts is the period of discretization.

clear all;
close all;
clc;

fs=2600;
fn=fs/2;
ws=[710 839]/fn;
wp=[752 787]/fn;
rp=0.1;
rs=25;
k=1024;

[n,wn]=cheb2ord(wp,ws,rp,rs);
[nz,dz] = cheby2(n,rs,wn,'bandpass');
[H,w]=freqz(nz,dz,k/2);
m=abs(H);
ma=-20*log10(m);
plot(((w*fs)/(2*pi)),ma,'g');

Answer 1

[Star Strider]

Easiest:

[ns,ds] = cheby2(n,rs,wn,'bandpass','s');

Then use freqs instead of freqz.