matlab总练习题Word文件下载.docx
- 文档编号:22255191
- 上传时间:2023-02-03
- 格式:DOCX
- 页数:25
- 大小:834.54KB
matlab总练习题Word文件下载.docx
《matlab总练习题Word文件下载.docx》由会员分享,可在线阅读,更多相关《matlab总练习题Word文件下载.docx(25页珍藏版)》请在冰豆网上搜索。
z=
alpha=x-z
alpha=
beta=y-z
beta=
alpha=[alpha0]
0
beta=[beta0]
cross(alpha,beta)
00
面积
a=11:
19;
b=a;
fork=1:
8
b=[b;
a+10*k];
end
rank(b)
2
a=vander(1:
9);
b=fliplr(a)
b=
Columns1through5
11111
124816
1392781
141664256
1525125625
16362161296
17493432401
18645124096
19817296561
Columns6through9
1111
3264128256
24372921876561
102440961638465536
31251562578125390625
7776466562799361679616
168071176498235435764801
327682621442097152
590495314414782969
det(b)
+15
方式一>
f=@(x,y)exp(x+y)+sin((x^2)+(y^2))
f=
@(x,y)exp(x+y)+sin((x^2)+(y^2))
f(1,2)
方式二
functionf=myfunfun(x,y)
f=exp(x+y)+sin((x^2)+(y^2));
myfunfun(1,2)
Char
a=ans;
sqrt2char(3-2)=a(3)
sqrt2char=
727
forx=1:
100
b(x)=str2num(sqrt2char(x))
sum(b)/100
f=@(x)(x^3)*sin(x)+(x^2)/3+x*cos(x)
@(x)(x^3)*sin(x)+(x^2)/3+x*cos(x)
ezplot(f,-2,1)
x0=fzero(f,-1)
x0=
另一根为0,是显然的
functiony=difun(x)
ifx<
-pi
y=-x-pi;
elseifx>
-pi&
x<
pi
y=sin(x);
else
y=(x-pi)/2;
y=[]
[]
forx=-6:
:
6
y=[ydifun(x)];
plot(x,y)
plot([-6:
6],y)
pi/4
%pi/4的理想值
矩形公式:
functiony=rectangle(n)
x=0:
1/n:
1;
a=1./(1+x.*x);
y=sum(a)*(1/n);
rectangle(1000)
rectangle(10000)
rectangle(100000)
梯形公式:
functiony=trapezoid(n)
begin=a
(1);
endd=a(n+1);
a
(1)=0;
a(n)=0;
y=sum(a)*(1/n)+begin*(1/n)*+endd*(1/n)*;
trapezoid(1000)
trapezoid(100)
Simpson公式
functiony=simpson(n)
a=thefun(x);
medium=[];
forx=1:
n
medium=[medium(a(x)+a(x+1))*];
y=begin*(1/n)*(1/6)+endd*(1/n)*(1/6)+sum(a)*(1/n)*(1/3)+sum(medium)*(1/n)*(1/6)*4;
functione=thefun(r)
e=1./(1+r.*r);
simpson(100)
simpson(10)
A=[621-1;
2410;
114-1;
-10-13];
b=[615-5]'
;
x=A\b
[diag(1:
4)eye(4)]
10001000
02000100
00300010
00040001
functionyh=yhsj(n)
yh=1;
disp
(1);
fork=2:
n
yh=[yh,0]+[0,yh];
disp(yh)
end
%n=11运行
x=sym('
x'
);
f=sqrt(1+((4/9)*x^(1/2))^2)
f=
((16*x)/81+1)^(1/2)
a=sym('
a'
b=sym('
b'
int(f,a,b)
(16*b+81)^(3/2)/216-(16*a+81)^(3/2)/216
t=sym('
t'
x=a*(t-sin(t))
a*(t-sin(t))
y=a*(1-cos(t))
-a*(cos(t)-1)
dx=diff(x)
dx=
dy=diff(y)
dy=
a*sin(t)
int(sqrt(dx^2+dy^2),0,2*pi)
8*(a^2)^(1/2)
p=polyfit([-pi-pi/20pi/2pi],[0-1010],5)
p=
Column6
plot([-pi:
pi/100:
pi],polyval(p,-pi:
pi));
holdon
pi],sin(-pi:
functiona=num2p(n)
strr=num2str(n);
index=size(strr);
indexx=index
(2);
a=[];
fork=indexx-1:
-1:
a=[arem(fix(n/(10^k)),10)];
num2p(494)
Columns1through8
95489298
Columns9through11
494
functionppi=ttry(n)
pointx=-1+2*rand(1,n);
pointy=-1+2*rand(1,n);
index=pointx.*pointx+pointy.*pointy;
indexx=index<
=1;
ppi=4*sum(indexx)/n;
ttry(100)
ttry(10000)
functiony=collatz(n)
ifn==1
y=1
return;
ifrem(n,2)==0
n=n/2;
else
n=3*n+1;
end
n
y=collatz(n);
collatz(12)
n=
6
3
10
5
16
8
4
1
f=@(x)1./(1+x.^2)
@(x)1./(1+x.^2)
y=f(x);
x=0:
p=polyfit(x,y,5)
intp=polyint(p)
intp=
Columns6through7
res1=polyval(intp,1)-polyval(intp,0)
res1=
%这是插值拟合解
atan
(1)
%actan
(1)公式解
fid=fopen('
C:
\陈民权的文档2\matlab作业\'
'
r'
fid=
b=fscanf(fid,'
%s'
fclose(fid)
b(1:
55)=[];
编辑“”文件:
y=1;
n=zeros(1,10)
flag=0;
1500000
ify==1000001
break
ifb(x)~='
['
&
flag==0&
isempty(str2num(b(x)))~=1
c=b(x);
d=str2double(c);
n(d+1)=n(d+1)+1;
y=y+1;
elseifb(x)=='
flag=1;
]'
flag=0;
givemerun%运行“givemerun”
0000000000
n=
Columns1through3
9995999758100026
Columns4through6
100229100230100359
Columns7through9
995489980099985
Column10
100106%分别为0,1,2,…,8,9在前一百万位出现的次
sum(n)
1000000
cc=num2str(n)
cc=
9995999758100026100229100230100359995489980099985100106>
fopen('
'
w'
fprintf(6,cc)
78
functionf=veryfun(n)
ifn==6174
return
index
(1)=(n-rem(n,1000))/1000;
index
(2)=fix(n/100)-index
(1)*10;
index(3)=rem(fix(n/10),10);
index(4)=rem(n,10);
index2=sort(index,'
ascend'
index3=sort(index,'
descend'
a=[1000100101];
max=index3*(a'
min=index2*(a'
max-min
veryfun(max-min);
veryfun(9864)
5175
5994
5355
1998
8082
8532
6174
functionm=tem(index,t)
m=index
(1).*exp(index
(2).*t)
[index,r,j,c,m,e]=nlinfit(t,m,@tem,[00])
index=
%在这里省略中间的迭代过程和表格变量输入过程
m=tem(index,[0:
25]);
plot([0:
25],m)
functionson=sonfind(A)
index=sum(A)/2;
fork=1:
5
b=nchoosek([1:
10],k);
sizee=size(b);
kk=sizee
(1);
kk
c=b(x,:
A_son=[];
fory=c
A_son=[A_sonA(y)];
ifsum(A_son)==index
son=A_son;
end
%集合元素有互异性,但不影响验证!
sonfind(11:
20)
sonfind([105510601010101010])
5560
sonfind([10101010901010101010])
90
%进一步验证(这里元素互异)
sonfind([21222324251112131465])
2122232425
sort([21222324251112131465])
11121314212223242565
sonfind(ans)
1112131465
201612
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- matlab 练习题