how to find mode shape and natural frequency from the below script.

f1=m1(:,1); r1=m1(:,2); i1=m1(:,3); mag1=sqrt(r1.^2+i1.^2); 
subplot(3, 3, 1);
plot(f1,mag1); 
f2=m2(:,1); r2=m2(:,2); i2=m2(:,3); abs2=sqrt(r2.^2+i2.^2); 
subplot(3, 3, 2);
plot(f2,abs2);
f3=m3(:,1); r3=m3(:,2); i3=m3(:,3); abs3=sqrt(r3.^2+i3.^2); 
subplot(3, 3, 3);
plot(f3,abs3);
f4=m4(:,1); r4=m4(:,2); i4=m4(:,3); abs4=sqrt(r4.^2+i4.^2); 
subplot(3, 3, 4);
plot(f4,abs4);
f5=m5(:,1); r5=m5(:,2); i5=m5(:,3); abs5=sqrt(r5.^2+i5.^2); 
subplot(3, 3, 5);
plot(f5,abs5);
f6=m6(:,1); r6=m6(:,2); i6=m6(:,3); abs6=sqrt(r6.^2+i6.^2); 
subplot(3, 3, 6);
plot(f6,abs6);
f7=m7(:,1); r7=m7(:,2); i7=m7(:,3); abs7=sqrt(r7.^2+i7.^2); 
subplot(3, 3, 7);
plot(f7,abs7);

[John D'Errico]

What do mode shape and natural frequency have to do with anything? We don't have your data, nor can see the figures that you plotted. So how can anyone else hazard a guess what to do from something that nobody but you can see?

[KSSV]

Mode shapes and natural frequency of what? what structure? What formulation you are using?