【转载】matlab 学习点滴

晓夜寒风

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

1. a=1:500;
   
2. b=a(find(rem(a,3)==2));
   
3. c=b(find(rem(b,5)==3));
   , d; M" G! l' ]8 ?! D  G
4. d=c(find(rem(c,7)==2));
   
5. disp('在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数:')
   6 I% J" X7 E: i5 n' V& A9 i2 F4 }9 ~
6. disp(d)

复制代码

在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数:
  23   128   233   338   443

晓夜寒风

matlab上标标注

1. subplot(1,1,1);
   . `: F+ [7 `# }' W
2. x=0:0.1:2*pi;
   
3. y1=sin(x);
   
4. y2=exp(-x);
   9 k5 d5 F  Y& @' F# X
5. plot(x,y1,'--*',x,y2,':o');
     `5 q# \0 L" b* k& R
6. xlabel('x=0 to 2\pi');
   : k4 K6 b7 c  p4 }2 E3 t5 ], g) H
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)
   : Y$ j; \& g& I: q7 f
2. 1)x=[-10:1:10];
   1 ~& L- l/ c6 u5 z5 U' p; ?
3. for i=1:21
   
4. f(i)=(x(i)-1)*(x(i)>1)+(x(i)+1)*(x(i)<-1);
   ! j( h' e' B2 V' n) t' [
5. end
   
6. plot(x,f,'r*')
   
7. 2)x=[-10:10];
   % A& E( q" {6 x7 z) e. N
8. f=zeros(size(x));
   
9. a=x>1;
   
10. b=x<-1;
    
11. f(a)=x(a)-1;
    ! i( Z0 n7 A) V
12. f(b)=x(b)+1;
    
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]'
   : i: ?- [4 z* o0 W
2. a=rand;
   
3. k=find(x>a);
   ' t+ B0 I1 c5 l$ `8 L
4. x(k(1))

复制代码

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

1. b=a(1:k);
   ) b# g: K/ y, b+ F
2. b=a(l:l+k); (l+k<=n)
   " i. m5 U# X% l1 _8 y) H

复制代码

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');
     @, ?9 U# E- m) v4 v) w
2. if double(x1)>1   
   
3.   x=1
   
4. else   
   ( U+ `, i4 U* c" r1 E
5.   x=0
   
6. end

复制代码

晓夜寒风

Q5: MATLAB中用pie命令画饼图

1. function PlotPie
   
2. x=[1 2 3 4 5 6 7 8 9 10];
   + d8 u. w# i" X' `  U) G
3. y=[40 50 60 70 80 90 100 110 120 130 ];
   9 W: H' ?- N7 l0 b/ S( V% b
4. str1={'A','B','C','D','E','F','G','H','I','J'};
   9 x  U8 x7 b3 ~5 A
5. por=y./sum(y);
   , d' r7 ?* X9 o0 @$ }( G
6. for i=1:length(por)   str{i}=[str1{i},' ',num2str(100*por(i)),'%'];
   7 \3 e0 [  B5 e; {) [" F
7. end
   
8. pie(x,y,str)

复制代码

晓夜寒风

Q6 非线性拟和y＝a*x^b+c的程序，求a,b,c.
  [: r# O/ v0 w$ K/ K0 tfunction FitPaperData
9 A% ]# v$ Z4 ]% nformat long
E1=rand(89,1);
E2=rand(89,1)+3;
F=inline('p(1).*x.^p(2)+p(3)','p','x');
7 r0 O* N9 q' g& H/ I- p[p,r,j]=nlinfit(E1,E2,F,[1,2,3]);
6 w3 S! ~1 D  [5 L+ b* wP

晓夜寒风

Q7 已知向量A=[1 7 6 3 2 9 8];B=[1 6 3 8] ;如果我想把向量A中与向量B中相同的元素清除掉，例如最后得到C=[7 2 9]
& ^* n/ n6 Z7 s2 w2 r  \8 E8 d  y应该怎么做
如果各个向量中均无重复元素，提供一个简便的方法
A=[1 7 6 3 2 9 8];
* k0 z9 f3 G7 g% {" nB=[1 6 3 8] ;
setdiff(A,B)

晓夜寒风

Q8 内联函数如何部分赋值传参

$ }' N; @& \9 e

1. tic
     X' g9 s* q1 x. Y2 e
2. a=rand(1,10);
   
3. for i=1:length(a)
   
4. ff=['x.^2+',num2str(a(i))];
   
5. f=inline(ff,'x');
   
6. d(i)=quadl(f,-1,1);
   # q1 x8 G  P2 m6 A
7. end
   
8. d
   
9. toc

复制代码

晓夜寒风

Q13 我想写一个函数，输入参数是一个字符串比如是＂１６０＂，函数主要目的是完成数据存储，我现在的问题是存储文件时，文件命名的问题，比如我想把输入的参数作为文件名的一部分，比如希望输出的文件名为"temp１６０",这个该怎样实现啊，谢谢！！！i = 160;
4 q) s1 V* l& R: f3 j8 F% _3 nname = ['temp' , num2str(i), '.txt']; %假设你想把当前工作区的数据存储为temp１６０.txt；
eval('save(name)');

具体的请看eval save fopen fread fwrite fprintf等函数的帮助

晓夜寒风

Q10 [1，2，3，NaN，5，NaN，7，8，9，10]
7 Q, |8 O* j% a- b' T: \' ?* Q如何将NaN以外的8个数加起来，然后求这8个数的平均数。
9 }9 j7 m( ]% A8 s-----------------------------------------------------------------
' j' u- l; K1 t# u% X
>> a=[1 2 3 nan 5 nan 7 8 9 10];
>> b=a;
3 P0 |" k5 l& u9 B& Y" r; {>> b(find(isnan(a)))=[]
* n& `. r& N" u- `
b =
9 v  G! o" y* O, T6 {$ l. M
' s& t7 V' V! b' j; E1 2 3 5 7 8 9 10

# t; j# x/ e# c+ U
& l# T, ?* i$ ^1 n) _. ja=[1 2 3 nan 5 nan 7 8 9 10];
! ~& x( W' u2 Z, v7 M简洁些可以这样:
1 U9 Q2 z# i( j; c  D  S( dmean_a=mean(a(~isnan(a)))