青岛大学MATLAB6、7章

第六章、第七章6-2
C=0.25*10^(-9);
L=0.1*10^(-6);
R=100;
f=1*10^5:10:1*10^7;
I=120./(R+j*2*pi*f*L-j./(2*pi*f*C)); magyI=abs(I);
angI=angle(I);
figure(1);
subplot(211)
plot(f,magyI);
title('幅频特性');
grid on;
subplot(212)
semilogx(f,magyI);
title('幅频特性（对数坐标）'); grid on;
figure(2);
subplot(211)
plot(f,angI);
title('相频特性');
012345678910
x 106
幅频特性
10
1010
幅频特性（对数坐标）
10
1010
幅频特性（对数坐标）
10
1010
相频特性（对数坐标）
grid on;
subplot(212)
semilogx(f,angI);
title('相频特性（对数坐标）');
grid on;
figure(3);
subplot(211)
semilogx(f,magyI);
title('幅频特性（对数坐标）'); grid on;
subplot(212)
semilogx(f,angI);
title('相频特性（对数坐标）'); grid on;
6-10
l=0:0.01:2*pi;
y=cos(l)+j*sin(l);
subplot(211)
plot(y);
x 106
相频特性
10
1010
相频特性（对数坐标）
realy=real(y); imagy=imag(y); subplot(212)
plot3(realy,imagy,l)
6-11
[x,y]=meshgrid(-1:.1:1,-2*pi:0.2:2*pi); z=exp(x+i*y);
contour (x,y,z);
[x,y]=meshgrid(-1:.1:1,-2*pi:0.2:2*pi); z=exp(x+i*y);
surf (x,y,z);
[x,y]=meshgrid(-1:.1:1,-2*pi:0.2:2*pi); z=exp(x+i*y);
mesh (x,y,z);
-1-0.500.51 -1
-0.5
0.5
1
1 1
5
10
s=input('Please input a string:');
l=size(s);
n=0;
for i=1:l(:,2)
if s(i)=='s'
n=n+1;
else
end
end
n
验证：
Please input a string:'SHE is busy.'
n =2
6-17
function str2=caps(str1)
str1(1)=upper(str1(1));
s=size(str1);
for ii=1:(s(:,2)-1)
if str1(ii)==' '
str1(ii+1)=upper(str1(ii+1));
else
str1(ii+1)=lower(str1(ii+1));
end
end
str1(end)=lower(str1(end));
str2=str1;
验证：
caps('I am a girl.')
ans =I Am A Girl.
7-7
function y=dist3(point1,point2)
msg=nargchk(2,2,2);
error(msg);
y=sqrt((point1.x-point2.x)^2+(point1.y-point2.y)^2); 验证：
point2=struct('x','2','y','0');
>> point1=struct('x','0','y','0');
>> dist3(point1,point2)
7-15
='john doe';
student.addr1='123 Main Street';
student.city='Anytown';
student.state='LA';
student.zip='71211';
student
student = name: 'john doe'
addr1: '123 Main Street'
city: 'Anytown'
state: 'LA'
zip: '71211'
7-16
tspan=0:0.1:10;
v0=0;
[t,v]=ode45(@fun,tspan,v0);
plot(t,v);
function dv=fun(t,v)
R=1*10^6;
C=1*10^-6;
if t>0&t<6
vin=(-1)^t;
else
vin=0;
end
dv=(-v+vin)/(R*C);
7-17
tspan=0:0.1:10;
v0=0;
[t,v]=ode45(@fun1,tspan,v0); plot(t,v);
function dy=fun1(t,y)
if t>=0&t<=5
a=t;
else
a=0;
end
dy=-y+a;。