【转载】matlab 学习点滴

晓夜寒风

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

1. a=1:500;
   3 F( q5 Q/ @% A: g, ^
2. b=a(find(rem(a,3)==2));
   
3. c=b(find(rem(b,5)==3));
   & }6 h6 B5 I7 b1 s! s
4. d=c(find(rem(c,7)==2));
   ) W! B( H" o, W3 C, U7 x+ ~( h% s( k
5. disp('在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数:')
   & I* U3 E$ g! O: L
6. disp(d)

复制代码

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

晓夜寒风

matlab上标标注

1. subplot(1,1,1);
   
2. x=0:0.1:2*pi;
   
3. y1=sin(x);
   8 G. H. C) U7 P' f0 W
4. y2=exp(-x);
   0 {% G% Y2 W  _2 l; @) T
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];
   7 L; ~- z: E! N# E* d* s% d
3. for i=1:21
   
4. f(i)=(x(i)-1)*(x(i)>1)+(x(i)+1)*(x(i)<-1);
   
5. end
   
6. plot(x,f,'r*')
   
7. 2)x=[-10:10];
   
8. f=zeros(size(x));
   1 ^& G) }9 N! \6 b' m4 ^
9. a=x>1;
   " C, S2 [) [0 o  W- R& B) h
10. b=x<-1;
    # U$ T2 L/ J4 R- c( I7 n" W
11. f(a)=x(a)-1;
    3 }9 [1 ]; H6 h* R7 K4 F/ _
12. f(b)=x(b)+1;
    * p! `5 Z0 t$ b8 p* }, ?, X
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]'
   ( O4 B2 {( d% _4 W
2. a=rand;
   
3. k=find(x>a);
   
4. x(k(1))

复制代码

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

1. b=a(1:k);
   : h2 k% d# u! g; ]
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');
   ! _3 i7 ]* m1 x1 ?( z9 P
2. if double(x1)>1   
   5 X, b7 s8 m% |- l1 r
3.   x=1
   
4. else   
   
5.   x=0
   
6. end

复制代码

晓夜寒风

Q5: MATLAB中用pie命令画饼图

1. function PlotPie
   : @1 ^5 C* B  b. }  B
2. x=[1 2 3 4 5 6 7 8 9 10];
   
3. y=[40 50 60 70 80 90 100 110 120 130 ];
   ; F! ~% x; m9 v& y& ^6 {* m1 {; J
4. str1={'A','B','C','D','E','F','G','H','I','J'};
   
5. por=y./sum(y);
   
6. for i=1:length(por)   str{i}=[str1{i},' ',num2str(100*por(i)),'%'];
   & `- w5 I, o) e+ e! h
7. end
   . v* t( V; D" @" p2 c# `3 E1 \
8. pie(x,y,str)

复制代码

晓夜寒风

Q6 非线性拟和y＝a*x^b+c的程序，求a,b,c.
function FitPaperData
5 X. z& d* x& r& `& O4 q( xformat long
# Y8 q) u  {& t0 e- nE1=rand(89,1);
E2=rand(89,1)+3;
1 P" v4 m. O# @% K' S- ^  }F=inline('p(1).*x.^p(2)+p(3)','p','x');
[p,r,j]=nlinfit(E1,E2,F,[1,2,3]);
: W' U# k" R& jP

晓夜寒风

Q7 已知向量A=[1 7 6 3 2 9 8];B=[1 6 3 8] ;如果我想把向量A中与向量B中相同的元素清除掉，例如最后得到C=[7 2 9]
应该怎么做
+ W7 E7 S  @! V如果各个向量中均无重复元素，提供一个简便的方法
6 |, h) Q1 n& j6 kA=[1 7 6 3 2 9 8];
B=[1 6 3 8] ;
setdiff(A,B)

晓夜寒风

Q8 内联函数如何部分赋值传参
' U. u7 t4 y, h- ?; Y

1. tic
   
2. a=rand(1,10);
   
3. for i=1:length(a)
   
4. ff=['x.^2+',num2str(a(i))];
   
5. f=inline(ff,'x');
   + J; C# u: L/ ^& h8 j
6. d(i)=quadl(f,-1,1);
   
7. end
   ; J' |2 {8 {1 G5 c: I
8. d
   
9. toc

复制代码

晓夜寒风

Q13 我想写一个函数，输入参数是一个字符串比如是＂１６０＂，函数主要目的是完成数据存储，我现在的问题是存储文件时，文件命名的问题，比如我想把输入的参数作为文件名的一部分，比如希望输出的文件名为"temp１６０",这个该怎样实现啊，谢谢！！！i = 160;
name = ['temp' , num2str(i), '.txt']; %假设你想把当前工作区的数据存储为temp１６０.txt；
5 @5 j+ Q9 {6 v, keval('save(name)');

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

晓夜寒风

Q10 [1，2，3，NaN，5，NaN，7，8，9，10]
如何将NaN以外的8个数加起来，然后求这8个数的平均数。
6 Z$ _6 e; e' {-----------------------------------------------------------------
& @- W3 ]/ G; r8 E. s, P; D
8 U6 p9 a0 J* l3 I! N* k% [3 \>> a=[1 2 3 nan 5 nan 7 8 9 10];
! L9 S- L' J5 w: l& h>> b=a;
5 }) x8 l4 J/ ^) f>> b(find(isnan(a)))=[]
% A7 I$ ^) h  }: x& X0 e
b =
1 X  a- `6 S0 ^1 t' D+ \& @. B
% E' r6 ]8 d, R( P/ N1 2 3 5 7 8 9 10
! m  ~4 p  \6 I" K

4 W1 @& N7 ^  l: sa=[1 2 3 nan 5 nan 7 8 9 10];
简洁些可以这样:
3 K- u% i7 W, qmean_a=mean(a(~isnan(a)))