【转载】matlab 学习点滴

晓夜寒风

在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数；

1. a=1:500;
   ' ~3 ~) y5 @: l1 v1 o. {
2. b=a(find(rem(a,3)==2));
   , i' G" u2 V7 L
3. c=b(find(rem(b,5)==3));
   7 {4 S  p. E' L4 T
4. d=c(find(rem(c,7)==2));
   
5. disp('在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数:')
   / X; U4 V9 z" I- s8 @( D5 q3 V
6. disp(d)

复制代码

在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数:
; e- J4 ?5 Q7 j; ?4 C3 e- r0 T  23   128   233   338   443

晓夜寒风

matlab上标标注

1. subplot(1,1,1);
   + O. N5 P+ K! k; q. l: X' _
2. x=0:0.1:2*pi;
   
3. y1=sin(x);
   * B& x5 q/ m3 @% Y5 t3 c* R5 k. b
4. y2=exp(-x);
   
5. plot(x,y1,'--*',x,y2,':o');
   
6. xlabel('x=0 to 2\pi');
   
7. ylabel('values of sin(x) and e^{-x}');
   
8. title('function plots fo sin(x) and e^{-x}');
   
9. legend('sin(x)','e^{-x}');

复制代码

晓夜寒风

请问在matlab里面怎么把一个三段的分段函数表示出来呢

1. f(x)=x-1 (x>1);f(x)=0 (abs(x)<1);f(x)=x+1 (x<-1)
   
2. 1)x=[-10:1:10];
   - r& N3 i- L. S. F
3. for i=1:21
   
4. f(i)=(x(i)-1)*(x(i)>1)+(x(i)+1)*(x(i)<-1);
   
5. end
   + F/ K3 [8 H" B
6. plot(x,f,'r*')
   
7. 2)x=[-10:10];
   
8. f=zeros(size(x));
   " j. d% M' g# q- F. ]- {- f
9. a=x>1;
   & }* j$ z5 U' j& D. }9 }
10. b=x<-1;
    " `# t3 f' Z% e4 p4 z. n5 U: C
11. f(a)=x(a)-1;
    
12. f(b)=x(b)+1;
    , R" q; B. v9 m. u0 T* u8 @# A# k
13. 3)f(x)=(x-1).*(x>1)+0.*(abs(x)<1)+(x+1).*(x<-1)

复制代码

晓夜寒风

Q1: 对于x = [0.1200 0.2400 0.3600 0.5800 0.6600 0.7400 0.8100 0.8700 0.9100 1.0000]' 随机产生一个数a = rand
怎么样把大于a的第一个x（i）选出来？

1. x = [0.1200 0.2400 0.3600 0.5800 0.6Q600 0.7400 0.8100 0.8700 0.9100 1.0000]'
   1 R, e8 V" S6 ~
2. a=rand;
   ; E: E* [/ m% C: y- Q" c
3. k=find(x>a);
   
4. x(k(1))

复制代码

Q2: 如何从列矩阵(n*1)中分段截取数据组成新的矩阵?比如用列矩阵a的前k个数据组成一个新的列矩阵b?中间k个呢?

1. b=a(1:k);
   ) c% [" U* w8 `. \6 a
2. b=a(l:l+k); (l+k<=n)
   

复制代码

Q3: [1，2，3，NaN，5，NaN，7，8，9，10]如何将NaN以外的8个数加起来，然后求这8个数的平均数。

1. mean_a=mean(a(~isnan(a)))

复制代码

Q4 matlab里面怎么将符号常数变为数值

1. syms x;x1=solve('x+1=2','x');
   
2. if double(x1)>1   
   
3.   x=1
   
4. else   
   + m! z% d; L+ a' G# ^
5.   x=0
   1 A) H% y& H' B4 w! [; ], O
6. end

复制代码

晓夜寒风

Q5: MATLAB中用pie命令画饼图

1. function PlotPie
   . R  A* a( Z$ t: A6 f
2. x=[1 2 3 4 5 6 7 8 9 10];
   
3. y=[40 50 60 70 80 90 100 110 120 130 ];
   
4. str1={'A','B','C','D','E','F','G','H','I','J'};
   ; G6 s4 X6 A. @3 N
5. por=y./sum(y);
   
6. for i=1:length(por)   str{i}=[str1{i},' ',num2str(100*por(i)),'%'];
   ( ^* _  o+ l! C4 O; Q# f
7. end
   
8. pie(x,y,str)

复制代码

晓夜寒风

Q6 非线性拟和y＝a*x^b+c的程序，求a,b,c.
function FitPaperData
format long
9 o% _% D, R6 W; EE1=rand(89,1);
E2=rand(89,1)+3;
( H& K2 o) e8 X) ~! c5 c& N6 _9 @F=inline('p(1).*x.^p(2)+p(3)','p','x');
3 |5 X% N5 n9 Y" ~[p,r,j]=nlinfit(E1,E2,F,[1,2,3]);
* l; H( _( t5 \/ FP

晓夜寒风

Q7 已知向量A=[1 7 6 3 2 9 8];B=[1 6 3 8] ;如果我想把向量A中与向量B中相同的元素清除掉，例如最后得到C=[7 2 9]
# A8 O( k* v  P0 U* z6 I; v$ c6 D应该怎么做
7 O8 V. L+ n4 H! t, \; f3 X: E2 w如果各个向量中均无重复元素，提供一个简便的方法
A=[1 7 6 3 2 9 8];
B=[1 6 3 8] ;
setdiff(A,B)

晓夜寒风

Q8 内联函数如何部分赋值传参
- P, v- I: q: e

1. tic
   5 N/ Z, r/ }8 Y* Q
2. a=rand(1,10);
   
3. for i=1:length(a)
   
4. ff=['x.^2+',num2str(a(i))];
   
5. f=inline(ff,'x');
   5 x3 K% F0 `3 o9 I  V% g
6. d(i)=quadl(f,-1,1);
   
7. end
   ( k" |0 X2 P% G. D7 f# }: z1 _
8. d
   . f! ^$ b6 Q- S0 L( O7 h1 D
9. toc

复制代码

晓夜寒风

Q13 我想写一个函数，输入参数是一个字符串比如是＂１６０＂，函数主要目的是完成数据存储，我现在的问题是存储文件时，文件命名的问题，比如我想把输入的参数作为文件名的一部分，比如希望输出的文件名为"temp１６０",这个该怎样实现啊，谢谢！！！i = 160;
$ z, E0 ~. d6 T( P! j3 [8 h! Vname = ['temp' , num2str(i), '.txt']; %假设你想把当前工作区的数据存储为temp１６０.txt；
eval('save(name)');
( @7 ?" |' P9 g! c2 A( ]  R
# u4 p, z' ]- {) _- {7 Y3 a( L具体的请看eval save fopen fread fwrite fprintf等函数的帮助

晓夜寒风

Q10 [1，2，3，NaN，5，NaN，7，8，9，10]
如何将NaN以外的8个数加起来，然后求这8个数的平均数。
-----------------------------------------------------------------

>> a=[1 2 3 nan 5 nan 7 8 9 10];
>> b=a;
4 |/ ]6 u9 ^. k* F  X5 Z>> b(find(isnan(a)))=[]
8 m5 I% [" b2 p8 ~1 ]9 s# n
b =

5 \5 l; T9 z+ Y4 G9 M1 2 3 5 7 8 9 10
% z) s! N! g2 X

- J. ]) G7 W; I8 I. w& d. Qa=[1 2 3 nan 5 nan 7 8 9 10];
简洁些可以这样:
mean_a=mean(a(~isnan(a)))