数模实例一电阻问题Word文档下载推荐.docx
《数模实例一电阻问题Word文档下载推荐.docx》由会员分享，可在线阅读，更多相关《数模实例一电阻问题Word文档下载推荐.docx（12页珍藏版）》请在冰点文库上搜索。
root（ans）
root（A）
helproot
polyval（p1,20）
roots（p1）
Dp=polyder（p1）
theta=0:
0.01:
6*pi;
rho=5*sin（4*theta/3）;
polar（theta,rho）figure;
polar（theta,rho）figure;
polar（theta,rho）
%--2012/5/1310:
11--%
figure;
rho=5*sin（theta/3）;
polyvalm（p1,20）
eig（compan（p1））
diff（x^3）
diff（x^4-（69*x^3）/10-（3863*x^2）/50-（8613*x）/100+12091/20）
symx
symsx
%--2012/6/1011:
49--%
A=[11121314;
21222324;
31323334;
41424344];
rank（A）
[l,u]=lu（A）
x=[-pi:
0.05:
pi];
y=sin（tan（x））-tan（sin（x））;
plot（x,y）
ezplot（'
x^2*sin（x+y^2）+y^2*exp（x+y）+5*cos（x^2+y）'
）
[X,Y]=meshgrid（-10:
0.5:
10）;
a=sqrt（X.^2+Y.^2）+eps;
Z=sin（a）./a;
mesh（X,Y,Z）;
symsxab;
f=x*（1+a/x）^x*sin（b/x）;
L=limit（f,x,inf）
symst;
f=t^2*exp（-2*t）*sin（t+pi）;
F=laplace（f）
symsabcx;
solve（'
a*x^2+b*x+c'
symx;
cos（2*x）+sin（x）'
A=[123;
-137;
903];
b=[147]'
;
x=A\b
x2=inv（A）*b
A=[123];
B=[3,4,5];
C=cross（A,B）
S=dot（A,cross（B,C））
symsx;
int（（x+sin（x））/（1+cos（x））,x,0,pi/2）
int（cos（x）^5*sin（x）,x,0,pi/2）
a=limit（x^2/（sin（x/3））^2,0）
b=limit（（tan（x）-sin（x））/（sin（x））^3,0）
d=[7408174405743737469174256739987400674099];
m=mean（d）
v=var（d）
xdwc=sqrt（v）/m
minZ=-6X1-4X2
a=[2,3;
4,2];
c=[-6,-4];
b=[100,120];
vlb=[0,0];
vub=[];
[x,lam]=lp（c,a,b,vlb,vub）
[x,lam]=lip（c,a,b,vlb,vub）
vub=[1010]
vub=[];
a=[330335340345350355360365370375385390405425];
b=[19792182226521992221226223022307236624292477261928433899];
plot（a,b,a,b,'
xlabel（'
电压（V）'
ylabel（'
计数'
title（'
坪特性曲线'
function[x,y]=jitu（t,j）
fori=1:
t
x=i;
y=t-x;
if2*x+4*y==j
break;
end
[x,y]=jitu（t,j）
symsxytj;
A={'
张三'
'
李四'
王五'
};
B={201001040081,201001040082,201001040083};
C={'
男'
女'
D={19840111,19850215,19870623};
E={'
语文'
数学'
英语'
物理'
F={91929394,89888786,79787776};
MM={A,B,C,D,E,F}
A.a1='
A.a2='
A.a3='
C.c1='
%--2012/7/1115:
vap（pi,300）
vap（pi,3）
symsvar_listvar_props
vap（pi,30）
vpa（pi,300）
A=[1,2,3;
45,6;
7,80]
7,80];
A=[[A;
[123]],[1;
2;
3;
4]];
B=[1+9i,2+8i,3+7j;
4+6j,5+5i,6+4i;
7+3i,8+2j1i]
v1=0:
0.2:
v2=0:
-0.1:
pi,v3=0:
pi,v4=pi:
-1:
B1=A（1:
2:
end,:
B2=A（[3,2,1],[1,1,1]）
B3=A（:
end;
1）
end;
end:
4,5,6;
7,8,0];
B=A.^A
C=A^（1/3）,e=norm（A-C^3）
C=A^（1/3）
e=norm（A-C^3）
C=A.^（1/3）
e=norm（A-C.^3）
j1=exp（sqrt（-1）*2*pi/3）;
A1=C*j1,A2=C*j1^2,norm（A-A1^3）,norm（A-A1^3）
find（A>
=5）'
[i,j]=find（A>
=5）
a1=all（A>
=5）,a2=any（A>
all（A（:
）>
symss;
P=（s+3）^2*（S^2+3*s+2）*（s^3+12*s^2+48*s+64）
P=（s+3）^2*（s^2+3*s+2）*（s^3+12*s^2+48*s+64）
P1=simple（P）
[P2,m]=simple（P）
P3=expand（P）
symssz;
P=（s+3）^2*（s^2+3*s+2）*（s^3+12*s^2+48*s+64）;
P1=simple（subs（P,s,（z-1）/（z+1）））,latex（P1）
helplatex
LATEX（P1）
P1=simple（subs（P,s,（z-1）/（z+1）））,latex（P1）
P1=simple（subs（P,s,（z-1）/（z+1）））,LaTeX（P1）
%--2012/7/128:
29--%
%--2012/7/1215:
31--%
A=[-0.2765,05772,1.4597,2.1091,1.191,-1.6187];
v1=floor（A）,v2=ceil（A）,v3=round（A）,v4=fix（）
v4=fix（A）
A=[-0.2765,0.5772,1.4597,2.1091,1.191,-1.6187];
v1=floor（A）,v2=ceil（A）,v3=round（A）,v4=fix（A）
A=hilb（3）
[n,d]=rat（A）
helprat
m=sym（1856120）;
n=sym（1483720）;
gcd（m,n）,lcm（m,n）,factor（lcm（n,m））
A=1:
1000;
B=A（isprime（A））
s1=0;
100,s1=s1+i;
end,s1
s2=0;
i=1;
while（i<
=100）,s2=s2+i;
i=i+1;
end,s2
s=0;
m=0;
while（s<
=10000）,m=m+1;
s=s+m;
end,s,m
tic,s=0;
100000,s=s+1/2^i+1/3^i;
toc
tic,i=1:
100000;
s=sum（1./2.^i+1./3.^i）;
10000,s=s+i;
ifs>
10000,break;
end,end
[m1,s1]=findsum（145323）
helpmyhilb
A1=myhilb（3,4）
A2=myhilb（4）
my_fact（11）
tic,my_fibo（25）,toc
tic,a=[1,1];
fork=3:
100,a（k）=a（k-1）+a（k+2）;
end,toc
P=[12405];
Q=[12];
F=[123];
D=convs（P,Q,F）
E=convs（convs（P,Q）,F）
G=convs（P,Q,F,[1,1],[1,3],[1,1]）
-1.8,-1.799:
.001:
-1.2,-1.2:
1.2,...1.201:
0.001:
1.8,1.81:
y=sin（tan（x））-tan（sin（x））;
s
ds
45
34
4e
1.2,1.201:
%--2012/7/1219:
23--%
x=[-2:
0.02:
2];
y=1.1*sign（x）.*（abs（x）>
1.1）+x.*（abs（x）<
=1.1）;
plot（[-2,-1.1,1.1,1.2],[-1.1,-1.1,1.1,1.1]）
t=0:
.2:
2*pi;
y=sin（t）;
subplot（2,2,1）,stairs（t,y）
subplot（2,2,2）,stem（t,y）
subplot（2,2,3）,bar（t,y）
subplot（2,2,4）,semilogx（t,y）
[-10,10]）
%--2012/7/1316:
18--%
t:
.12:
x=t.^3.*sin（3*t）.*exp（-t）;
y=t.^3.*cos（3*t）.*exp（-t）;
z=t.^2;
plot3（x,y,z）,grid
stem3（x,y,z）;
[x,y]=meshgrid（-3:
3,-2:
2）;
z=（x.^2-2*x）.*exp（-x.^2-y.^2-x.*y）;
mesh（x,y,z）
surf（x,y,z）
shadingflat
shadinginterp
waterfall（x,y,z）
contour3（x,y,z）
[x,y]=meshgrd（-2:
.1:
[x,y]=meshgrid（-2:
z=1./（sqrt（（1-x）.^2+y.^2））+1./（sqrt（（1+x）.^2+y.^2））;
surf（x,y,z）,shadingflat
xx=[-2:
-1.2,-1.1:
-0.9,-0.8:
0.8,0.9:
1.1,1.2:
yy=[-1:
-0.2,-0.1:
0.1,0.2:
1];
[x,y]=meshgrid（xx,yy）;
zlim（[0,15]）
[x,y]=meshgrid（-1.5:
1.5,-2:
z=0.5457*exp（-0.75*y.^2-3.75*x.^2-1.5*x）.*（x+y>
1）0.7575*exp（-y.^2-6*x.^2）.*（（x+y>
-1）&
（x+y<
=1））+0.5457*exp（-0.75*y.^2-3.75*x.^2+1.5*x）.*（x+y<
=-1）;
1）+0.7575*exp（-y.^2-6*x.^2）.*（（x+y>
surf（x,y,z）,xlim（[-1.51.5]）;
view（80,10）,xlim（[-1.51.5]）
1.2）;
subplot（224）,surf（x,y,z）
subplot（221）,surf（x,y,z）,view（0,90）
subplot（222）,surf（x,y,z）,view（90,0）
subplot（223）,surf（x,y,z）,view（0,0）
W=imread（'
baibing-02.jpg'
imtool（W）
W1=rgb2gray（W）;
W2=edge（W1）
figure
imtool（~W2）
f=x*（1+a/x）^x*sin（b/x）;
%--2012/7/1321:
37--%
limit（（exp（x^3）-1）/（1-cos（sqrt（x-sin（x））））,x,0,'
right'
x=-0.1:
0.1;
y=（exp（x.^3）-1）./（1-cos（sqrt（x-sin（x））））;
plot（x,y,'
-'
[0],[12],'
o'
limit（（exp（x^3）-1）/（1-cos（sqrt（x-sin（x））））,x,0）
f=tan（t）;
L1=limit（f,t,pi/2,'
left'
）,L2=limit（f,t,pi/2,'
symsxya;
f=exp（-1/（y^2+x^2））*sin（x）^2/x^2*（1+1/y^2）^（x+a^2*y^2）;
L=limit（limit（f,x,1/sqrt（y））,y,inf）
f=sin（x）/（x^2+4*x+3）;
f1=diff（f）
x1=0:
5;
y=subs（f,x,x1）;
y1=subs（f1,x,x1）;
ploy（x1,y,x1,y1,'
:
'
plot（x1,y,x1,y1,'
f4=diff（f,x,4）
latex（f4）
tic,diff（f,x,100）;
y=sym（'
f（t）'
G=simple（diff（t^2*sin（t）*y,t,3））
simple（subs（G,y,exp（-t）））,simple（diff（t^2*sin（t）*exp（-t）,3）-ans）
%--2012/7/1415:
38--%
%--2012/7/1421:
42--%
F=[（290,290,290）（85,90,100）（85,90,100）（75,80,100）（75,80,95）（75,80,85）;
（280,280,280）（50,60,65）（75,80,85）（85,90,100）（65,70,75）;
（278,278,278）（50,60,65）（65,70,75）（75,80,85）（85,95,100）;
（277,277,277）（85,90,100）（75,80,85）（65,70,75）（85,90,100）;
（275,275,275）（75,80,85）（65,70,75）（50,60,65）（85,90,100）;
（275,275,275）（50,60,65）（75,80,85）（85,90,100）（75,80,85）;
（274,274,274）（85,90,100）（75,80,85）（65,70,75）（75,80,85）;
（273,273,273）（75,80,85）（85,90,100）（75,80,85）（65,75,75）]
%--2012/7/1518:
39--%
a=polyfit（t,r,1）
升级会员 