disp('newNums1='); disp(newNums1);
x2=find(mod(find(nums2),2)==1);%向量2 [values2]=nums2(x2); newNums2=values2; disp('nums2='); disp(nums2); disp('newNums2='); disp(newNums2);
x3=find(mod(find(nums3),2)==1);%向量3 [values3]=nums3(x3); newNums3=values3; disp('nums3='); disp(nums3); disp('newNums3='); disp(newNums3);
Matlab第四章
1. 有如下数据：
X y 1 1.00000 1.1 1.23368 1.2 1.55271 1.3 1.99372 1.4 2.61170
利用本章介绍的几种插值方法对其进行插值，得到每隔0.05的结果。
%4种插值方法
x=[1 1.1 1.2 1.3 1.4]; %x与y的值 y=[1.00000 1.23368 1.55271 1.99372 2.61170]; length_of_x=length(x);
scalar_x=[x(1):0.05:x(length_of_x)]; length_of_sx=length(scalar_x);
for i=1:length_of_sx 法
y_nearest(i)=interp1(x,y,scalar_x(i),'nearest'); y_liner(i)=interp1(x,y,scalar_x(i),'liner'); y_spline(i)=interp1(x,y,scalar_x(i),'spline'); y_cubic(i)=interp1(x,y,scalar_x(i),'cubic'); end
subplot(2,2,1),plot(x,y,'*'),hold
on,plot(scalar_x,y_nearest),title('nearest'); subplot(2,2,2),plot(x,y,'*'),hold
on,plot(scalar_x,y_liner),title('linear'); subplot(2,2,3),plot(x,y,'*'),hold
on,plot(scalar_x,y_spline),title('spline'); subplot(2,2,4),plot(x,y,'*'),hold
on,plot(scalar_x,y_cubic),title('cubic');
%定义4种方%输出
2. 2. 求下列函数的解，并绘制图形。 (1) y=ex-x5，初始点为x=8. (2) y=xsin(x)
>> y=@(x)exp(x)-x^5; >> x=fzero(y,8) x =
12.7132
>> fplot(y,[10,14]); >> hold on
>> plot(x,y(x),'r*') >>y=@(x)x*sin(x); >> x=fzero(y,0) x = 0
3.求下列函数的极值。 (1) z=x2-(y-1)2 (2) z=(x-y+1)2
>> z=@(x)x(1)^2-(x(2)-1)^2;
>> [x,fvalue]=fminsearch(z,[-1,1])
Exiting: Maximum number of function evaluations has been exceeded - increase MaxFunEvals option. Current function value:
-29278786696692084000000000000000000000000000000000000000000000000000000000000000000000.000000 x =
1.0e+043 *
1.1138 1.2383 fvalue =
-2.9279e+085
>> z = @(x)(x(1)-x(2)+1)^2; >> [x,fvalue]=fminsearch(z,[-1,1]) x =
-1.2073 -0.2073 fvalue =
1.9045e-022
4.计算下列积分。 (1)
(2)
>> f = @(x)x+x.^3+x.^5; >> q = quad(f,-1,1) q =
2.2204e-016
>> f= @(x,y)sin(y).*(x+y)./(x.^2+4); >> q = dblquad(f,1,10,1,10) q =
5.5254
Matlab第五章
1.编制一个脚本，查找给定字符串中指定字符出现的次数和位置.
%查找给定字符串中指定字符出现的次数和位置
letter='a'; %指定字符 string='I have a long string'; %给定字符串 places=findstr(string,letter); %位置 times=length(places); %次数 disp('places are :'); disp(places); disp('times are :'); disp(times);
2.编写一个脚本，判断输入字符串中每个单词的首字母是否为大写，若不是则将其修改为大写，其他字母为小写.
%将每个单词首字母变为大写
str = 'I have a long string'; disp(str);
nlength = length(str); for k=1:nlength
if (k==1 || isspace(str(k-1))) && (str(k)<='z' && str(k)>='a') str(k) = char(double(str(k)) - 32); end end disp(str);
3.创建 2×2 单元数组，第 1、2 个元素为字符串，第三个元素为整型变量，第四个元素为双精度（double）类型，并将其用图形表示.
>> A=cell(2,2) A =
[] [] [] []
>> A(1,1)={'string1'}; >> A(1,2)={'string2'}; >> A(2,1)={uint8(100)}; >> A(2,2)={double(60)};