【转载】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));
   . k3 s# Y2 D* m' V+ U" T
4. d=c(find(rem(c,7)==2));
   
5. disp('在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数:')
   
6. disp(d)

复制代码

在1——500中，找出能同时满足用3除余2，用5除余3，用7除余2的所有整数:
2 \, ]$ ]/ k( C  23   128   233   338   443

晓夜寒风

matlab上标标注

1. subplot(1,1,1);
   7 X6 [0 F- f/ m# D# }: a4 j
2. x=0:0.1:2*pi;
   
3. y1=sin(x);
   " Y+ y6 g$ U9 N' W# D4 S
4. y2=exp(-x);
   
5. plot(x,y1,'--*',x,y2,':o');
   6 y, r7 Y  o  j5 m+ \
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)
   $ x$ e1 x; \4 C# Q0 \3 V
2. 1)x=[-10:1:10];
   ; ]9 ?. ?9 k6 y5 I5 H1 j7 |9 l; Z
3. for i=1:21
   
4. f(i)=(x(i)-1)*(x(i)>1)+(x(i)+1)*(x(i)<-1);
   1 ^! ]6 @, Y( s
5. end
   
6. plot(x,f,'r*')
   & E' l' j9 \; }
7. 2)x=[-10:10];
   
8. f=zeros(size(x));
   
9. a=x>1;
   
10. b=x<-1;
    
11. f(a)=x(a)-1;
    
12. f(b)=x(b)+1;
    3 h/ q& _% O& z# f
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
! X  g' [% }  r- z怎么样把大于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 a7 ]7 ?3 H2 k1 p  h/ O* A
2. a=rand;
   
3. k=find(x>a);
   
4. x(k(1))

复制代码

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

1. b=a(1:k);
   ' K# [3 c# |+ @* R1 L0 t$ N
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');
   + Z- P8 E8 w& S9 Y! ?7 v
2. if double(x1)>1   
   & \) R5 m8 U; z0 u9 w; S0 z5 `
3.   x=1
   ; }9 [, Y! F7 o
4. else   
   9 K; ^* |) o; L5 I
5.   x=0
   ! E9 r" h. j% ^* ^1 s
6. end

复制代码

晓夜寒风

Q5: MATLAB中用pie命令画饼图

1. function PlotPie
   
2. x=[1 2 3 4 5 6 7 8 9 10];
   " Q9 t7 t# d* }: R
3. y=[40 50 60 70 80 90 100 110 120 130 ];
   
4. str1={'A','B','C','D','E','F','G','H','I','J'};
   1 |8 R9 u. j6 N' z3 K
5. por=y./sum(y);
   
6. for i=1:length(por)   str{i}=[str1{i},' ',num2str(100*por(i)),'%'];
   
7. end
   ( I8 j6 t8 l5 ^
8. pie(x,y,str)

复制代码

晓夜寒风

Q6 非线性拟和y＝a*x^b+c的程序，求a,b,c.
  A8 m, w( [+ l) lfunction FitPaperData
) x* o3 H2 N% l0 Dformat long
  r$ a4 ~( P6 Z( K, Z6 W# J" \E1=rand(89,1);
E2=rand(89,1)+3;
F=inline('p(1).*x.^p(2)+p(3)','p','x');
6 F( I2 K9 B7 b! z! {9 X! T[p,r,j]=nlinfit(E1,E2,F,[1,2,3]);
P

晓夜寒风

Q7 已知向量A=[1 7 6 3 2 9 8];B=[1 6 3 8] ;如果我想把向量A中与向量B中相同的元素清除掉，例如最后得到C=[7 2 9]
应该怎么做
0 _# B- N& v: V( ]/ l' A2 ?6 {* H如果各个向量中均无重复元素，提供一个简便的方法
; ~; _5 C/ u) N. G7 ?+ z, W% nA=[1 7 6 3 2 9 8];
. p* D1 l! p, y1 c" R1 N$ U+ Y2 t' A/ @B=[1 6 3 8] ;
setdiff(A,B)

晓夜寒风

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

1. tic
   2 m- a0 K4 [4 X, u: a! f
2. a=rand(1,10);
   . K4 b( O- j4 K- o' o" d& Q
3. for i=1:length(a)
   
4. ff=['x.^2+',num2str(a(i))];
   / ]- f+ y; W# i, q; v4 H* |! P
5. f=inline(ff,'x');
   
6. d(i)=quadl(f,-1,1);
   
7. end
   + d' N) b8 d$ V8 K
8. d
   & _' g2 O: h/ k% C2 B" s
9. toc

复制代码

晓夜寒风

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

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

晓夜寒风

Q10 [1，2，3，NaN，5，NaN，7，8，9，10]
9 [/ L! k+ _, ]& d7 _* V/ a如何将NaN以外的8个数加起来，然后求这8个数的平均数。
3 T: r# K9 W3 ~+ ?0 y+ {-----------------------------------------------------------------
$ Z0 d: m0 r, [6 H6 P$ z
>> a=[1 2 3 nan 5 nan 7 8 9 10];
>> b=a;
>> b(find(isnan(a)))=[]

b =

1 2 3 5 7 8 9 10


. T# S$ o6 N. h8 B7 M5 @1 Ha=[1 2 3 nan 5 nan 7 8 9 10];
) T9 g% y  y/ k0 p$ S& Y) ~6 R& u简洁些可以这样:
mean_a=mean(a(~isnan(a)))